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ABSTRACT 

A perturbative renormalization group is formulated for the study of Hamiltonian light- 
front field theory near a critical Gaussian fixed point. The only light- front renormalization 
group transformations found here that can be approximated by dropping irrelevant oper- 
ators and using perturbation theory near Gaussian fixed volumes, employ invariant- mass 
cutoffs. These cutoffs violate covariance and cluster decomposition, and allow functions of 
longitudinal momenta to appear in all relevant, marginal, and irrelevant operators. These 
functions can be determined by insisting that the Hamiltonian display a coupling constant 
coherence, with the number of couplings that explicitly run with the cutoff scale being lim- 
ited and all other couplings depending on this scale only through their dependence on the 
running couplings. Examples are given that show how coupling coherence restores Lorentz 
covariance and cluster decomposition, as recently speculated by Wilson and the author. 
The ultimate goal of this work is a practical Lorentz metric version of the renormalization 
group, and practical renormalization techniques for light-front quantum chromodynamics. 
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I. Introduction 

In a series of remarkable papers Wilson reformulated the original renormalization group 
approach to relativistic field theory [1-3], initially developing the modern renormalization 
group as a tool for the study of the strong interaction in Minkowski space [4-6]. He was 
later diverted to the study of Euclidean field theory [7] and statistical field theory [8- 
11], where it was possible to implement perturbative and numerical renormalization group 
transformations for theories of physical interest. He has written a number of reviews of 
this work [12-14] and two simple introductions [15,16]. This paper relies heavily on Ref. 
13, and also on Wegner's formulation of the perturbative renormalization group [17-19]. 
Of course, all of this work rests on the early development of the renormalization group, 
especially on the work of Gell-Mann and Low [2] ; as well as on the ideas of Kadanoff that 
inspired the modern renormalization group [20]. 

The most promising area for application of the renormalization group in the study of the 
strong interaction is currently provided by the Euclidean lattice formulation of quantum 
chromo dynamics (QCD) [21], but here the nonperturbative renormalization group has 
found limited application and one is still forced to directly include short distance scales 
in large numerical calculations. The lattice itself introduces numerical complications that 
are not easily overcome, and alternative nonperturbative tools should be developed. At 
this point there is no serious challenger to lattice field theory, as all other nonperturbative 
algorithms rely on uncontrolled 'approximations.' 

The most significant barrier to the application of the renormalization group is algebraic 
complexity. This complexity stems partially from the general nature of the renormaliza- 
tion group as formulated by Wilson. The original renormalization group formalism devel- 
oped by Stueckelberg and Petermann [1], and Gell-Mann and Low [2], as well as direct 
derivatives such as that of Callan and Symanzik [22,23], are tailored to the problem of 
renormalizing canonical field theories, and take advantage of tools that have been devel- 
oped for Feynman perturbation theory [24] . These versions are of limited utility for some 
problems, particularly those that cannot be adequately solved with Feynman perturbation 
theory and those in which one needs to remove degrees of freedom with explicit cutoffs. In 
light-front field theory we encounter a problem that requires the more general renormal- 
ization group, because we need to use cutoffs that reduce the size of Fock space to attack 
nonperturbative problems; and because all cutoffs at our disposal violate symmetries such 
as Lorentz covariance and gauge invariance. Furthermore, we do not want to include a 
complicated vacuum in our state vectors, so we need the more general formulation of the 
renormalization group to allow interactions induced by the vacuum to directly enter our 
Hamiltonians. 

The modern renormalization group is a pragmatic approach to any problem that in- 
volves very many degrees of freedom that can be profitably divided according to a distance 
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or energy scale. Instead of trying to solve such problems by considering all scales at once, 
which usually fails even in perturbation theory, one breaks the problem into pieces, 'solving' 
each scale in sequence. There are many numerical advantages to this approach. 

As currently formulated the nonperturbative renormalization group [13] is reminiscent 
of the calculus as applied by Newton [25] . Most physicists would find it impossible to read 
the Principia, and few would recognize the calculus in the form that Newton and his peers 
employed it. Fortunately, most of us can avoid the Principia; but while there are many 
good introductions to the perturbative renormalization group [26-31], the nonperturbative 
renormalization group has not yet been developed to the point where simple introductions 
exist. This paper follows Wilson's remarkable review article on the renormalization group 
and its application to the Kondo problem [13]. I concentrate only on the development 
of a perturbative light-front renormalization group, and do not discuss possible nonper- 
turbative renormalization groups. One of the chief purposes of this paper, however, is to 
outline some of the problems one must face when developing a nonperturbative light-front 
renormalization group. 

Dirac formulated light-front field theory [32] during his unsuccessful search for a rea- 
sonable Hamiltonian formulation of relativistic field theory [33]. For the most complete 
set of references available on light-front physics, see Ref. 34. Unfortunately, Dirac did not 
follow through with his initial development of light-front field theory and it was largely 
ignored until Weinberg developed the closely related infinite momentum frame formalism 
[35]. The principal advantages of the light- front formalism are that boost invariance is 
kinematic, and that the bare vacuum mixes only with modes that have identically zero 
longitudinal momentum. The first advantage allows one to factor center-of-mass momenta 
from the equations of motion, which may be extremely important if it proves possible to 
formulate any relativistic problem so that it becomes a few body problem; a fact that 
can be appreciated after a study of the few-body problem in nonrelativistic quantum me- 
chanics. The second advantage is sometimes misrepresented to mean that the vacuum is 
trivial in light-front field theory. What is actually implied is that one can isolate all modes 
that mix with the trivial vacuum to form the physical vacuum, and then note that all 
of these have infinite energy in light-front field theory (ignoring a set of measure zero in 
theories with massless particles; i.e., states with identically zero transverse momentum). 
This observation may indicate that it is possible to replace the problem of building the 
physical vacuum with the problem of renormalizing the light-front Hamiltonian. Using a 
boost-invariant renormalization group one may be able to embed the vacuum problem into 
the larger problem of using the light-front renormalization group to remove high energy 
degrees of freedom. 

The light-front formalism does not automatically solve the physical problems that 
force one to actually build the vacuum in Euclidean field theory. It merely allows us to 
reformulate these problems in what will hopefully prove to be a more tractable form. If 
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a symmetry is broken by the vacuum, we are not allowed to assume that symmetry when 
using the renormalization group to construct the renormalized Hamiltonian. In other 
words, we are not allowed to use a broken symmetry to restrict the space of Hamiltonians, 
even when the symmetry is broken by the vacuum. For example, if we want to study 
(/> 4 theory in 1+1 dimensions beyond the critical coupling (i.e., in the symmetry broken 
phase) , we must allow the space of Hamiltonians to include symmetry breaking interactions 
(e.g., a </> 3 interaction). If such interactions are not allowed, states with imaginary mass 
will typically appear in the spectrum [36-38]. The renormalization group allows one to 
isolate relevant and marginal symmetry breaking interactions that might be tuned to 
reproduce the vacuum effects. A simple example involving spontaneous symmetry breaking 
has recently been provided by Wilson and the author [39]. 

While the infinite momentum frame formalism has been widely used, especially in the 
study of current algebra [40,41,34] and the parton model [42-44,34], little formal work has 
been completed on renormalization in light-front field theory, and almost all early work 
concentrates primarily on developing a map between light-front field theory and equal- 
time field theory in perturbation theory [45-52]. More recently a number of theorists 
have begun to study light-front perturbation theory directly [53-59], and especially the 
renormalization of light-front field theory after a Tamm-Dancoff [60-62] truncation is made 
[63-71] . The only formalism currently available for nonperturbative renormalization is the 
renormalization group, and work on developing the renormalization group for light-front 
field theory is in its infancy [73-75,59,39]. In my opinion, without a light-front version of the 
renormalization group, light-front field theory may be relegated to being a tool of last choice 
for doing perturbative calculations in 3+1 dimensions. Of course, in superrenormalizable 
theories in 1+1 dimensions light-front field theory has already proven to be very successful 
[76-82] ; and it can be argued that light-front field theory is a much more powerful tool for 
many nonperturbative calculations in 1+1 dimensions than equal-time field theory. 

In 3+1 dimensions we are in a situation where, because of serious renormalization 
difficulties, we neither know the correct light-front Hamiltonians to study, nor can we 
compute the physical ground states to which these unknown Hamiltonians lead. This is 
exactly the type of problem for which the modern renormalization group is suited. 

The primary purpose of this paper is to develop a perturbative light-front renormal- 
ization group as a tool for the study of light-front Hamiltonian field theory, with the hope 
that this work may aid the development of a practical nonperturbative light-front renor- 
malization group. The most interesting theory that one may be able to study with a 
perturbative light-front renormalization group is QCD; however, the algebraic complexity 
of QCD makes it a poor development ground. Following tradition, I use scalar field the- 
ory for all of my examples, as it is straightforward (i.e., difficult but not impossible) to 
generalize the formalism to other theories. 

In the remainder of this Introduction I outline the rest of the paper. In the process I use 
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both light-front and renormalization group jargon, often without offering any definition. I 
have tried to carefully define most of the renormalization group jargon in the text, so the 
reader who is unfamiliar with the modern renormalization group may want to read this 
Introduction again after reading the rest of the text. There are a number of articles that 
introduce or review most of the basic aspects of light-front field theory required in this 
work, and the reader unfamiliar with this formalism may want to consult one or more of 
these [46-51,83-86,65]. 

To implement a renormalization group calculation one must first delineate a space of 
Hamiltonians, and then define a renormalization group transformation [20] that maps this 
space into itself. The renormalization group transformation must be carefully formulated 
because it is central to the whole approach. Several renormalization group transformations 
are given for light-front field theory, including some boost-invariant transformations. These 
latter transformations allow one to impose the constraint of boost invariance directly on 
the space of Hamiltonians. The next step in a renormalization group study is to explore the 
topology of the Hamiltonian space, first searching for fixed points of the transformation 
(i.e., Hamiltonians that remain fixed under the action of the transformation) and then 
studying the trajectories of Hamiltonians near these fixed points. The Hamiltonian itself 
was originally formulated as a powerful tool for the study of physical trajectories found in 
nature. The renormalization group can be considered to be a generalization of this idea, 
in which the renormalization group transformation is used to derive physical Hamiltonians 
found in nature. As such, it is an alternative to the canonical procedure that starts with 
classical equations of motion. 

Almost all analytic work concentrates on Gaussian fixed points (i.e., Hamiltonians with 
no interactions) and near-Gaussian fixed points where perturbation theory can be used to 
approximate the renormalization group transformation itself. All of the examples in this 
paper are of this type, and the entire investigation is directed toward the development of 
transformations near Gaussian fixed points. Actually, there are fixed volumes rather than 
isolated fixed points, but I usually refer to fixed points rather than fixed volumes. Such 
examples may be of immediate relevance for QCD; and they illustrate the basic light-front 
renormalization group machinery. 

In Section II I provide a brief summary of the modern renormalization group formal- 
ism and the generalizations required when there are an infinite number of relevant and 
marginal operators. The generalizations are not easily appreciated until one has studied 
the entire paper and understood why they are required. This Section is a poor substitute 
for Wilson's review article [13], but I have attempted to introduce the most important 
concepts required by a perturbative renormalization group. I have also attempted to 
make the differences between a perturbative renormalization group and a nonperturbative 
renormalization group clear, focusing on the former. The modern renormalization group 
formalism may be unfamiliar to many students of relativistic field theory, being employed 
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primarily in the study of critical phenomena. With a few notable exceptions [29,30], most 
field theory textbooks deal exclusively with the original renormalization group formalism 
and its modern descendant, the Callan-Symanzik formalism [22,23]. Perhaps the most 
important difference between these two types of renormalization group is that the original 
renormalization group does not actually remove any degrees of freedom, being concerned 
primarily with the problem of divergences in perturbation theory and techniques for al- 
lowing all cutoffs to be removed. Here one typically uses either Pauli-Villars regularization 
[87], which actually increases the number of degrees of freedom, or dimensional regular- 
ization [88,89], which retains all degrees of freedom while analytically continuing in the 
number of dimensions. Neither of these regulators is well-suited to many nonperturbative 
calculations. 

In order to reduce the number of degrees of freedom, one must introduce a real cutoff, 
such as a momentum cutoff or a lattice cutoff. The cutoff is an artifice, so results should 
not depend on its particular value and it should be possible to change the cutoff without 
changing any physical result; i.e., without changing the matrix elements of observables 
between cutoff physical states. (Later I often refer to such matrix elements as observables, 
since the only operator that is discussed in any detail is the Hamiltonian.) This can be 
accomplished by making the observables, and in particular the Hamiltonian, depend on 
the cutoff in precisely the manner required to yield cutoff independent matrix elements. 
The modern renormalization group is designed to achieve this goal. 

In Section III I illustrate what is meant by a space of Hamiltonians, and provide 
several light-front renormalization group transformations. I begin by discussing transfor- 
mations that resemble those developed by Wilson for Euclidean field theory; however, these 
transformations later lead to pathologies because they remove states of lower energy than 
some that are retained. Almost all light-front transformations suffer from these patholo- 
gies in perturbation theory, and in the end I am forced to consider renormalization groups 
with some unusual properties to obtain a transformation that may be approximated by 
discarding at least some irrelevant operators in perturbation theory. These transforma- 
tions employ invariant-mass cutoffs, so I refer to them as invariant-mass transformations. 
The restriction to transformations that may be approximated perturbatively is an extreme 
limitation, and it is not even clear that the invariant-mass transformations can be ap- 
proximated in perturbation theory with controlled errors. We will find that couplings 
depend on longitudinal momentum fractions when one uses a boost-invariant cutoff, and 
that corrections to the Hamiltonian diverge logarithmically for states containing particles 
with arbitrarily small longitudinal momentum fractions, and for interactions involving ar- 
bitrarily small longitudinal momentum transfer. Thus, ultimately it is not clear that the 
invariant-mass transformations are 'better' than other transformations one may use; and 
one should certainly consider other transformations when developing a nonperturbative 
renormalization group. This article does not attack these nonperturbative problems, even 



6 



though they may be of more interest than the perturbative results obtained here; how- 
ever, I try to clarify some of the nonperturbative renormalization problems light-front field 
theorists should be attacking. 

An essential restriction in all Euclidean renormalization group calculations is that long 
range interactions are excluded. As Wilson notes [13], this is one of the most tenuous 
assumptions of the renormalization group approach. This locality assumption must be 
altered in light-front field theory where inverse powers of longitudinal derivatives are re- 
quired already at the Gaussian fixed points of interest; i.e., in free field theory. Allowing 
inverse powers of longitudinal derivatives may seem to be a rather minor conceptual mod- 
ification of the Euclidean version of the renormalization group; but the recognition that 
there are separate renormalization group transformations that act on the longitudinal and 
transverse directions may be profound. This generalizes the fact that there are separate 
power counting analyses for longitudinal and transverse dimensions. The analysis of rel- 
evant, marginal and irrelevant longitudinal operators indicates that it is inverse powers 
of longitudinal derivatives that arise as irrelevant operators when one of the light-front 
renormalization group transformations is applied to a Hamiltonian. Having allowed in- 
verse longitudinal derivatives, one naturally considers the possibility that inverse powers 
of transverse derivatives occur, a possibility that is especially intriguing for the study of 
QCD; however, it is difficult to introduce such operators in a controlled manner, and I 
avoid their introduction in this article. 

I assume that transverse interactions are local, or at least short range, and refer to this 
assumption as transverse locality. This restriction is not merely a technical convenience, 
because inverse transverse derivatives typically lead to an infinite number of relevant op- 
erators near critical Gaussian fixed points, including products of arbitrarily large numbers 
of field operators. Near fixed points that contain interactions, such operators may ap- 
pear without causing trouble, but I discuss only Gaussian fixed points. Once the space of 
Hamiltonians is restricted, the renormalization group transformation may produce Hamil- 
tonians that lie outside the space. We are only interested in trajectories of Hamiltonians 
generated by repeated application of the transformation that remain inside the restricted 
space. Strict transverse locality is violated by the step function cutoffs I employ; however, 
these violations appear to be controllable. Moreover, inverse transverse derivatives are 
generated by the transformations; but they are accompanied by cutoffs and I show that 
the resultant distributions do not produce long range transverse interactions, and that they 
do not introduce relevant operators. In all cases that I have discovered where a product of 
inverse transverse derivatives and cutoff functions arise, the product can be shown to be an 
irrelevant operator with respect to transverse scaling; although the operator may contain 
delta functions or derivatives of delta functions of longitudinal momentum fractions. 

Each of the light-front renormalization group transformations I consider consists of 
two steps. In the first step one alters a cutoff {e.g., lowers a cutoff on transverse momenta) 
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to remove degrees of freedom, and computes a renormalized Hamiltonian that produces 
the same eigenvalues and suitably ortho normalized, projected eigenstates, in the remain- 
ing subspace. All operators that correspond to observables are renormalized, not just the 
Hamiltonian, but in this paper I focus only on the Hamiltonian. Of particular interest for 
future work is the renormalization of other Poincare generators and various current oper- 
ators. In the second step of each transformation the variables (e.g., transverse momenta) 
are rescaled to their original range, and the fields and Hamiltonian are rescaled. The most 
difficult part of this procedure is the construction of an effective Hamiltonian, which is 
analogous to the development of a block spin Hamiltonian in simple spin systems. In Sec- 
tion IV I discuss two related procedures for accomplishing this task; the first pioneered 
by Bloch and Horowitz [90], and the second by Bloch [91,92]. The second method was 
employed by Wilson in his first serious numerical renormalization group study [5]. 

In Section V I turn to the study of fixed points in perturbation theory and linearized 
behavior of the transformations near these fixed points [17-19], developing simple examples 
that hopefully clarify the procedure. The study of linearized behavior near critical Gaussian 
fixed points in light-front field theory is dimensional analysis, as usual. The fact that 
longitudinal boost invariance corresponds to an exact scale invariance [32,44,93] leads to 
the conclusion that all physical Hamiltonians are fixed points with respect to a longitudinal 
light-front renormalization group transformation. This should be true order-by-order in 
any perturbative expansion of the Hamiltonian in powers of a parameter upon which it 
depends analytically, and it should be an extremely powerful tool for the analysis of fixed 
points. However, the pathologies mentioned above make it difficult to find applications. 
Transformations that scale only the transverse momenta produce an infinite number of 
relevant and marginal operators, in addition to the familiar infinity of irrelevant operators. 
This happens because entire functions of longitudinal momentum fractions may appear in 
any given operator without affecting the linear analysis that determines this classification. 
Transformations that scale only longitudinal momenta also lead to an infinite number of 
relevant and marginal operators because entire functions of transverse momenta are allowed 
to appear. The appearance of entire functions drastically complicates the renormalization 
group analysis, but it may also eventually lead to tremendous power in the application of 
the renormalization group if one can learn how to accurately approximate these functions. 

In Section VI I study second-order perturbations about the critical Gaussian fixed 
point. I concentrate on Hamiltonians near the canonical massless (p 4 Hamiltonian, and 
first show that several candidate transformations lead to divergences at second-order. I 
turn to a boost-invariant transformation and concentrate primarily on a few marginal and 
relevant operators. I show that it is possible to find a closed set of two marginal operators 
and one relevant operator in a second-order analysis, despite the possible appearance of an 
infinite number of such operators. This simple analysis illustrates many of the features of 
a full transformation. I then study Lorentz covariance and cluster decomposition, showing 
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that both are violated by the invariant-mass cutoff in second-order perturbation theory. 
These properties must be restored by counterterms, and the formalism should be able to 
produce these counterterms without referring to covariant results. 

Since there are an infinite number of relevant and marginal operators in the light- 
front renormalization group, and a simple perturbative analysis indicates that Lorentz 
covariance and cluster decomposition cannot be restored without actually employing an 
infinite number of such operators, one must worry that the light-front renormalization 
group analysis requires one to adjust an infinite number of independent variables. Wilson 
and I have recently proposed conditions under which a finite number of variables are 
actually independent [39] . These conditions turn out to be a generalization of the coupling 
reduction conditions first developed by Oehme, Sibold, and Zimmerman [94-100]. 

Simply stated, we have proposed that one should seek solutions to the renormalization 
group equations in which a finite number of specified variables are allowed to be inde- 
pendent functions of the cutoff. While there are an infinite number of relevant, marginal, 
and irrelevant variables, all but a finite number of variables depend on the cutoff only 
through their dependence on these independent variables. This condition could merely 
be a re-parameterization of the cutoff dependence, but we also add the constraint that 
all dependent variables must vanish when the independent variables are zero. This last 
constraint is motivated by considering the dependent variables to be counterterms. This 
means that as one changes the cutoffs, all couplings (including masses) evolve coherently; 
whereas the general solution to the renormalization group equations might allow much 
more complicated behavior, at least for the relevant and marginal variables. It is for this 
reason that we call our conditions coupling constant coherence conditions, or more briefly, 
coupling coherence. 

In Section VI I provide part of the demonstration that coupling coherence fixes the 
strengths of all operators to second order in the canonical coupling, 0(g 2 ), and that the 
resultant strengths are precisely the values required to restore Lorentz covariance and 
cluster decomposition in second-order perturbation theory. I show that coupling coher- 
ence uniquely fixes the relevant mass operator, and the dispersion relation associated with 
the bare mass term is not that of a free massive particle. I determine the complete set of 
irrelevant four-boson couplings, and show that a third-order (i.e., two loop) calculation is 
required to fix the marginal four-boson couplings. I do not compute all couplings to 0(g 2 ), 
but it is straightforward to complete the analysis for the terms I do not evaluate. I am pri- 
marily interested in the marginal four-boson couplings, because these indicate how one can 
expect couplings to run with longitudinal momenta in a light-front renormalization group 
analysis. In the scalar theory, couplings decrease in strength as the longitudinal momen- 
tum fraction carried by the bosons decreases, and as the longitudinal momentum fraction 
transferred through the vertex decreases. In QCD one expects the opposite behavior. 

There are two essential expansions that are made in a perturbative renormalization 
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group analysis. The first is the perturbative expansion of the transformation itself, which 
may converge sufficiently near the Gaussian fixed point. The second is the expansion of the 
transformed Hamiltonian in terms of relevant, marginal, and irrelevant operators. Each 
expansion must normally be truncated at some finite order, and one should try to show that 
each truncation leads to controllable errors. Much of the discussion of the errors introduced 
by truncating the perturbative expansion is delayed to Section VII, although some of the 
most important sources of errors are already evident in a second-order analysis and are 
discussed in Section VI. Many of the errors that arise when one truncates the expansion of 
the Hamiltonian in terms of relevant, marginal and irrelevant operators by discarding all 
or most of the irrelevant operators can be studied using the second-order approximation 
of the transformation. Dropping some irrelevant operators is essential to the program, 
because the transformations become too complicated if one must follow the evolution of 
too many operators. 

Most transformations one might construct lead to uncontrollable errors even in a 
second-order analysis for a simple reason. When a transformation is applied once, irrele- 
vant operators are produced; and when the transformation is applied again these irrelevant 
operators produce divergences. The pathology is easily analyzed. Simple transformations 
attempt to remove degrees of freedom with much lower free energy than some of the degrees 
of freedom retained. As one step in making the expansion in terms of relevant, marginal 
and irrelevant operators, one typically expands every energy denominator encountered in 
the perturbative expansion of the transformation. This expansion is in powers of one of 
the energies of a state that is retained after the transformation, and if this state has a 
higher energy than one of the states removed, the expansion of the energy denominator 
fails to converge. Simply stated, l/(Ei n — Ei oop ) cannot be expanded in powers of E{ n when 
Ei n is larger than Ei oop . I believe that there are few solutions to this problem. One can 
abandon the perturbative expansion of the transformation, one can abandon the renor- 
malization group approach or drastically modify it, or one can design the transformation 
so that Ei n is always less than Ei oop . In this paper, I choose the final option. The simplest 
boost-invariant transformations suffer from this same problem, and only the invariant-mass 
transformation escapes. 

In Section VII I discuss third- and higher-order corrections to an invariant-mass trans- 
formation near the critical Gaussian fixed point. I first show that the third-order analysis 
introduces new marginal operators that contain logarithms of longitudinal momentum 
fractions. I also show that when one insists that these new marginal operators depend 
on the cutoff only through their dependence on the original marginal </> 4 interaction, their 
strength is precisely that required to restore Lorentz covariance and cluster decomposition 
to the boson-boson scattering amplitude computed in second-order perturbation theory. I 
then turn to a discussion of the errors introduced by various approximations, showing that 
these errors may be quite large. 
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Wilson has provided a thorough discussion of perturbative renormalization group equa- 
tions and the errors that result from various approximations, so I focus on new approxi- 
mations that must be made in a light-front renormalization group analysis. In a typical 
Euclidean renormalization group there are a finite number of relevant and marginal op- 
erators and one can accurately tune their boundary values at the lowest cutoff. In the 
light-front renormalization group functions of longitudinal momenta appear in each rele- 
vant and marginal operator, and these functions must be approximated. This means that 
errors are made in the relevant and marginal operators themselves, and one must determine 
whether these errors can be controlled. 

These new approximations are extremely problematic, because one expects that the 
strength of a relevant operator will grow exponentially as the cutoff is lowered; which means 
that any error made in the approximation of a relevant operator will grow exponentially. 
Of course, one partial solution to this problem is to approximate the relevant operators at 
the lowest cutoff, and solve the renormalization group equations for the relevant operators 
in the direction of increasing cutoff, as is usually done [13]; because the strengths of 
relevant operators decrease exponentially when the cutoff is increased. Marginal operators 
present a more difficult problem, as usual; because errors tend to grow linearly regardless 
of the direction in which one solves the equations for marginal operators. Moreover, we 
will see that operators develop logarithmic singularities when one uses the invariant-mass 
transformation. These problems require careful study, and this paper barely initiates such 
a study. No convincing solution to these problems is proposed in this paper, but I have tried 
to present the problems in a clear manner; because any attempt to use Hamiltonian light- 
front field theory to perform nonperturbative calculations must address such problems. 
I offer some speculation on how one might approach this task, and argue that the same 
problems that complicate the perturbative analysis may actually lead to a simplification 
of the nonperturbative analysis. 

The renormalization group analysis is drastically simplified when a Tamm-Dancoff 
truncation [60,61] is made and the light-front Tamm-Dancoff (LFTD) approximation [63- 
71] is used, so I occasionally mention important aspects of a renormalization group analysis 
of LFTD; however, I do not consider such examples in this paper. The Tamm-Dancoff 
truncation can be simply included in the initial definition of the space of Hamiltonians, 
after which the analysis proceeds exactly as when no truncation is made. The truncation 
preserves boost invariance, and thereby allows the boost-invariant renormalization group 
transformations to be illustrated. More importantly, it drastically simplifies the operators 
that are included in the space of Hamiltonians, even allowing one to solve some simple ex- 
amples analytically. Sector-dependent renormalization, in which parameters appearing in 
the Hamiltonian and other observables are allowed to depend on the Fock space sector(s) 
within or between which they act, arises in a natural manner; and the light-front renormal- 
ization group drastically improves the discussion of renormalization in LFTD. However, if 
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one wants to use LFTD to study physical theories, it is probably necessary to remove some 
of the restrictions on the space of allowed Hamiltonians that I assume, in particular re- 
strictions associated with transverse locality. I do not yet know how one can introduce the 
required nonlocalities and still control the number of operators required; however, simple 
perturbation theory arguments show that nonlocal transverse operators inevitably arise if 
one derives a LFTD Hamiltonian by eliminating states with extra particles. These issues 
are not discussed in this paper. 

In the conclusion I discuss the possible relevance of this work for the study of Hamil- 
tonian light- front QCD [101-104,83,85,86], indicating some of the difficult problems that I 
carefully avoid with my simple examples in this paper. 

II. The Renormalization Group 

In classical mechanics the state of a system is completely specified by a fixed number 
of coordinates and momenta. The objective of classical mechanics is to compute the state 
as a function of time, given initial conditions. The state is not regarded as fundamental; 
rather a Hamiltonian that governs the time evolution of the state is regarded as funda- 
mental. In nonrelativistic quantum mechanics, one must generalize the definition of the 
state, so that it is specified by a ket in a state space, and one must drastically alter the 
theory of measurement; but it remains possible to specify a Hamiltonian that governs the 
time evolution of the state. In both cases the time evolution of the state is a trajectory in 
a Hilbert space, and the trajectory is determined by a Hamiltonian that must be discov- 
ered by fitting data. In principle, one would like to further generalize this procedure for 
relativistic field theory; however, any straightforward generalization that maintains local- 
ity leads to divergences that produce mathematical ambiguities. To make mathematically 
meaningful statements we must introduce an ad hoc regulator, to which I refer as a cutoff; 
so that physical results can be derived as limits of sequences of finite quantities. The 
renormalization group provides methods for constructing such limits that are much more 
powerful than standard perturbation theory. 

Even if divergences did not signal the need for a cutoff in field theory, we would be 
forced to introduce a cutoff in some form to obtain finite dimensional approximations 
for the state vector and Hamiltonian. Fock space is an infinite dimensional sum of cross 
products of infinite dimensional Hilbert spaces, and this is not a convenient starting point 
for most numerical work. If all interactions are weak and of nearly constant strength 
over the entire range of scales that affect an observable, we can use standard perturbation 
theory to compute the observable; however, if either of these conditions is not met, we 
cannot directly compute observables with realistic Hamiltonians. This problem is easily 
appreciated by considering a simple spin system in which 1024 spins are each allowed 
to take two values. The Hamiltonian for this system is a 2 1024 x 2 1024 matrix, and this 
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matrix cannot generally be diagonalized directly. This matrix is infinitely smaller than 
the matrices we must consider when solving a relativistic field theory. If the interactions 
remain weak over all scales of interest, but change in strength significantly, we can use the 
perturbative renormalization group. If the interactions become strong over a large number 
of scales of interest, a nonperturbative renormalization group must be developed. A final 
possibility is that the interactions are weak over almost all scales, becoming strong only 
over a few scales of interest. In this case, the perturbative renormalization group can be 
used to eliminate the perturbative scales; after which one can use some other method to 
solve the remaining Hamiltonian. 

The introduction of an ad hoc cutoff in field theory complicates the basic algorithm 
for computing the time evolution of a state, because one must somehow remove any de- 
pendence on the cutoff from physical matrix elements. This complication is so severe that 
it has caused field theorists to essentially abandon many of the most powerful tools em- 
ployed in nonrelativistic quantum mechanics (e.g., the Schrodinger picture). How can we 
make reasonable estimates in relativistic quantum mechanics? How can we guarantee that 
results are independent of the cutoff? How can we find a sequence of Hamiltonians that 
depends on the cutoff in a manner that leads to correct results as the cutoff approaches 
its limit? These are the type of questions that led Wilson to completely reformulate the 
original Gell-Mann-Low renormalization group formalism. 

Wilson adopted the same general strategy familiar from the study of the time evolution 
of states, adding a layer of abstraction to the original classical mechanics problem to 
compute 'Hamiltonian trajectories'. The strategy is universal in physics, but the layer of 
abstraction leads to a great deal of confusion. In analogy to a formalism that yields the 
evolution of a state as time changes, he developed a formalism that yields the evolution 
of a Hamiltonian as the cutoff changes. In quantum mechanics a state is represented 
by an infinite number of coordinates in a Hilbert space, and the Hamiltonian is a linear 
operator that generates the time evolution of these coordinates. In the renormalization 
group formalism, the existence of a space in which the Hamiltonian can be represented by 
an infinite number of coordinates is assumed, and the cutoff evolution of these coordinates 
is given by the renormalization group transformation. The Hamiltonian is less fundamental 
than the renormalization group transformation, which can be used to construct trajectories 
of Hamiltonians. Typically there are restrictions placed on the Hamiltonians (e.g., no 
long-range interactions) that make it possible for these trajectories to leave the space. 
Trajectories of renormalized Hamiltonians remain in the space of Hamiltonians, and are 
roughly analogous to physical trajectories in classical mechanics that meet some additional 
requirement (e.g., do not leave a specified finite volume). Transformations that change the 
cutoff by different amounts are members of a renormalization 'group', which is actually a 
semi-group since the transformations cannot be inverted. The fact that inverses do not 
exist is obvious because the transformations reduce the number of degrees of freedom. 
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A more thorough, although by no means complete, discussion of the space of Hamilto- 
nians is given in Section III. Here I will simply state that each term in a Hamiltonian can 
typically be written as a spatial integral whose integrand is a product of derivatives and 
field operators. The definition of the Hamiltonian space might be a set of rules that show 
how to construct all allowed operators. These operators should be thought of as unit vec- 
tors, and the coefficients in front of these operators as coordinates. This type of operator 
is not usually bounded, and this is a source of divergences in field theory. To regulate these 
divergences the cutoff is included directly in the definition of the space of Hamiltonians. 
The cutoff one chooses has drastic effects on the renormalization group. While we will 
see several examples in Section III, one familiar example of a cutoff is the lattice, which 
replaces spatial integrals by sums over discrete points. The facts that the Hamiltonian can 
be represented by coordinates that correspond to the strengths of specific operators, and 
that these operators are all regulated by a cutoff that is part of the definition of the space, 
is all that one needs to appreciate at this point. 

Given a space of cutoff Hamiltonians, the next step is to construct a suitable trans- 
formation. This is slightly subtle and is usually the most difficult conceptual step in a 
renormalization group analysis, because one must find a transformation that manages to 
alter the cutoff without changing the space in which the Hamiltonian lies. These two 
requirements seem mutually contradictory at first. An additional problem for relativistic 
field theory is that all transformations one can construct change the cutoff in the wrong 
direction. 

To see how these difficulties are averted, let me again use the lattice as an example. 
A typical lattice transformation consists of two steps. In the first step one reduces the 
number of lattice points, typically by grouping them into blocks [20] and thereby doubling 
the lattice spacing; and one computes a new effective Hamiltonian on the new lattice. I 
do not discuss how this effective Hamiltonian is constructed for a lattice, but this issue is 
carefully discussed for light-front Hamiltonians in later Sections. At this point the lattice 
has changed, so the space in which the Hamiltonian lies has changed. The second step in 
the transformation is to rescale distances using a change of variables, so that the lattice 
spacing is returned to its initial value, while one or more distance units are changed. After 
both steps are completed the lattice itself remains unchanged, if it has an infinite volume, 
but the Hamiltonian changes. This shows how one can alter the cutoff without leaving 
the initial space of Hamiltonians. Numerically the cutoff does not change, but the units 
in which the cutoff is measured change. 

We want to study the limit in which the lattice spacing is taken to zero, but the 
transformation increases the lattice spacing as measured in the original distance units. 
While no inverse transformation that directly decreases the lattice spacing exists, we can 
obtain the limit in which the lattice spacing goes to zero by considering a sequence of 
increasingly long trajectories. Instead of fixing the initial lattice spacing, we fix the lattice 
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spacing at the end of the trajectory and we construct a sequence of trajectories in which 
the initial lattice spacing is steadily decreased. We then directly seek a limit for the last 
part of the trajectory as it becomes infinitely long, by studying the sequence of increasingly 
long trajectories. 

This procedure is illustrated by Wilson's triangle of renormalization, which is briefly 
discussed below. One must employ Wilson's algorithm to perform a nonperturbative renor- 
malization group analysis; however, it is possible to study the cutoff limit more directly 
when a reasonable perturbative approximation exists. In this case, the renormalization 
group transformation can be approximated by an infinite number of coupled equations 
for the evolution of a subset of coordinates that are asymptotically complete, and these 
equations can be inverted to allow direct study of the Hamiltonian trajectory as the cutoff 
increases or decreases [17-19]. If it can be shown that all but a finite number of coor- 
dinates remain smaller than some chosen magnitude, it may be possible to approximate 
the trajectory by simply ignoring the small coordinates, retaining an increasing number 
of coordinates only as one increases the accuracy of the approximation. In this case the 
task of approximating a trajectory of renormalized Hamiltonians is reduced to the task of 
solving a finite number of coupled nonlinear difference equations. The primary goal of this 
paper is the development of a perturbative light-front renormalization group. 

Given a transformation T that maps a subspace of Hamiltonians into the space of 
Hamiltonians, with the possibility that some Hamiltonians are mapped to Hamiltonians 
outside the original space, we study T[H]. We can apply the transformation repeatedly, 
and construct a trajectory of Hamiltonians, with the Z-th point on the trajectory being 



Any infinitely long trajectory that remains inside the space is called a trajectory of renor- 
malized Hamiltonians. The motivation for this definition of renormalization is clarified 
further below. It is assumed that the trajectory is completely determined by the initial 
Hamiltonian, H$, and T; however, the dependence on Hq is usually not explicitly indi- 
cated. Moreover, we will see later that boundary conditions may be specified in a much 
more general fashion. 

Any renormalization group analysis begins with the identification of at least one fixed 
point, H*. A fixed point is defined to be any Hamiltonian that satisfies the condition 



For perturbative renormalization groups the search for such fixed points is relatively easy, 
as we will see in Section V; however, in nonperturbative studies such a search typically 



Hi = T l [H ]. 



(2.1) 



H* = T[H*] . 



(2.2) 
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involves a difficult numerical trial and error calculation [5,7,13]. If H* contains no interac- 
tions (i.e., no terms with a product of more than two field operators), it is called Gaussian. 
If H* has a massless eigenstate, it is called critical. If a Gaussian fixed point has no mass 
term, it is a critical Gaussian fixed point. If it has a mass term, this mass must typically 
be infinite, in which case it is a trivial Gaussian fixed point. In lattice QCD the trajectory 
of renormalized Hamiltonians stays near a critical Gaussian fixed point until the lattice 
spacing becomes sufficiently large that a transition to strong-coupling behavior occurs. If 
H* contains only weak interactions, it is called near- Gaussian, and one may be able to 
use perturbation theory both to identify H* and to accurately approximate trajectories 
of Hamiltonians near H* [12]. Of course, once the trajectory leaves the region of H* it 
is generally necessary to switch to a nonperturbative calculation of subsequent evolution. 
If H* contains a strong interaction, one must use nonperturbative techniques to find H*, 
but it may still be possible to produce trajectories near the fixed point using perturbation 
theory. The perturbative analysis in this case includes the interactions in H* to all orders, 
treating only the deviations from these interactions in perturbation theory. 

Consider the immediate neighborhood of the fixed point, and assume that the trajec- 
tory remains in this neighborhood. This assumption must be justified a posteriori, but if 
it is true we should write 



H^H' + SHt, (2.3) 

and consider the trajectory of small deviations 5 Hi. 

As long as 5 Hi is 'sufficiently small', we can use a perturbative expansion in powers of 
6 Hi, which leads us to consider 

5Hi +1 = L-5Hi+N[5Hi}. (2.4) 

Here L is the linear approximation of the full transformation in the neighborhood of the 
fixed point, and N{5H[\ contains all contributions to 5Hi + i of 0(5Hf) and higher. 

The object of the renormalization group calculation is to compute trajectories and this 
requires a representation for 5Hi. The problem of computing trajectories is one of the most 
common in physics, and a convenient basis for the representation of 5Hi is provided by 
the eigenoperators of L, since L dominates the transformation near the fixed point. These 
eigenoperators and their eigenvalues are found by solving 



L ■ O m — X m O m . (2-5) 

If H* is Gaussian or near-Gaussian it is usually straightforward to find L, and its eigenop- 
erators and eigenvalues. This is not typically true if H* contains strong interactions, and 
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in much of the remaining discussion I focus on formalism that is primarily useful for the 
study of trajectories in the neighborhood of Gaussian and near-Gaussian fixed points. In 
Section V the linear approximations of several simple light-front renormalization group 
transformations about a critical Gaussian fixed point are derived, and their eigenoperators 
and eigenvalues are computed. 

Using the eigenoperators of L as a basis we can represent 5H[, 

SHi = ^ fJmiOm + ^2 9™i° m + W ™i°™ ■ (2.6) 
meR rneM m€l 

Here the operators O m with m G R are relevant {i.e., X m > 1), the operators O m with 
m G M are marginal (i.e., A m = 1), and the operators with m & I are either irrelevant 
(i.e., X m < 1) or become irrelevant after many applications of the transformation. The 
motivation behind this nomenclature is made clear by considering repeated application 
of L, which causes the relevant operators to grow exponentially, the marginal operators 
to remain unchanged in strength, and the irrelevant operators to decrease in magnitude 
exponentially. There are technical difficulties associated with the symmetry of L and 
the completeness of the eigenoperators that I ignore [17-19]. I occasionally refer to the 
relevant variables as masses, and the marginal and irrelevant variables as couplings; but I 
also occasionally refer to all variables, including relevant variables, as couplings. What is 
meant should be clear from context. 

L depends both on the transformation and the fixed point, but there are always an 
infinite number of irrelevant operators. On the other hand, transformations of interest for 
Euclidean lattice field theory typically lead to a finite number of relevant and marginal 
operators. One of the most serious problems for a perturbative light-front renormalization 
group is that an infinite number of relevant and marginal operators are required. This 
result is derived in Section V, and I discuss some of the consequences below. In the case 
of scalar field theory, an infinite number of relevant and marginal operators arise because 
the light-front cutoffs violate Lorentz covariance and cluster decomposition. These are 
continuous symmetries, and their violation leads to an infinite number of constraints on 
the Hamiltonian. The key to showing that the light-front renormalization group may not be 
rendered useless by an infinite number of relevant and marginal operators is the observation 
that both the strength and the evolution of all but a finite number of relevant and marginal 
operators are fixed by Lorentz covariance and cluster decomposition. However, one does 
not want to employ either of these properties directly in the construction of Hamiltonians, 
because they are never explicit in the renormalization group calculation of a Hamiltonian 
trajectory. Lorentz covariance and cluster decomposition are properties of observations 
that are obtained using the full Hamiltonian, and one does not want to solve problems 
that require the entire Hamiltonian to compute the Hamiltonian itself. The alternative 
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that Wilson and I have proposed [39,94-100] is to insist that the new relevant and marginal 
variables are not independent functions of the cutoff, but depend only on the cutoff through 
their dependence on canonical variables. While this requirement obviously fixes the manner 
in which the new variables evolve with the cutoff, it also fixes their value at all cutoffs once 
the values of the canonical variables are chosen. The remarkable feature of this procedure 
is that the value it gives to the new variables is precisely the value required to restore 
Lorentz covariance and cluster decomposition. This conclusion is not proven to all orders 
in perturbation theory, but it is illustrated by a nontrivial second order example in this 
paper. 

To simplify subsequent discussion, the statement that 5Hi is small is assumed to mean 
that all masses and couplings in the expansion of 5Hi are small. The analysis itself should 
signal when this assumption is naive. A rigorous discussion would require consideration of 
the spectra of the eigenoperators. In Section VI I show that several candidate light-front 
renormalization group transformations lead to unbounded operators, including unbounded 
irrelevant operators, even though there are cutoffs. In some of these cases, corrections that 
are 0(5H 2 ) are shown to be infinite; i.e., not small. If the coefficient of a single operator 
(e.g., a single mass) becomes large, it may be straightforward to alter the analysis so that 
this coefficient is included to all orders in an approximation of the transformation, so that 
one perturbs only in the small coefficients; but this possibility is not pursued in this paper. 

For the purpose of illustration, let me assume that A m = 4 for all relevant operators, 
and A m = 1/4 for all irrelevant operators. The transformation can be represented by an 
infinite number of coupled, nonlinear difference equations: 

mi ) 9mi 5 w mi ] , (2-7) 



9m l+1 =g mi +Ng m [n 

mi j 9mi , w mi ] , (2-8) 



W mi+i = 1 W mi +N Wm [li 

mi i 9 mi i w mi ] • (2-9) 

Sufficiently near a critical Gaussian fixed point, the functions N^ m , N 9m , and N Wm should be 
adequately approximated by an expansion in powers of fx mi , g m[ , and w mi . The assumption 
that the Hamiltonian remains in the neighborhood of the fixed point, so that all \i mi , 
g mi , and w mi remain small must be justified a posteriori. Any precise definition of the 
neighborhood of the fixed point within which all approximations are valid must also be 
provided a posteriori. 
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Wilson has given a general discussion of how these equations are solved [13], and 
I repeat only the most important points. In perturbation theory these equations are 
equivalent to an infinite number of coupled, first-order, nonlinear differential equations. 
To solve them we must specify 'boundary' values for every variable, possibly at different /, 
and then employ a stable numerical algorithm to find the variables at all other values of I 
for which the trajectory remains near the fixed point. We want to apply the transformation 
M times, letting M — > oo, and adjusting the initial Hamiltonian so that this limit exists. 
Eq. (2.7) must be solved by 'integrating' in the exponentially stable direction of decreasing 
/ (i.e., typically toward larger cutoffs), while Eq. (2.9) must be solved in the direction of 
increasing /. Eq. (2.8) is linearly unstable in either direction. The coupled equations 
must be solved using an iterative algorithm. Such systems of coupled difference equations 
and the algorithms required for their solution are familiar in numerical analysis. In this 
context the need for renormalization can be understood by considering the fact that the 
renormalization group difference equations need to be solved over an infinite number of 
scales in principle. 

The final output of the renormalization group analysis is the cutoff Hamiltonian Hj^. 
If this Hamiltonian is the final point in an infinitely long trajectory of Hamiltonians, it will 
yield the same observables below the final cutoff as Ho; but for an infinitely long trajectory 
Hq contains no cutoff, so Hj^f will yield results that do not depend on the cutoff. It is for this 
reason that Hj^ and all other Hamiltonians on any infinitely long trajectory are referred to 
as renormalized Hamiltonians. How one solves the final cutoff Hamiltonian problem using 
Hj\f depends on the theory. For the scalar theories used as examples in this paper I assume 
that perturbation theory can be used to predict observables. For QCD, even if Hj^f can 
be derived by purely perturbative techniques, it will have to be solved nonperturbatively 
because of confinement. In either case, we must have an accurate approximation for Hfj; 
however, we do not necessarily need to explicitly construct accurate approximations for all 
Hi. 

The boundary values for the irrelevant variables should be set at / = 0, because we 
need to solve Eq. (2.9) in the direction of increasing /. At large / all variables are ex- 
ponentially insensitive to the irrelevant boundary values. Therefore, they can be chosen 
arbitrarily (universality); and the values of the irrelevant variables at / = M are output 
by the modern renormalization group. This is one of the crucial differences between the 
modern renormalization group and the Gell- Mann-Low renormalization group in which 
irrelevant variables are not treated. Irrelevant operators are important in Hj^ unless the 
final cutoff is much larger than the scale of physical interest. The fact that they are ir- 
relevant implies that their final values are exponentially insensitive to their initial values; 
and it implies that they are driven at an exponential rate toward a function of the relevant 
and marginal variables, as discussed below. The fact that they are irrelevant does not 
necessarily imply that they are unimportant. This depends on how sensitive the physical 
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observables in which one is interested are to physics near the scale of the cutoff. 

The boundary values required by Eqs. (2.7) and (2.8) can be given at / = Af. Suffi- 
ciently far from / = the irrelevant variables are exponentially driven to maintain poly- 
nomial dependence on relevant and marginal variables, and sufficiently far from / = M 
the relevant variables are exponentially driven toward similar polynomial dependence on 
the marginal variables. While the calculation of transient behavior near / = and / = M 
usually requires a numerical computation, the relevant and irrelevant variables are read- 
ily approximated by polynomials that involve only marginal variables in the intermediate 
region. These polynomials are determined by the expansions of A^ m , N gm , and N Wm ; and 
they can be fed back into Eq. (2.8) to find an approximate equation for the marginal 
variables that requires direct knowledge only of the marginal variables themselves. 

These points may be confusing, so let me consider a simple example. Consider coupled 
differential equations for a relevant variable m, a marginal variable g, and an irrelevant 
variable w, 



— = -2m + cis 2 + c 2 gw , (2.10) 
at 



02 = -C30 3 + C4g 2 m + c 5 w 2 , (2.11) 



^ = 2w + c 6 g 2 . (2.12) 

Here the cutoff increases as t increases. We want to fix the boundary condition for w at 
t — > oo, and the boundary condition for m at t = 0. We are only interested in the solution 
near t = 0, and this should not depend on the boundary condition for w; as long as w(t) 
remains finite as t — > oo. Let m(0) = mo, and g(0) = go. 

To satisfy the boundary condition for w we must have 

w (t) = -—g 2 (t) + O(cubic) , (2.13) 

for all finite values of t, where all cubic and higher order terms in g and m are readily 
computed. We can substitute this result in Eqs. (2.10) and (2.11) and we see that the 
irrelevant variable has no effect to leading order in an expansion in powers of g. For very 
large, but finite values of t, we find 
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m(t) = -jg 2 (t) + 0(g 3 ) 



(2.14) 



Therefore, the marginal variable must satisfy 



^ = - C3 g* + 0(g 4 ), (2.15) 

for large but finite values of t. Eq. (2.15) can now be solved easily to obtain an accurate 
approximation for g{t) for large t, 



g*(t) = ^(0) (216) 

To obtain an accurate approximation for small t, we can continue to use Eq. (2.13), 
but we need to use an iterative algorithm to improve our estimate of g(t) and m(t). This 
is done by integrating Eqs. (2.10)-(2.11) near t = repeatedly, using the estimates from 
one iteration on the right-hand-sides of these equations to generate a subsequent estimate. 
This process is repeated until a convergence criterion is met. The initial seed is given by 
Eq. (2.13), the solution of Eq. (2.15) for g(t), and 



m(t) = m e- 2t + | (g 2 {t) - g 2 e' 2t ) , (2.17) 

for example. The transient behavior near t = in this approximation of m(t) is wrong, 
but any guess that is sufficiently near the solution should lead to convergence. After this 
iterative process converges, the desired result is obtained. In this case, the only output is 
w(0), because m(0) and g(0) are input. 

In the simplest case there is only one marginal variable and a finite number of rele- 
vant variables. It is assumed that all variables are small near the critical Gaussian fixed 
point, and in particular it is assumed that the marginal variable is small. There are an 
infinite number of irrelevant variables, but one can classify these variables according to 
the eigenvalues in Eq. (2.5) and according to their magnitude in terms of the marginal 
variable. One can first replace the irrelevant variables with an appropriate polynomial 
involving marginal and relevant variables. To leading order these are given by the zeroes 
of the polynomials on the right-hand-side of Eq. (2.9). Substituting these results in Eq. 
(2.7) one can next determine the strength of each relevant variable in terms of the single 
marginal variable. The appropriate polynomials are given to leading order by the zeroes of 
the polynomials on the right-hand-side of Eq. (2.7). After this one also has expansions for 
the irrelevant variables in terms of the single marginal variable. Every irrelevant variable 
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has a leading term of 0(g Pm ), where g is the single marginal variable. There are typically 
a finite number of irrelevant variables with an eigenvalue X m from Eq. (2.5) greater than 
any given value and with any given value of p m , and one can construct a first approxima- 
tion by keeping only the 'leading' irrelevant operators. In other words, one constructs a 
perturbative approximation for the trajectory in which the order of perturbation theory is 
determined by the single marginal variable g. Generalization of this procedure to the case 
where there is any finite number of marginal and relevant variables is straightforward. 

Next consider the case where there are an infinite number of relevant variables, in 
addition to an infinite number of irrelevant variables; but there are only a finite number 
of marginal variables. The irrelevant variables are treated exactly as above, except now 
the evolution of the irrelevant variables may be sensitive to an infinite number of relevant 
variables. There are an infinite number of boundary values that must be specified at / = Af 
in principle, and an infinite number of polynomials that must be considered in Eqs. (2.7) 
and (2.9). One can simultaneously classify all of the irrelevant and relevant operators in 
terms of the leading powers of the marginal variables that appear in the zeroes of these 
polynomials, just as the irrelevant operators were classified above. To leading order, one 
can replace the irrelevant and relevant variables in the equations for the marginal variables 
using these zeroes, and one can construct a perturbative approximation for the marginal 
variables for <C / <C Af by dropping sub-leading irrelevant and relevant variables. This 
part of the procedure is a straightforward generalization of the procedure used to handle 
an infinite number of irrelevant variables; however, there is a crucial difference between 
relevant and irrelevant variables. As discussed above, the boundary values chosen for the 
irrelevant variables are arbitrary. The boundary values chosen for the relevant variables do 
not affect the trajectory for / <C Af, because their effects are exponentially suppressed as 
/ decreases in Eq. (2.7); however, we need to construct Hjj, and the boundary values for 
the relevant variables may have important effects on all irrelevant, marginal, and relevant 
variables near / = Af. 

Remember that the boundary values of the relevant and marginal variables are input, 
while the values of the irrelevant variables at / = Af are output by the renormalization 
group equations. If there are an infinite number of relevant variables, we are forced to 
fix an infinite number of boundary conditions. Moreover, even if we want to compute 
a finite number of 'leading' irrelevant variables at / = Af, in principle we must approx- 
imate the evolution of an infinite number of relevant variables near / = Af, because all 
of the relevant variables affect the evolution of each irrelevant operator; and near / = Af 
transient behavior may prevent us from replacing the relevant variables with functions of 
the marginal variables. Does this render the perturbative renormalization group useless? 
Hopefully not, for two reasons. First, there are many problems in physics where an infi- 
nite number of boundary conditions must be fixed (e.g., the value of a field on a surface). 
There are also many problems in which the evolution of an infinite number of variables 
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must be computed (e.g., the value of a field in a volume). The key to the solution to 
such problems is to show that it is possible to approximate an infinite number of variables 
with a few well-chosen variables (e.g., parameters in functions of the original variables). 
In the case of light- front field theory we encounter an infinite number of relevant operators 
because functions of longitudinal momenta appear in these operators. However, a finite 
number of functions appear, and specifying the infinite number of boundary conditions is 
accomplished by specifying these functions. 

Further consideration of the renormalization group equations reveals a second reason 
that the perturbative renormalization group may survive, even though there are an infinite 
number of relevant operators. It may be natural for only a few of the relevant variables 
to be 'independent'. Suppose that we independently fix the boundary values of an infinite 
number of relevant variables and then solve Eqs. (2.7)-(2.9). Since we assume that Hj\f is 
near the Gaussian fixed point, all relevant variables are small at / = M, and we expect all 
relevant variables to approach polynomials of the marginal variables at an exponential rate. 
It is only possible for these variables to deviate significantly from these polynomials near the 
end of the trajectory, so on any trajectory of renormalized Hamiltonians near a Gaussian 
fixed point all relevant variables are exponentially near these polynomials over an infinite 
number of cutoff scales. Thus, it is natural to speculate that almost all of these variables 
track these polynomials exactly, never deviating from this behavior. If only a finite number 
of relevant variables depart from these polynomials, we can approximate the trajectory by 
dropping all sub-leading relevant variables and numerically computing the behavior only 
of those relevant variables that depart from these polynomials. I call a relevant variable 
whose evolution is exactly given by a polynomial a dependent relevant variable. A relevant 
variable whose value departs from such a polynomial is called an independent relevant 
variable. I simply postulate that there are only a finite number of independent relevant 
variables in theories of physical interest (e.g., Lorentz covariant theories). The boundary 
values of dependent relevant variables cannot be adjusted independently because they are 
determined by the polynomials; and this places severe constraints on the theories that 
satisfy this postulate. These are the coupling constant coherence conditions, and Wilson 
and I have shown that they arise naturally when there is a hidden symmetry [39]. 

If there are an infinite number of marginal variables, in addition to an infinite number 
of relevant and irrelevant variables, their behavior is unpredictable without further assump- 
tions. I postulate that there are a finite number of independent marginal variables, and 
that the infinite number of remaining dependent marginal variables can each be replaced 
by a polynomial expansion in powers of the independent marginal variables. 

To see that these postulates are reasonable, one must try to understand why light- 
front renormalization group transformations lead to an infinite number of relevant and 
marginal variables. We will see in Section III that every cutoff that is run by a light-front 
renormalization group transformation breaks Lorentz covariance, and most violate cluster 
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decomposition. Moreover, in gauge theories these cutoffs violate gauge invariance. The 
price one pays for violating a continuous symmetry in this case is the appearance of an 
infinite number of relevant and marginal variables. This price should not be surprising, 
because these symmetries must be restored to all physical observables and this requirement 
imposes an infinite number of conditions on the renormalized Hamiltonian. These condi- 
tions relate the strengths of operators in the Hamiltonian, providing an infinite number of 
relationships between the relevant and marginal operators at / = M. Thus, it is expected 
that the appearance of an infinite number of relevant and marginal variables does not imply 
that there are an infinite number of free parameters in the theory. There should be ex- 
actly as many independent relevant and marginal variables in a light-front renormalization 
group analysis as there are in a Euclidean renormalization group analysis, and the strength 
of all dependent relevant and marginal variables should be fixed by the requirement that 
the broken symmetries are restored. These symmetries can be restored by the adjustment 
of the strength of the dependent variables at the boundary, and the requirement that the 
symmetries be maintained for all other cutoffs places an infinite number of conditions on 
the renormalization group. These points are illustrated by examples in Sections VI and 
VII. 

In order to summarize the most important aspects of a perturbative renormalization 
group analysis and clarify the difference between perturbative and nonperturbative anal- 
yses, I introduce Wilson's triangle of renormalization, shown in figure 1. The triangle 
displays a sequence of renormalization group trajectories of increasing length, Af. We can 
label the Hamiltonian using a superscript to denote the absolute cutoff and a subscript to 
denote the effective cutoff, H^°. Assume that the cutoff is a cutoff on energy. The object 
of the renormalization group is to make it possible to let Ao — > oo while keeping A^ at 
some fixed value, say 2 GeV. The subscript M indicates how many times the transforma- 
tion must be applied to the original Hamiltonian to lower the effective cutoff to its final 
value, so one has Aj\f = Aq/(2^) for example. To fix Aj\f and let Ao — > oo, we must also let 
A/"— > oo. 

The renormalization group enables one to compute renormalized Hamiltonians, shown 
as the right- most column in figure 1, by providing an operational definition of renormal- 
ization. In the perturbative renormalization group this task is reduced to solving a finite 
number of difference or differential equations, as shown above. At each stage of a non- 
perturbative renormalization group calculation one selects a cutoff Hamiltonian H^°, and 

applies the transformation, T, M times to generate the Hamiltonian . In a successive 
stage one selects a new Hamiltonian and increments J\f by one. The sequence of initial 
Hamiltonians are related in a manner that must be determined as part of an algorithm 
tailored to the specific theory. In a nonperturbative calculation one probably must con- 
struct the triangle of Hamiltonians directly, being satisfied with numerical evidence that 
the limiting trajectory of renormalized Hamiltonians exists. 



24 



If H%° lies near a Gaussian fixed point we have seen that the irrelevant variables in 

should have an exponentially small effect on most of the trajectory. We have also 
seen that we want to fix the strength of some operators at the end of the trajectory, not 
at the beginning. Similar features should appear in a nonperturbative analysis, but there 
is no general procedure to identify and order operators in this case. We may have to 
simply search on initial values of parameters that are identified as important, to find 
that yields desired operators in H^. After convincing oneself that an arbitrarily long 
trajectory can be constructed if is adjusted to sufficient precision, in practice one 
hopefully needs to explicitly construct only the last part of a finite trajectory to obtain an 
accurate approximation of a renormalized Hamiltonian, H^. 

In light-front renormalization groups there are an infinite number of relevant and 
marginal operators, because undetermined functions of longitudinal momenta appear in 
a finite number of operators that are relevant or marginal according to transverse power 
counting. In a perturbative analysis one can use continuous symmetries such as Lorentz 
covariance to fix these functions in each order of an expansion in terms of a single coupling; 
however, in a nonperturbative analysis one may have to parameterize each function and 
seek an approximation of the renormalization group transformation that is represented by 
a set of coupled equations for the evolution of these parameters. 

Wilson has given an excellent discussion of the relationship between divergences in 
standard perturbation theory {e.g., Feynman perturbation theory) and the perturbative 
renormalization group [13], and I close this Section by repeating the most salient points. 
There are usually no divergences encountered when one applies the renormalization group 
transformation once; however, divergences can arise in the form of powers of / and expo- 
nents containing /, when T is applied a large number of times; and these divergences are 
directly related to the divergences in Feynman perturbation theory. There are no diver- 
gences apparent when one solves the perturbative renormalization group equations using 
a stable numerical algorithm; however, if one attempts to expand a coupling g\ in powers 
of go, for example, powers of I appear. As / — > oo these powers of / lead to the divergences 
familiar in Feynman perturbation theory. One can see an example of this by studying Eq. 
(2.16), where t represents the logarithm of a ratio of cutoffs. If the right-hand side of Eq. 
(2.16) is expanded in powers of g(0) : each term diverges like a power of t; so this expansion 
is numerically useless. On the other hand, g(t) is perfectly well-behaved if g(0) is small, 
and the divergences result from the fact that g(t) and g(0) differ by orders-of-magnitude 
for sufficiently large t. 

Eq. (2.16) illustrates the significant improvement over standard perturbation theory 
offered by the perturbative renormalization group. In standard perturbation theory one 
deals only with bare and renormalized parameters. This is analogous to dealing only with 
the parameters in and H^., without ever encountering separate parameters for every 
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Hamiltonian in the trajectory. Except in super-renormalizable theories, the ratio of bare 
and renormalized parameters goes to infinity (or zero). If a perturbative expansion of an 
observable in powers of the renormalized parameters converges, the expansion for the same 
observable in terms of the bare parameters cannot converge. This leads to some interesting 
departures from logic in standard perturbation theory [105]. A small contribution of the 
renormalization group is that logic may sometimes be restored to perturbation theory. 

III. Light-Front Hamiltonians and Renormalization Group Transformations 

The first step in defining a renormalization group transformation is to define the space 
of Hamiltonians upon which this operator acts. I give no precise definition of this space, 
partially because it must usually be defined after studying a transformation, not before. 
I restrict myself to scalar field theory, as it is straightforward but tedious to generalize 
the discussion to theories that include more complicated fields. I indicate what kind 
of operators are allowed in the Hamiltonians by example, and I display these operators 
in several forms that prove useful. Ultimately the Hamiltonians must be expressed in 
terms of Fock space eigenstates of the Gaussian fixed point Hamiltonian (i.e., in terms of 
projection operators) if one wants to use an invariant-mass transformation, so much of the 
early discussion is schematic. A brief summary of canonical light-front scalar field theory 
is given in Appendix A. 

The operators and constants with which Hamiltonians can be formed in a 3+1 dimen- 
sional scalar field theory, and their naive engineering dimensions, are 
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(3.1) 



I should also note that the Fourier transform of the field operator, 4>(q), has the dimension 
[x 1 -]. I work with a metric in which x ± = x° ± x 3 . In addition to derivative operators 
and the scalar field operator, I indicate that there may be a cutoff with the dimension of 
transverse momentum (A) that can be used, and there may be a cutoff with the dimension 
of longitudinal momentum (e) that can be used. All masses are expressed as dimensionless 
constants multiplying A. In general there may be many cutoffs (e.g., different cutoffs in 
different sectors of Fock space), but all of them can be expressed in terms of A and e. 

Perhaps the most important feature of Eq. (3.1) is that transverse and longitudinal 
dimensions are treated separately [106], just as one treats time and space differently in 
nonrelativistic physics. There is no analog of physical mass with the dimensions of longi- 
tudinal momentum instead of transverse momentum, because longitudinal boost invariance 
is a scale invariance [44,93], and physical masses (not necessarily bare masses) violate scale 
invariance. The cutoff e is the only constant with the dimensions of longitudinal momen- 
tum that can enter the definition of the Hamiltonian, and it must enter in a manner that 
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restores boost invariance to observables despite the violation of explicit boost invariance 
caused by the cutoff itself. In general one cannot be sure that naive engineering dimensions 
are significant in an interacting theory; however, near a Gaussian fixed point naive power 
counting is appropriate for the same reasons it is appropriate in standard perturbation 
theory. This is explicitly shown in Section V for light-front transformations. 

The assumption of transverse locality naively means that no inverse powers of d 1 - are 
allowed. Restrictions on inverse powers of d 1 - are clarified in Section VI where they first 
appear in the second-order behavior of the light-front transformations. 

Masses appear in the Hamiltonian as dimensionless constants multiplying A. I always 
assume that physical masses are much smaller than A, and I make an operational distinc- 
tion between physical masses and mass counterterms. Mass counterterms are present even 
when the physical mass is zero unless a symmetry protects the mass operator, and I know 
of no examples in cutoff light-front field theory where this occurs. Mass counterterms can 
have a very different dependence on longitudinal momenta than physical mass terms, as 
shown in Section VI. Except in Appendix C I usually assume that the physical mass is 
zero and focus on the critical theory. 

I make no initial restriction on the manner in which longitudinal derivatives appear. 
In canonical scalar field theory longitudinal derivatives appear only as inverse powers (see 
Appendix A); however, in canonical light- front QCD in light-cone gauge [107,101] one finds 
both inverse powers and powers of the longitudinal derivative in the three-gluon vertex 
and in the exchange of an instantaneous gluon between a quark and gluon or between two 
gluons. 

The Hamiltonian, H, is the integral of the Hamiltonian density, H = J dx^d 2 x^7i. 
Their dimensions are easily derived in canonical free field theory, and here I simply take 
them as given to be 
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(3.2) 



Given the catalog of operators from which the Hamiltonian can be formed, the space of 
Hamiltonians that I initially consider consists of all operators that can be formed from 
the basic catalog and that have the appropriate engineering dimension. Inverse powers 
of the transverse derivative operator are excluded initially and inverse powers of the field 
operator are always forbidden. Furthermore, cutoffs must be imposed to complete the 
definition. 

I work in momentum space rather than position space, and the Hamiltonian can be 
written schematically as 
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H = - dqi dq 2 (16tt )S (91 + 92) ^2(91* 92) 0(9i)0(92) 
+ 41 / d ^ d ^ 2 d ^ 3 dd - A ( 167I_3 ) 53 (^i + 92 + 93 + 94) 



(3 3) 
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+ oT J dqi " ' ( 1671 " 3 ^ 3 ^ 1 H h 56 ) u e(9i, 96) 0(9i) • • • 0(96) 



+ 



where, 



dq+d 2 q ± 

dq= m^+- (3 - 4) 



2 2 

We will see that this leads to a free energy (q^ + m?)/q + when U2 = + m 2 . In terms 
of plane wave creation and annihilation operators, 



0(g) =0(9) (9 + >0); 0(9) = ~a\-q) (q+ < 0) . (3.5) 



The above restrictions on transverse derivatives become restrictions on the functions 
U2,U4, etc. At this point the integrals include both positive and negative longitudinal 
momenta. The next step toward an expression for the Hamiltonian that can be directly 
manipulated involves replacing the field operators in Eq. (3.3) with creation and anni- 
hilation operators, normal-ordering the Hamiltonian and changing variables so that only 
positive longitudinal momenta appear. There are no modes with zero longitudinal mo- 
mentum. This complicates the Hamiltonian algebraically, but the advantages far outweigh 
this complication. I should mention that there is no need to define the normal-ordering 
operation until after the cutoffs required by the transformation are implemented, and that 
after this no divergences are encountered in the normal-ordering procedure. The initial 
Hamiltonian is simply assumed to be normal-ordered. The schematic Hamiltonian becomes 
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Since modes with identically zero longitudinal momentum are not allowed, there are 
no operators in this light-front Hamiltonian that contain only creation operators or only 
annihilation operators. This is only natural in light-front field theory. From this point for- 
ward it is assumed that all longitudinal momenta are positive. The Hamiltonian is simply 
assumed to be normal-ordered, and it is assumed as part of the definition of each trans- 
formation that the transformed Hamiltonian is normal-ordered. This is readily insured 
in perturbation theory, where one only encounters products of operators that are easily 
normal-ordered, as shown below. 

When one studies a boost-invariant renormalization group transformation, the Hamil- 
tonians must be written in terms of projection and transition operators that are constructed 
from the free many-body states generated by products of the above creation and annihila- 
tion operators acting on the vacuum. This introduces severe notational complications, but 
such representations are trivially constructed and manipulated. All boost-invariant cutoffs 
involve the total longitudinal and transverse momenta of a state, and not simply a few 
momenta carried by individual particles. When the cutoff depends on extensive quantities 
such as the total momentum, spectator-dependence inevitably appears in the operators, as 
we will find in Section VI. Such spectator-dependence has been studied recently in LFTD 
[63-71], but I want to emphasize that it is required even if one does not make a Tamm- 
Dancoff truncation. Since interactions depend on spectators, they also depend on the Fock 
space sector(s) in or between which they act [63,65,68]. In this case, one might worry that 
the distinction between the functions u 2 , u 4 , etc. might become altered or blurred. How- 
ever, these distinctions are easily maintained by considering how many individual particle 
momenta are altered by the operator. Examples below should clarify these points. 
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In order to develop simple examples of the light-front renormalization group one can 
use a Tamm-Dancoff truncation [60,61], which simply limits the number of particles and 
introduces an additional source of sector-dependence in the operators [63]. In LFTD the 
Hamiltonian again must be written in terms of projection operators [68]. Let me emphasize 
that some of the most interesting results one can find in a renormalization group analysis 
are apparently lost when a Tamm-Dancoff truncation is made, and I do not use LFTD in 
any examples in this paper. However, to illustrate what a Hamiltonian written in terms 
of projection and transition operators looks like, let me truncate Fock space to allow only 
the one-boson and two-boson sectors. In this case the complete Hamiltonian is 
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(3.7) 



Note that \ki,k 2 ) = a*(ki)aJ(k 2 )\0), etc. More complicated examples are readily con- 
structed. A superscript is added to u 2 , u 3 , etc. to indicate the Fock space sector (s) within 
which or between which the operator acts. In this example only the superscript on u 2 is 
required. The Hamiltonians that are actually manipulated in the light-front renormaliza- 
tion group are similar to this last expression, not the expressions in Eqs. (3.3) and (3.6). 
However, to see the connection with the more familiar canonical formalism, it is necessary 
to start with the latter expressions. 

As mentioned above, I have simply dropped all terms that explicitly involve zero modes 
only; i.e., terms containing creation operators only or annihilation operators only. The 
motivation for dropping such terms was briefly discussed in the Introduction, but here I 
simply assume this restriction is placed on the space of Hamiltonians ab initio . I do not 
question whether this assumption is reasonable and I do not believe that it is possible to 
answer such questions with a perturbative renormalization group analysis. The only issue 
at present is whether this restriction can be maintained or whether the transformations 
themselves regenerate pure creation or pure annihilation operators. To insure that this 
does not occur I simply drop the zero modes in every term in the Hamiltonian, and assume 
that this does not affect any momentum integral because a set of measure zero is being 
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subtracted. The zero modes are automatically removed by several of the cutoffs studied 
below, and there is no need to be careful about how the zero modes are removed in any 
of the examples considered in this paper. This may not be the case in QCD, where severe 
divergences associated with small longitudinal momenta force one to be more careful. I 
have nothing more to say about this issue until the Conclusion. 

Up to this point the operators act in an infinite dimensional Fock space unless otherwise 
specified, and momenta are left unconstrained. The definition of the space of Hamiltonians 
is not complete until the momentum cutoffs that enter the renormalization group transfor- 
mation are added. Let me begin by introducing the two simplest light-front renormalization 
group transformations, T 1 - and T + . 

In the renormalization group one starts with degrees of freedom that are already re- 
stricted by momentum cutoffs. There is never any point at which one explicitly considers 
Hamiltonians that contain no cutoffs, although renormalized Hamiltonians are obtained 
by considering what happens when the initial cutoff is taken to its limit. After introducing 
a cutoff one removes additional degrees of freedom and computes an effective Hamiltonian 
that is required, for example, to reproduce all of the low-lying eigenvalues and 'accurately' 
approximate the low-lying eigenstates. This step is radically different in the light-front 
renormalization group from the integration over large momentum components in a path 
integral [12] employed in a Euclidean renormalization group, and details are not discussed 
until Section IV. The final step in a light-front renormalization group transformation is to 
rescale all momenta so that their original numerical range is restored, to rescale the field 
operators, and to rescale the Hamiltonian itself. At this point the Hamiltonian is identical 
in form to the original Hamiltonian, but the functions U2, 114, etc. may all change. The 
aim of the analysis is to understand how these functions change. In general the analysis is 
most useful if one can reduce the task of following these functions to the infinitely simpler 
task of following a few constants. Naturally, the most difficult step is demonstrating that 
this simplification is a legitimate approximation. 

For the light-front renormalization group transformation T 1 - one begins by cutting off 
all transverse momenta, so that 

< p ±2 < Aq . (3.8) 
One then removes any state in which one or more momenta lie in the range 

A? < p^ 2 < A 2 . (3.9) 

Typically Ai = Aq/2 or Ai differs by an infinitesimal amount from Ao- In most examples 
in this paper I use Ai = Aq/2. A new Hamiltonian must be found that is able to reproduce 
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the eigenvalues and approximate the eigenstates of the original Hamiltonian without the 
degrees of freedom that have been removed. Techniques for computing such 'effective' 
Hamiltonians are discussed in the next Section. This new Hamiltonian must be written 
in the same form in which the initial Hamiltonian is written, and in particular it cannot 
include any explicit energy dependence. Let the initial Hamiltonian be called H^ , with 
the new effective Hamiltonian being . 

The renormalization group transformation T 1 - is completed by changing variables to 



p^fy, (3.io) 



scaling the field operator by a constant so that 



^(p 1 ,P + )=(V(p 1 'p + ), (3-11) 

and multiplying the entire Hamiltonian by a constant z^. The constants an d z 1 - 
are introduced so that fixed points can exist. They are part of the definition of the 
transformation and can be chosen freely. The only practical restriction on the choice of 
and z 1 - is that the choice should lead to a fixed point of physical interest, z 1 - is not 
essential, but it simplifies the task of comparing Hamiltonians, because it allows the range 
of eigenvalues of the transformed Hamiltonian to be the same as the range of eigenvalues 
of the initial Hamiltonian. The price paid for introducing z 1 - is that one must multiply the 
Hamiltonian by 1/z 1 - to obtain eigenvalues in the original units. If Af transformations are 
made, one must multiply the resultant Hamiltonian by (1/ ' z^)^ '. This point is clarified by 
examples below. The cutoffs in the original Hamiltonian introduce step functions into the 
momentum integrals in Eq. (3.3), for example. The elimination of the degrees of freedom 
with momenta in Eq. (3.9) leads to a new set of step functions, and the rescaling leads back 
to the original step functions. Each term in changes in a simple manner as a result 
of these various rescaling operations, and the final result is the transformed Hamiltonian, 
T [H]. The only difficult step is constructing H\ 1 . 

I do not discuss any examples in detail until Section V, where Gaussian fixed points 
and linear approximations of the transformations are constructed. However, to orient the 
reader I consider the simple Hamiltonian needed to find the Gaussian fixed point of T 1 in 
Section V. If there are no interactions the Hamiltonian is 



and when the cutoff changes the effective Hamiltonian is simply 
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ffA,= /«»(A?-q^( 4 )!i^. (3.13) 

The complicated procedure for producing effective Hamiltonians developed in Section IV is 
not needed when there are no interactions, because the original Hamiltonian projected onto 
the subspace of allowed states exactly reproduces all of the eigenvalues and eigenstates in 
this subspace. Next we change variables according to Eqs. (3.10) and (3.11), and multiply 
the Hamiltonian by a constant. The final transformed Hamiltonian is 

TVaJ= .V^Y^-q^Wfy.^)^^. (314) 

Naturally, when interactions are included the transformation is far more complicated. In 
fact, we will see that perturbation theory generally cannot be used to approximate T [H]. 

For the light-front renormalization group transformation T + , one begins by cutting off 
all longitudinal momenta, so that 

e <p + <oo. (3.15) 

It is immediately evident that this transformation is going to be radically different from 
the transformation on transverse momenta or any transformation considered in Euclidean 
field theory, because small momenta are removed instead of large momenta. This is an 
appropriate cutoff because the energy is (p 1 - + m 2 )/p + near the Gaussian fixed point of 
interest. States with small longitudinal momenta correspond to high energy states unless 
both p 1 and m vanish, leading to Eq. (3.15). 

The first step in the transformation is to remove all states in which one or more particle 
momenta lie in the range 

eo<P + <ei. (3.16) 

Typically e\ = 2eo or e\ differs by an infinitesimal amount from eo- Again a new Hamilto- 
nian must be found that reproduces the original low- lying eigenvalues and approximates 
the low-lying eigenstates of the original Hamiltonian without the degrees of freedom that 
have been removed. The procedures for computing such Hamiltonians are identical to 
those used for T^. However, because longitudinal boosts are scale transformations, as 
seen in Eq. (3.19) below, it is possible to say a great deal about the transformed Hamilto- 
nian without going through an explicit construction. In Section V I prove that physical 
Hamiltonians are fixed points of the transformation T + . 
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The light-front renormalization group transformation T + is completed by changing 
variables to 



P + ' = V , (3.17) 



scaling the field operator by a constant so that 



^/IKVIP 1 ,/), (3-18) 

and multiplying the entire Hamiltonian by a constant z + . Again the constants £ + and z + 
are introduced so that a fixed point can exist. I do not discuss T + further in this Section, 
but perturbation theory generally cannot be used to approximate T + either. 

At this point a general procedure for inventing light-front renormalization group trans- 
formations should be apparent, although many important details may be obscure. I close 
this Section by discussing boost-invariant light-front renormalization group transforma- 
tions. Both T 1 - and T + break manifest boost invariance because they employ cutoffs 
that are not boost-invariant. Under a longitudinal boost the Hamiltonian and all particle 
momenta are transformed [44,93], 



p+_ >e "p+ j p-L_> p -L ? H^e~ v H, (3.19) 
and under a transverse boost, 

P + ^P + , P 1 ^P 1 +p + v L , H -> H + 2P X • v x + V + v ±2 . (3.20) 

In Eq. (3.19) v is an arbitrary real number, and in Eq. (3.20) v 1 - is an arbitrary velocity, 
while V + is the total longitudinal momentum and P -1- is the total perpendicular momen- 
tum. Since the cutoffs are not changed by these transformations, these transformations 
place severe restrictions on the possible form of physical Hamiltonians that are not easily 
satisfied. Let me mention that this is not the chief problem one encounters when studying 
T 1 and T+. The chief problem, as mentioned above, is that no perturbative expansion 
exists for these transformations in general. When a perturbative expansion exists for a 
transformation, it should be possible to implement Lorentz covariance order-by-order in 
perturbation theory. 

Physical Hamiltonians must transform according to Eqs. (3.19) and (3.20), and if 
possible one would like to build these constraints into the space of Hamiltonians ab initio 
and construct transformations that automatically maintain these constraints. Boosts are 
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part of the kinematic subgroup of Poincare transformations in light-front field theory, and 
this should allow one to make the kinematic invariances manifest by choosing the correct 
variables. This kinematic subgroup is isomorphic to the two-dimensional Galilean group, 
and the use of appropriate variables resembles the separation of center-of-mass motion in 
nonrelativistic quantum mechanics [44,93]. The appropriate variables are 

x=|+, K ± = p 1 -xP ± . (3.21) 

One can easily verify that x and k 1 - are invariant under the above boosts. Therefore 
any cutoff constructed from these variables does not interfere with our ability to maintain 
manifest boost invariance. In particular, if we use a cutoff formed from these variables we 
may simply fix V + and choose P 1 - = as part of the definition of the space of Hamiltonians, 
without loss of generality. 

An obvious feature of such cutoffs is their use of the total momenta V + and P , 
which are themselves extensive quantities. The momentum available to a system in one 
room, which determines what free states can be superposed to form physical states, may 
depend on the state of a system in the next room; and this may introduce nonlocalities 
into the Hamiltonian when it is constrained to produce exact physical results despite such 
nonlocal constraints on phase space. This could be a very severe price to pay for manifest 
boost invariance, but it does not seem likely that one can avoid paying this price without 
sacrificing the possibility of using perturbation theory if cutoffs that remove degrees of 
freedom are used, as is shown in Section VI. The important point for a perturbative 
renormalization group analysis is to show that such cutoff nonlocalities introduce new 
irrelevant, marginal, and relevant operators that are readily computed in perturbation 
theory. It is essential to show that cutoff nonlocalities do not lead to inverse powers of 
transverse momenta that produce long-range interactions. We do encounter inverse powers 
of transverse momenta that do not produce long-range interactions because of the cutoffs 
that accompany them. 

There are many boost-invariant renormalization group transformations one can intro- 
duce, and I introduce three. All of these transformations are identical in form to T and 
T + , so in each case it is sufficient to specify the cutoffs. There is no reason to consider 
a cutoff on longitudinal momentum fractions alone, because if we change such a cutoff 
we cannot rescale momenta to recover their original range. When individual longitudinal 
momenta are rescaled, the total longitudinal momentum is rescaled by the same amount, 
and the longitudinal momentum fractions are invariant. In other words, longitudinal boost 
invariance is a scale invariance. The first cutoff one might consider starts by cutting off 
the relative transverse momenta defined in Eq. (3.21), so that 
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< K ±2 < An • 



(3.22) 



One then proceeds by lowering this cutoff, computing a new Hamiltonian, and then by 
rescaling the transverse momenta exactly as in Eq. (3.10). This cutoff should violate 
locality in a minimal fashion, but we shall see that perturbation theory again cannot be 
used to approximate this transformation. 

The second boost-invariant transformation begins with the introduction of the cutoff 



By expanding this sum one easily demonstrates that this cutoff acts only on the relative 
transverse momenta defined in Eq. (3.21). This cutoff is the invariant-mass cutoff for 
massless theories, and again the remaining steps in the transformation involve lowering 
the cutoff, computing a new Hamiltonian and rescaling transverse momenta and fields. 
We will find that this transformation can be approximated perturbatively if there are 
no physical mass terms in v>2- A variety of 'mass' terms arise in light-front field theory, 
because any function of longitudinal momentum fractions can accompany a mass. Physical 
mass terms appear in U2 with no additional function of longitudinal momenta; e.g., U2 = 

2 2 

p- 1 + m 2 results in a free energy of the form (p 1 - + m 2 )/p + . We will find in Section 
VI that the invariant-mass transformation leads to mass counterterms that appear in U2 
in the form # 2 Aq p + /P + , where g is a coupling constant and P + is the total longitudinal 
momentum of a state. This type of mass counterterm can be treated perturbatively, 
whereas a perturbative treatment of the physical mass terms leads to divergent longitudinal 
momentum integrals. 

The only transformation that I have been able to approximate using perturbation 
theory when there are physical mass terms in U2 begins with the cutoff 



Here the cutoff masses rrii that appear must be specified, and the appropriate values can 
only be chosen after an analysis of the transformation. In Wilson's triangle of renormaliza- 
tion discussed in Section II, one must consider the limit in which the initial cutoff An — > oo. 
When taking this limit one should consider the masses, m^, to be fixed. Clearly, if this 
limit can actually be taken, the specific values of the cutoff masses should enter primarily 
in the justification of perturbation theory and are adjusted for that purpose. 




(3.23) 




(3.24) 
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I go through each step of this transformation, because cutoff masses introduce signifi- 
cant new features. The first step in the transformation is to remove all states in which 



P ±2 + A 2 ^ p-L? + m? P ±2 + A 2 



< 



£^i<^Pa. (3.25) 



V+ ^ p+ v+ 

After computing the effective Hamiltonian one rescales all transverse momenta by changing 
variables as in Eq. (3.10). However, after this change of variables the states satisfy the 
constraint 



P ±2 i + (Ap/Ai) 2 m?^ P ±2 + A 2 



< " ' < • (3.26) 



The momenta do not have the same range after rescaling as they initially had, because 
the cutoff masses change. We will see that masses in the Hamiltonian rescale in exactly this 
manner. If one applies the transformation a large number of times, eventually the factor 
(Ao/Ai) 2 m? becomes large and all of phase space is eliminated by the transformation. 
Clearly the transformation must be highly nonperturbative at this point; however, there 
is no reason to believe that one can ever lower the cutoff on transverse momenta to a scale 
comparable to physical mass scales in the problem without the transformation becoming 
highly nonperturbative. One should consider the cutoff masses to be of the same scale as 
physical masses. This transformation still defines a semi-group, and I show in Appendix 
C that all of the basic features of the renormalization group apparently survive when 
one allows such a transformation. In particular, one can define the functions 112, 114, 
etc. independently of the cutoffs, and study their evolution as analytic functions of their 
arguments over the original range of momenta. This discussion is clarified by examples in 
Appendix C. 



IV. Perturbation Formulae for the Effective Hamiltonian 

I discuss two related methods for computing the effective Hamiltonian after the cutoff 
is changed, one developed originally by Bloch and Horowitz [90] and a second developed 
by Bloch [91]. I call the first Hamiltonian a Bloch- Horowitz Hamiltonian and the second 
simply a Bloch Hamiltonian. The Bloch Hamiltonian was used by Wilson in early work on 
the renormalization group [5] . More sophisticated methods must be employed if one wants 
to work far from the Gaussian fixed point, but I concentrate only on the development of 
the effective Hamiltonian in perturbation theory. 
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The primary requirement is that the effective Hamiltonian produce the same low-lying 
eigenvalues in the second cutoff subspace that the initial Hamiltonian produces in the initial 
cutoff subspace. A second requirement is that the eigenstates of the effective Hamiltonian 
be orthonormalized projections of the original eigenstates. The Bloch Hamiltonian is de- 
signed to satisfy these requirements, while the Bloch-Horowitz Hamiltonian is not. One of 
the primary reasons that one wants to preserve the eigenstates in addition to the eigen- 
values is to compute measurable quantities in addition to the energy. Each measurable 
quantity corresponds to the absolute value of a matrix element, and when one introduces 
cutoffs all operators must be properly renormalized so that their matrix elements in the 
cutoff space are independent of the cutoffs. An intermediate step in the construction of 
the Bloch Hamiltonian is the construction of an operator 1Z. 1Z can be used to renormalize 
all observables. 

I add two additional requirements that must be met by the effective Hamiltonian. 
First, it must be manifestly Hermitian. Second, its representation can only involve the 
unperturbed (i.e., free) energies in any denominators that occur. The second property is a 
severe limitation that limits the utility of the Hamiltonian to studies of Gaussian and near- 
Gaussian fixed points. This property is not desirable a priori and is adopted only because 
I do not know of a general procedure for constructing Hamiltonians outside of perturbation 
theory. This limitation does not imply that one is always limited to the use of perturbation 
theory, after one uses a perturbative construction of the light-front renormalization group 
Hamiltonian. The perturbative construction only needs to be valid near the Gaussian 
fixed point. After one has used the transformation many times to reduce the cutoff to 
an acceptable value, the final Hamiltonian may have been accurately computed using 
a perturbative renormalization group transformation, but its diagonalization is typically 
nonperturbative. I am particularly interested in QCD, where one can hopefully use a light- 
front renormalization group to remove high energy partons through sequential application 
of a perturbative renormalization group transformation, justified by asymptotic freedom 
[108,109]; after which one must employ suitable nonperturbative techniques to diagonalize 
the final cutoff Hamiltonian and obtain low energy observables. 

I begin with the Bloch-Horowitz Hamiltonian because it is most easily derived. In fact, 
the primary reason I include a discussion of the Bloch-Horowitz Hamiltonian is because its 
development is particularly transparent and allows a reader unfamiliar with such many- 
body techniques to readily grasp the main ideas. Let the operator that projects onto all 
of the states removed when the cutoff is lowered be <2, and let V = 1 — Q. Of course, 
Q 2 = Q and V 2 = V. I occasionally refer to a subspace using the appropriate projector. 
With these projectors, Schrodinger's equation can be divided into two parts, 

VHV\^>) + VHQ\^) = EV\V) , (4.1) 
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QHV\V) + QHQ\^) = EQ\V) . (4.2) 

Solving Eq. (4.2) for Q\*&) and substituting the result into Eq. (4.1), one finds an operator, 
H 1 , that produces the same eigenvalue E in the subspace V as the original Hamiltonian, 

H\E)=VHV + VHQ—^——QHV. (4.3) 

h, — SzitL y 

This is the Bloch-Horowitz effective Hamiltonian. It is not satisfactory, however, because 
it contains the eigenvalue, E; and it is not Hermitian unless E is held fixed. 

The development of the Bloch Hamiltonian begins with the definition of the same 
projection operators Q and V used above. After defining these operators one looks for an 
operator TZ that satisfies 

Q\V)=KV\V) , (4.4) 

for all eigenstates of the Hamiltonian that have support in the subspace V. To construct 
this operator, multiply Eq. (4.1) by TZ and replace Q|^) in Eq. (4.2) with TZV\^), and 
subtract the resultant equations to obtain 

TZHpp - HqqTZ + TZHpgTZ - Hq P = . (4.5) 

Here I have introduced the notation Hpp = VHV, etc. This is the fundamental equation 
for TZ and one of the most difficult steps in constructing the Bloch Hamiltonian is solving 
this equation. I am only interested in the perturbative solution, and one can already see 
that TZ should be proportional to HQ-p. 

For notational convenience I let H = h + v for the Bloch Hamiltonian, where h is a 
'free' Hamiltonian and I assume [h, Q] = 0. Since VQ = 0, Eq. (4.5) can also be written as 

IZhpp — IiqqIZ — vqp + IZvpp — vqqTZ + IZvpqTZ = . (4.6) 

This equation is now in a form that can be solved using an expansion in powers of v, and 
it is apparent that TZ starts at first order in v. Before solving this equation, let me write 
the effective Hamiltonian and the eigenstates in terms of TZ. 

The states V\^) are not orthonormal. However, if we assume that no two eigenstates 
of H have the same projection in the subspace V (i.e., that TZ is single- valued) , then one 
can readily check that 
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H'\$)=E\&) , 



(4.7) 



when, 

1$) = y^i + n^n v |*) , (4.8) 

and, 

H'= - 1 =(V + rt)H(V + n)^=L==. (4.9) 

The states |$) are orthonormalized projections of the original eigenstates |^), and the 
manifestly Hermitian effective Hamiltonian H' yields the same eigenvalues as the original 
Hamiltonian H. H' is the Bloch Hamiltonian. 

To construct the Bloch Hamiltonian in perturbation theory, one first solves Eq. (4.6) 
in perturbation theory. Since we cannot assume that [1Z, h-p-p] or [1Z, Hqq\ are zero, we 
must employ the eigenstates 



hj>p\a) = e a \a) , h,QQ\i) = €i\i) , (4.10) 

to develop algebraic equations for the matrix elements of 71. I use \a), \b) , ... to indicate 
free eigenstates in V; and \i), ... to indicate free eigenstates in Q. In order to expand 
H' through third order in the interaction, we only need the first two terms in an expansion 
of 1Z, and these are 



(z |72.i | a) = 



(i\v \a) 



(4.11) 



( ., |fl) = mmv\a) _ y, <?\v\b)(b\v\a) 

There are a few general features of each term in this expansion that I want to note. 
First, note that TZ only has nonzero matrix elements between a bra in Q and a ket in 
V. Second, every energy denominator involves a difference between a free energy in V 
and a free energy in Q. Ultimately the convergence of this expansion is going to rest 
not only on the weakness of v, but also on the fact that high energy states are being 
eliminated so that these denominators are large throughout the most important regions 
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of phase space. We are ultimately interested only in the very low energy eigenstates and 
eigenvalues, and almost all of the states we eliminate to obtain the effective Hamiltonian 
for these states are extremely far off shell. The energy denominators may vanish, but this 
should only happen for a set of measure zero if the expansion is to converge. One can 
avoid all but accidental degeneracies by putting the system in a box, but I do not need to 
do this for the examples considered. Serious problems should only arise when this entire 
series needs to be re-summed because a particular nonperturbative effect must be properly 
included. Unfortunately, such serious problems are common. Higher order terms are easily 
constructed using the recursion relation 



n-2 

Tl n h vv - h QQ K n + K n -iv vv - v QQ K n -i + 1Z m v V QJl n - m -i = . (4.13) 

m—l 

Clearly the number of terms in each successive order grows exponentially, and it is likely 
that the expansion is at best asymptotic. 

Given a perturbative expansion for 1Z, one can use Eq. (4.9) to derive a perturbative 
expansion for the Bloch Hamiltonian. After some simple algebra, one finds that through 
third order in v, 



i \tti\i\ i i7 i7\ 1 / (a\v\i) (i\v\b) (a\v\i) (i\v\b) 
(a\H'\b) ={a\h + v\b) + -J2( LU1M + Ml I, 



1 (a\v\i)(i\v\j)(j\v\b) (a\v\i)(i\v\j)(j\v\b) > 

(e a -e t )(e a -e 3 ) + (e b - Ci)(e 6 - I ( - U ' > 

'•>3 



£ WWFAffl , (a\v\c)(c\v\i)(i\v\b) \ 4 
( eb - ei )(e c - ei ) + (e -ei)(e c -Ci) ' : ' ' 



This expression can be used to compute the renormalization of the quark-gluon vertex 
through third order in the bare coupling, for example [59]. It is quite similar to expres- 
sions one encounters in time-ordered perturbation theory for off-shell Green's functions, 
but there are some distinct differences. Given Eq. (4.14), or any extension of higher order, 
it is straightforward to derive a set of diagrammatic rules that allow one to summarize 
the operator algebra in fairly simple diagrams. These diagrams are similar to Goldstone 
diagrams [110], familiar from many-body quantum mechanics, but they require one to 
display the energy denominators. I refer to them as Hamiltonian diagrams. In a standard 
time-ordered diagram, there is a factor !/(£' — e^) for every intermediate state, and 
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the energy E is the same in every denominator. In the Hamiltonian diagrams there are a 
wide variety of denominators, always involving differences of eigenvalues of Hq; however, 
these eigenvalues may be associated with widely separated states in the Hamiltonian di- 
agram. While the above process of generating a perturbative expansion for the effective 
Hamiltonian is easily mechanized, I have not found a simple set of diagrammatic rules that 
summarize this process to all orders. It is quite possible that such rules already exist in 
the literature, but for the purposes of this paper it is only necessary to understand Eq. 
(4.14) and to appreciate the fact that higher order terms are readily generated. 

In figure 2 I show a few typical Hamiltonian diagrams that occur when H contains a 
interaction. Energy denominators are denoted by lines connecting the relevant states, 
with the arrow pointing toward the state whose eigenvalue occurs last in Eq. (4.14). 
One can infer from these lines which states lie in V and which states lie in Q, because 
energy denominators always involve differences of energies in V and energies in Q and the 
arrow points toward a state in Q. External states always lie in V of course. In a realistic 
calculation there are many different vertices, corresponding to the many different functions 
«4, uq, etc. in the Hamiltonian. Sometimes U2 determines h; however, in perturbation 
theory one usually needs to include part of 112 in v. This issue is clarified in the next 
Section. 

The division of the Hamiltonian into a 'free' and 'interaction' part above is arbitrary, 
except for the requirement that [h, Q] = 0. Let us suppose that this requirement is met, 
but some of the eigenvalues of h in V are larger than some of the eigenvalues of h in Q. 
The energy denominators in Eq. (4.14) pass through zero when this happens and one must 
question the convergence of the expansion in Eq. (4.14). The fact that the denominators 
may vanish does not automatically imply that the expansion does not converge, but at the 
minimum it forces one to carefully consider the boundary conditions required to construct 
the Green's functions in Eq. (4.14). However, even if the expansion in Eq. (4.14) does 
converge, it may prove useless for a renormalization group study. In general, each term 
in Eq. (4.14) leads to an infinite number of operators in the transformed Hamiltonian, 
because of the potentially complicated dependence each term may have on the momenta 
of the 'incoming' and 'outgoing' states. We will see in the next Section that when some 
of the eigenvalues of h in V are larger than some of the eigenvalues of h in Q, operators 
that are classified as irrelevant in the linearized renormalization group analysis occur with 
large and sometimes infinite coefficients that render this classification scheme useless. 

V. Critical Gaussian Fixed Points and Linearized Behavior Near Critical Gaus- 
sian Fixed Points 

A fixed point is defined to be any Hamiltonian that satisfies 
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T[H*] = H 



(5.1) 



The fixed point is central to the renormalization group analysis, and unless one has found 
a fixed point it is unlikely that anything can be accomplished with perturbation theory. 
The simplest example is the Gaussian fixed point, one for which = 0, uq = 0, etc. In 
this case the transformation that generates is trivial. Let me begin by discussing the 
Gaussian fixed point for the transformation T 1 . To find the Gaussian fixed point we need 
to consider 

ff A „= J d4 9(M- q >- i )v a ( q )*p*>, (5.2) 
and when the cutoff changes the effective Hamiltonian is simply 

H Al = y <U»(A? - q^ 2 ) M«) ■ (5-3) 

Next we change variables according to Eqs. (3.10) and (3.11), and multiply the Hamilto- 
nian by a constant. The final transformed Hamiltonian is 

T[H Aa] = *V a (£)7 " M^f) • 

The factor of (Ai/Ao) 2 in front of the integral arises because the Hamiltonian itself 
has the dimension given in Eq. (3.2). I absorb this overall rescaling by setting 

^ -(£)'• <"> 

z 1 - is introduced so that the transformed Hamiltonian can be directly compared with the 
initial Hamiltonian. To obtain eigenvalues in the initial units, one must multiply by l/z^. 
With z 1 - determined, the Gaussian fixed point is found by insisting that 

«2(q^ + ) = C ±2 «2(fV,<7 + ). (5.6) 

The general solution to this equation is familiar from Euclidean field theory [12], being a 
monomial in q^~. The general solution is 
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_L\n 



(5.7) 



One must decide ab initio which Gaussian fixed point is to be studied by fixing the disper- 
sion relation in the free theory, and I am only interested in the case n = 2. Since there is 
no mass term in this fixed point, the correlation length is infinite and it is called a critical 
fixed point. 

The fixed point of T 1 - contains an undetermined function of longitudinal momentum, 
f(q + ). This has important implications for the light-front renormalization group, which 
will become clearer below. 

One can look for fixed points that contain interactions, but no such fixed point has 
been found for a scalar field theory in 3 + 1 dimensions. It is possible to change the 
number of transverse dimensions [10] to 2 — e and construct an analog of the analysis 
found, for example, in the review article of Wilson and Kogut [12]. I do not pursue this 
idea, and turn to the next step in the analysis of T" 1 , the construction of the linearized 
transformation about the Gaussian fixed point. To construct the linearized transformation, 
consider Hamiltonians that are almost Gaussian, 



(5.8) 



H l = H* + 5Hi . (5.9) 

Here the subscript / labels the number of times the transformation has been applied, so 
that a sequence of Hamiltonians can be studied. I assume that 5Hi is 'small'. Applying 
the full transformation we find that 



6H l+1 = T[H* + 6Hi\-H 

= L H * ■ 5Hi + 0{8Hf) 



2 . (5.10) 



This equation defines the linear operator L#*, which depends explicitly on the fixed point. 
I typically drop the subscript on Ljj* and simply refer to L. 

In general the construction of the operator L is complicated; however, when H* is a 
Gaussian fixed point, L is easily constructed. In the first step of a renormalization group 
transformation a cutoff is changed and degrees of freedom are removed. Effective inter- 
actions result when an incoming state experiences an interaction, so that some momenta 
are altered and fall into the range being removed. However, both the incoming and out- 
going states must fall in the sector retained by the transformation, so a second interaction 
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is always required to return the state to the allowed sector. The only exception to this 
occurs when zero modes are allowed and one can form a 'loop' with a single interaction. I 
drop zero modes, so I can ignore this possibility and conclude that near the Gaussian fixed 
point all new interactions and changes to initial interactions in the effective Hamiltonian 
are 0(5H 2 ). This means that the only part of the full transformation that affects the 
linear operator is the rescaling. 

We immediately conclude that L for the Gaussian fixed point of T 1 is given by the 
rescaling operations in Eqs. (3.10) and (3.11), along with the overall multiplicative trans- 
formation of the Hamiltonian. z 1 - is given by Eq. (5.5), and C 1 - is given by Eq. (5.8), so 
one readily finds that 



/AaM) Ai ± I (5 - 11} 

Subsequent analysis of the full transformation should be simplified by identifying the 
eigenoperators and eigenvalues of this linearized transformation. We seek solutions to the 
equation 



L-0 = \0, (5.12) 

and the solutions are readily found to be homogeneous polynomials of transverse momenta. 
We can label O with a subscript that displays the number of field operators in O and a 
superscript that displays the number of powers of transverse momenta. For example, 

0^/d 5 -«(A§- q ^)M 4 ^^M. (5.13) 

In general more labels are required because there may be many possible polynomials of 
the same degree; however, the eigenvalue is determined by these two labels. Applying L 
to any operator O™ one finds 

/ Ai \ (m+n-4) 

L-0™= (^) . (5.14) 

Inspection of this result reveals that the eigenvalue is determined by the naive transverse 
engineering dimension of the operator. This is not surprising since the linear approximation 
of the transformation is given by the rescaling operations alone. 
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Relevant operators are defined to be operators for which A > 1; and since Ao > Ai, 
relevant operators have m + n — 4 < 0. Transverse locality implies that m > 0, and we 
always require n > 2, so the only relevant operator is one satisfying m = 0, n = 2, if we 
assume that a <p — > —<p symmetry is maintained. This is a mass term. There is no a 
priori reason to assume that such a symmetry can be maintained, and if it is violated a 
n = 3 relevant operator will appear. The m = 1, n = 2 operator is ruled out by rotational 
invariance about the z-axis. If one includes a power of a transverse momentum, the Lorentz 
indices on this momentum must be contracted with a second Lorentz index. In the absence 
of spin the only transverse Lorentz indices are carried by transverse momenta, and therefore 
single powers of transverse momenta cannot occur. Note that the dependence of O™ on 
longitudinal momenta does not affect this analysis; and in this sense there are an infinite 
number of relevant operators if there is one. 

Marginal operators are defined to be operators for which A = 1. Allowing — > — cj) 
symmetry to be broken, there appear to be (n = 2,m = 2), (n = 3,m = 1), and (n = 
4, m = 0) marginal operators. The (n = 3,m = 1) operator is excluded by rotational 
symmetry about the z-axis. Again, any dependence on longitudinal momenta does not 
affect this analysis. 

Irrelevant operators are defined to be operators for which A < 1, and one sees that al- 
most all operators are irrelevant by this definition. If inverse powers of transverse momenta 
are allowed, this conclusion is destroyed, and there are an infinite number of relevant and 
marginal operators in addition to those discussed above. 

This analysis is itself relevant only if there is a region near the fixed point in which 
the linearized transformation is a reasonable approximation of the full transformation. In 
this case there may be a convergent or asymptotic perturbative expansion for the trans- 
formation in a region near the fixed point. A perturbative analysis is probably useful 
only for theories that display asymptotic freedom. In this paper I simply assume that all 
couplings are small, even though we know that this condition cannot be satisfied for an 
interacting symmetric scalar theory if the initial cutoff approaches infinity, because it is 
not asymptotically free. 

Before discussing T + , I would like to explicitly show what happens when a physical 
mass is added to u\- In this case we have 

ffA„ = / ^«(A§ - q^) + m 2 ] . (5-15) 

One can readily verify that after T 1 - is applied M times the result is 

T^[H Ao ] = / ^(A?, - q^ 2 ) [q^ 2 + *V] ■ < 5 - 16 > 
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where I assume Ai = Aq/2. As N — > oo the Hamiltonian approaches another Gaussian fixed 
point, the so-called trivial fixed point at which the mass is infinitely large in comparison 
to the kinetic energy. Note that the original energy spectrum can always be obtained by 
multiplying the Hamiltonian by '. 

The analysis of the Gaussian fixed point and linearized approximation for T + closely 
parallels that of T 1 . To find the Gaussian fixed point we must analyze 

H eo = Jdqe( q +-e Q )u 2 ( q )^^-, (5.17) 
and when the cutoff changes the effective Hamiltonian is 

H„= J JWf (5.18) 

Next we change variables according to Eqs. (3.17) and (3.18), and multiply by an overall 
constant to obtain 



T[H n ] = ,+ C + 2 p) f - eo) urfqS V) =^ . (5-19) 

J eo q^ 

Here the factor of (eo/ei) in front of the integral arises because the Hamiltonian has 
the dimension given in Eq. (3.2), and I absorb this overall rescaling by setting 

z+ = (^) . (5.20) 
Given z + , the Gaussian fixed point satisfies 

^(q ± ,g + ) = C +2 «2(q ± , V)- (5-21) 
The general solution is a monomial in q + , 

«2(q ± ,? + ) = /(q ± )(^) n , (5-22) 

t + =(-Y- (5-23) 
Veo/ 

Again, one must decide ab initio which Gaussian fixed point to investigate, and I am 
interested in the case n = 0. 
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T 1 - and T + both have Gaussian fixed volumes rather than fixed points, because neither 
scales all momentum components simultaneously. The single fixed point that both share, 
with the definitions of the scaling constants above, is 



«2(q^ + )=q ±2 . (5.24) 

The linearized approximation for T + is determined by the rescaling operations in Eqs. 
(3.17) and (3.18). Combined with Eq. (5.23) these yield 

u n (q\, q+, ...) -> ^+C +n (-) UnCq-S, •••) 
^u n {q i,— ?i ,...) • 

The eigenoperators of this transformation are homogeneous polynomials of inverse longitu- 
dinal momenta. I choose to employ inverse derivatives for reasons that will become clear. 
For T + the solutions need not be labelled by the number of field operators, because the 
scalar field is longitudinally dimensionless, as in Eq.(3.1). However, for convenience I use 
the same notation found in Eq. (5.13), but with the superscript now indicating the power 
of inverse longitudinal momenta. In this case one finds that 

Oj= /W- £ o)[^l 4 ^^, (5.26) 



q 



q 



and 



L '°" = (|) m °"- (5 - 27) 

I have not included the inverse power of q + found in the measure dq, nor the inverse power 
of q + that accompanies any product of creation and annihilation operators found in the 
Hamiltonian, in the index m. This is clear in Eq. (5.26). The eigenvalue in this case is 
determined by the naive longitudinal engineering dimension of the operator. 

eo < 6\, so longitudinal relevant operators must have m < 0; longitudinal marginal 
operators have m = and irrelevant operators have m > 0. This explains the choice of 
inverse derivatives, because powers of inverse derivatives produce irrelevant operators. The 
same type of simplification that occurs in a Euclidean renormalization group when locality 
is assumed, might occur when one formulates a principle of longitudinal nonlocality for 
the longitudinal transformation. When we study second-order corrections in Section VI, 
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we will find that T + leads to problems that prevent us from exploiting any simplifications 
from longitudinal nonlocality, unfortunately. 

This completes the analysis of T + to leading order near the fixed point. The analysis 
of the next order shows that T + cannot be approximated in perturbation theory unless one 
allows arbitrary functions of transverse momenta, including functions that violate trans- 
verse locality. Worse, one encounters arbitrarily large coefficients of irrelevant operators 
in perturbation theory, so these operators cannot be dropped. 

It is fairly easy to see that Hamiltonians of physical interest are fixed points of T+, 
using Eq. (3.19). I demonstrate this and discuss some implications before closing this 
Section. In the first step of the transformation T + , we must remove degrees of freedom 
following Eq. (3.16), and find a new Hamiltonian that yields the same eigenvalues and 
properly renormalized eigenstate projections. A Hamiltonian that meets these criteria is 
readily found by applying a longitudinal boost, so that 



p+ ^ fi p + . (5.28) 

After this boost all longitudinal momenta satisfy the constraint e\ < p + < oo. Under this 
same boost we know that 



H^—H, (5.29) 

which reveals that when we multiply the boosted Hamiltonian by the constant (ei/eo), 
we retrieve the original Hamiltonian, but with momenta restricted to a new range, 
is completed by rescaling the variables according to Eqs. (3.17) and (3.18), and with 
z + = (ei/eo) and ( + = 1, this merely returns the longitudinal momenta to their original 
values before the boost and has no other effect on the Hamiltonian. In other words, as 
long as Eq. (5.29) holds, we have 



T + [H]=H, (5.30) 

so that H is a fixed point of the full transformation T+. 

This may be an important result. However, fixed points of T + that contain interac- 
tions invariably contain operators that violate transverse locality and leave one with no 
perturbative means to study the dependence of the Hamiltonian on the transverse scale. 
At this point I have discovered no benefit to implementing longitudinal boost-invariance 
by seeking fixed points of T+. 
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All of the boost-invariant renormalization group transformations have the same Gaus- 
sian fixed points as T^, and the same linear analysis as T^. Therefore, for each of these 
transformations, operators are classified according to their transverse dimensionality. Lon- 
gitudinal dimensionality plays no role whatsoever, and the linear analysis cannot be used 
to control the longitudinal functions in the Hamiltonian. There is no need to repeat the 
discussion leading up to Eq. (5.14). The only caveat is that strictly speaking the invariant- 
mass transformation that employs the cutoff in Eq. (3.24) has no fixed points. As discussed 
in Section III, after a finite number of transformations all of phase space is eliminated and 
the transformed Hamiltonian approaches zero, because of the masses in this cutoff; un- 
less the initial cutoff is allowed to approach infinity. However, using analytic continuation 
one can extend the functions 112, 114, etc. outside the range of the invariant mass cutoff 
and define the fixed point using the analytically continued Hamiltonian. In this case, the 
Gaussian fixed points are again identical to the Gaussian fixed points of T^. 

VI. Second-Order Behavior Near Critical Gaussian Fixed Points 

The perturbative renormalization group analysis becomes much more interesting when 
second-order corrections to the Hamiltonian are included. In fact, almost all features of a 
complete perturbative analysis are displayed in a second-order analysis. In this Section I 
begin by adding one simple interaction to the Hamiltonian and studying the interactions 
it generates to second-order. This analysis demonstrates that most candidate transforma- 
tions produce irrelevant interactions with divergent coefficients. I next concentrate on the 
invariant-mass transformation that runs the cutoff in Eq. (3.23). This transformation is 
used throughout the rest of the paper, except in Appendix C, where the effects of physical 
masses are briefly studied and a transformation that runs the cutoff in Eq. (3.24) must 
be used. It is shown that if one drops all irrelevant operators, a complete sequence of 
second-order transformations may be summarized by a few coupled difference equations. 
These equations are solved, and the physics implied by their solution is discussed. 

The final subject in this Section is the most important. I first show that when one 
uses the canonical </> 4 interaction and the invariant-mass cutoff, the boson mass shift and 
the boson-boson scattering amplitude violate both Lorentz covariance and cluster decom- 
position. The calculation of these expectation values reveals the 'counterterms' required 
to restore Lorentz covariance and cluster decomposition; and these counterterms include 
an infinite number of irrelevant operators, as well as a marginal operator that contains 
logarithms of longitudinal momentum fractions. This result is not surprising, and the 
important question is how to find a Hamiltonian that restores these properties to the 
mass shift and the scattering amplitude. Of course, one method is to simply compute 
amplitudes in perturbation theory, and add to the Hamiltonian the difference between the 
result desired and the result obtained. This is not very satisfactory, because it requires 
reference to a separate calculation, so I ask a question that Wilson and I have recently 
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posed [39]. What happens if one allows the Hamiltonian to become arbitrarily compli- 
cated, but one insists that a single coupling is allowed to explicitly depend on the cutoff, 
with all other interactions depending on the cutoff only through their dependence on this 
coupling? The answer may depend on the coupling one allows to run, and I select the 4 
canonical coupling, which a manifestly covariant analysis would show is the only coupling 
that must be allowed to run if — > — symmetry is maintained. In this Section I show that 
coupling coherence uniquely fixes all relevant and irrelevant operators to second-order in 
the </> 4 coupling constant; and that these operators restore Lorentz covariance and cluster 
decomposition to the boson mass shift, which is a relevant operator, and to the irrelevant 
part of the boson-boson scattering amplitude. I then show that a third-order analysis is 
required to determine the marginal operators to second-order in the canonical coupling. 
The third-order analysis is completed in Section VII. 

The second term on the right-hand side of Eq. (4.14) can be used to determine the 
second-order behavior of any transformation near a critical Gaussian fixed point. The fixed 
point Hamiltonian is h in Eq. (4.14), and all deviations from the fixed point are part of v. 
I should note that deviations of U2 from Ur> are produced by the second-order part of the 
transformation, but they do not directly affect the subsequent second-order behavior of the 
transformations. This is because we need two interactions to first produce an intermediate 
state with an energy above the new cutoff and then produce a final state with an energy 
below the new cutoff, and U2 does not change the state. 

If we consider the second-order behavior of a transformation when acting on a Hamil- 
tonian in which is the only nonzero interaction, for example, we encounter the Hamil- 
tonian diagrams in figure 3. A detailed derivation of the expression for the Hamiltonian 
diagram in figure 3a is given in Appendix B as an illustrative example, but throughout the 
text I simply give the appropriate expressions without explicit derivation. The diagram- 
matic rules for Hamiltonian diagrams are almost identical to the diagrammatic rules for 
time-ordered Green's functions given in Appendix A. The principal differences are that an 
infinite number of vertex rules are required, each of which is easily determined from the 
related interaction; and the rules for energy denominators differ, as discussed in Section 
IV. 

I have not bothered to indicate the energy denominators in figure 3, because as one 
sees in Eq. (4.14), there are only two choices in the second-order term. The denominator 
contains either a difference of the incoming state energy and the intermediate state energy 
or the outgoing state energy and the intermediate state energy. A transformation is com- 
pleted by evaluating these Hamiltonian diagrams with the appropriate cutoffs, and then 
rescaling the external variables. 

Given the transformations in Section III, and the techniques for constructing effective 
Hamiltonians in Section IV, we can readily compute the second-order behavior of the 
transformations 'near' the critical Gaussian fixed point in Eq. (5.24). The first step in this 
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process is to choose an initial Hamiltonian, H%°. If we choose a Hamiltonian for which 
only u 2 is nonzero, the linear analysis is exact, so we need to add at least one interaction 
to study second-order behavior. The simplest examples should result from adding only 
one interaction, and the strength of this interaction must be sufficiently weak that 
is near the Gaussian fixed point. Regardless of what interaction we add to the fixed 
point Hamiltonian to create H^°, we will find that contains an infinite number of 
interactions. Any interaction we add leads to an infinite number of irrelevant interactions 
at least, so I first consider the relevant or marginal interactions that can be added to the 
Gaussian fixed point. 

Let me specialize to the discussion of T 1 - and the boost-invariant transformations. 
For these transformations, there are relevant and marginal interactions with 113 nonzero, 
and there are marginal operators with nonzero. I assume that the symmetry — > — </> 
is maintained manifestly, so that only 114 is nonzero. If spontaneous symmetry breaking 
occurs, so that the 4> — > — <fi symmetry is hidden, u% must be allowed to appear in the 
Hamiltonian; however, in this situation it is a function of the variables appearing in the 
manifestly symmetric theory, as Wilson and I have discussed [39] . There can be no powers 
of transverse momenta in the marginal part of W4, because these lead only to irrelevant op- 
erators according to the analysis of Section V. On the other hand, functions of longitudinal 
momenta have no effect on the linear analysis, so we should consider 



H^; = H*+5H U (6.1) 



with, 



5H = ^ J dqi dq 2 dq s <% (16n 3 )5 3 (qi + q 2 + qs - 94) 

g(-Qi, -Q2 » -93" » it) at (?i) at (ffi) at (ffl)a(?4) 
+ ^Jdqi dq 2 dq 3 dq 4 (167r 3 )5 3 (gi + q 2 - qs - 94) ^ 

g(-Qi, -92", 93" , q! ) ^(91)^(92)0(93)0(94) 
+ \ J rf 9i dq 2 dq 3 dq^ (16n 3 )5 s (qi - q 2 - qs - q 4 ) 

K-Qi^Q^Qs^t) 0^91)0(92)0(93)0(94) , 

where g{qf) is the marginal part of U4. 

This is still too general for an initial discussion, because of the function g. The only 
constraints on this function are that it be symmetric under the interchange of any two 
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arguments with the same sign and that it be dimensionless, so that ratios of momenta 
must be employed unless there is a cutoff on longitudinal momenta that can be used to 
form g. I assume that such a cutoff does not exist, because one must typically employ T + 
to control a cutoff on longitudinal momenta, and I show below that perturbation theory 
cannot be used to construct T + . We will see below that spectator momenta also enter g 
for some transformations. Ultimately one would like to use Lorentz covariance to place 
restrictions on U4, but boost invariance alone places no further restrictions. Let me first 
concentrate on the simplest choice, 



9(Qi,Q2^3^t) = 9- (6-3) 

This leads to the familiar canonical </> 4 Hamiltonian, with no mass term. We will see that 
a mass term is always generated by the transformation. 

Figure 3 shows the second-order corrections to the effective Hamiltonian that arise 
when we use T 1 - or any other transformation. Let me begin by studying one of the 
corrections to 114, the first Hamiltonian diagram in figure 3b. The analytic expression 
corresponding to this diagram is determined by the second-order term in Eq. (4.14), and I 
further simplify the calculation by considering only the part of the expression that involves 
the energy of the incoming state. The relevant momenta are shown in the figure, but the 
analytic expression is drastically simplified if one chooses the Jacobi variables, 



(pl, p\) = (yP + , yP 1 + r x ) , ipt > pS) = ((1 - y)P + > (i - v)p ± - r x ) , 

(k+, k\) = (xP + , xP 1 + s x ) , (k+, k x 2 ) = ((1 - x)P+, (1 - x)P ± - s 1 ) . 



(6.4) 



Note that all longitudinal momentum fractions, such as x and y in Eq. (6.4), range from 
to 1. The total momentum is P. Using these variables one finds that the correction is 



Sv^ = - 
2 



«-9 



cPs^dx 



167r 3 x(l — x) Ly(l — y) x(l — x) 
9(Kl - (xP 1 + s ± ) 2 )0((xP ± + s x ) 2 - A 2 ) 
6(A 2 - ((1 - x)P ± - s ± ) 2 )6'(((l - x)P ± - s 1 ) 2 - A 2 ) 
9(Aj - (yP^ + r^)>(A 2 - ((1 - y)P^ - r^) 2 ) . 



(6.5) 



I use the notation 5v^ to indicate that the rescaling operations in T have not been 
completed. The first factor of 1/2 is seen in Eq. (4.14), while the second factor of 1/2 
is a combinatoric factor that can be seen in Eq. (A. 16) in Appendix A. The final step 
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functions are associated with the incoming and outgoing particles, whose momenta must 
satisfy < p 1 - < A 2 before the rescaling operations. Usually I do not display the step 
function cutoffs associated with the incoming and outgoing states, because they are always 
understood to be present. After completing the rescaling operations of Eqs. (3.10) and 
(3.11), and multiplying by z^, one obtains one of the second-order contributions to 8114. 

There are three important features of Eq. (6.5) upon which I want to focus. First, 
the cutoffs employed in T 1 - lead to a somewhat complicated analytical analysis, because of 
their dependence on the total transverse momentum of the incoming particles. However, 
one might normally be willing to live with this difficulty because it affects only irrelevant 
operators in 114. Such difficulties also occur for marginal operators in U2- The second feature 
of Eq. (6.5) is its dependence on and y. We started with a very simple expression for 
U4 given in Eq. (6.3), and after one application of T 1 - we generate a complicated U4. 
The renormalization group analysis is useful if we can show that most of this complexity is 
associated with irrelevant operators. Remembering that any powers of transverse momenta 
that appear in U4 are irrelevant operators, we should expand the integrand in powers of 
r 1 . This leads to 

g 2 fd 2 s ± dx 1 ^.vx(l - x) r ±2 i n 

^ 4 = _ T J ^ l^VyXTTy)^ 

n=0 6.6 
6>(Aq - (iP 1 + s- L ) 2 )6'((xP- L + s- 1 ) 2 - Af) 

6{k 2 - ((1 - x)P ± - s ± ) 2 )^(((l - x)P ± - s 1 ) 2 - A 2 ) . 
I do not display the cutoffs associated with the incoming and outgoing particles. There 

2 2 

may be little problem with the convergence of this expansion as far as the ratio r 1 - /s 1 - 
is concerned, because the cutoffs insure that this ratio is usually less than one. It is 
not guaranteed that this ratio is always less than one; however, we will find a far more 
serious problem. I choose Ai = Ao/2 in remaining calculations, unless specified otherwise. 
Consider the simplest case in which P 1 - = 0, so that 



c 9 2 J, (t> ,,^ ^ r(n + l) l~(l/2) 2 " r r^ 2 y\ 

6V4 = ~32^ \ H2) + ^ 2 2 "+i E(n + 3/2) 2n l y(l - y)A 2 \ ]' ^ 

After rescaling the momenta according to Eq. (3.10) we obtain a Hamiltonian that contains 



9 2 f, , n r(n+l) l-(l/2) 2n r r ±2 i - , 
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The first term in this expansion is identical in form to the interaction with which 
we began, and it causes the strength of this marginal operator to decrease as the cutoff 
decreases. If we add the correction involving the outgoing energy, the first Hamiltonian 
diagram in figure 3b alters the marginal part of 114 by a factor of — /n(2)/(167r 2 ) g 2 . The 
remaining terms all correspond to irrelevant operators. Before we can drop such irrelevant 
operators as a first approximation, however, we should show that not only do they occur 
with small coefficients, they also continue to lead to small corrections. In order to see 
that this is not the case, let me simplify the problem by considering a Hamiltonian that 
contains only the irrelevant interaction 



.I 2 „±- 2 



«4(<?i,<?2,<?3,<?4) =h (qf + <?2")(-|-^2 + ~Ta^) + permutations] . (6.9) 

While it may not be obvious, after a change of variables this interaction leads to the first 
irrelevant operator in Eq. (6.8). When we apply the transformation T 1 - to the Hamiltonian 
containing this irrelevant interaction we encounter Hamiltonian diagrams identical in form 
to those shown in figure 3; however, instead of the vertex in Eq. (6.3) we have the vertex in 
Eq. (6.9). Following steps analogous to those leading up to Eq. (6.6) one set of resultant 
terms is 



h 2 f d 2 s^dx 1 v^rxfl — x)r^ i n 



W4 = ~T J 16^(1 - x) tf^yO^yjs^ 

6(A 2 - (xP 1 - + s ± ) 2 )^((xP ± + s^) 2 - A 2 ) 

6(A 2 - ((1 - x)P± - s^)>(((l - x)P^ - s^) 2 - A 2 ) 



(6.10) 



While this is almost identical to Eq. (6.6), the n = term is infinite. This problem 
inevitably results when one applies T^; but it is much worse than this of course, because 
one generates arbitrarily high inverse powers of longitudinal momentum fractions, as is 
apparent in Eq. (6.7). The cutoff's do not prevent these fractions from becoming arbitrarily 
small, and divergences inevitably result. Let me note that T 1 is the type of transformation 
one must employ when transverse lattice regularization is used and one wants to vary the 
transverse lattice spacing [102,103,111]. 

Placing cutoffs on the longitudinal momentum, as required if we use T + , alters this 
problem slightly; but it does not cure the problem. Once we have placed a cutoff on 
longitudinal momentum, eo, we are free to consider the effective four-point interaction 
for particles that have much larger longitudinal momentum, P + , than this cutoff. The 
longitudinal momentum cutoff then shows up in Eq. (6.10) as a small cutoff on the x 
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integration, leading to a factor of ln(P + /eo) instead of oo. This is an improvement, but the 
analysis of higher-order irrelevant operators in Eq. (6.17), for example, leads to arbitrarily 
large powers of _P + /eo, and we have no way to prevent this ratio from producing arbitrarily 
large coefficients. The sceptical reader should explore this problem much further, but I do 
not believe that there are any simple solutions short of abandoning perturbation theory. 

All of the above problems can be traced to the fact that some of the energies of states 
retained by T- 1 , as determined by the fixed point Hamiltonian, are larger than some of 
the energies of states removed. The energy denominators must be expanded in powers of 
the energy of incoming and outgoing states to separate relevant, marginal, and irrelevant 
operators; and this expansion does not converge. The problem is extremely severe for T 1 - 
because states of infinite energy are retained, leading to infinite errors when irrelevant 
operators are dropped. Similar problems are encountered when one studies T + . The 
cutoffs in T + require one to integrate over transverse momenta from zero to the transverse 
momentum cutoff. For massless Hamiltonians this leads to transverse infrared divergences, 
and even for massive Hamiltonians one again finds that the irrelevant longitudinal operators 
end up producing arbitrarily large effects if the external transverse momenta are large. I do 
not detail this problem further, because the calculations are quite similar to those above, 
with the only changes being in the cutoffs that occur inside the integrals. 

If we consider a boost-invariant transformation that places any cutoff on transverse 
momenta without limiting small longitudinal momenta, we encounter exactly the same 
problems discussed above. Thus, there is no reason to consider the transformation that 
utilizes the cutoff in Eq. (3.22). The cutoffs in Eqs. (3.23) and (3.24) limit the longi- 
tudinal momenta, and I show that the transformations associated with them apparently 
avoid the above problems. However, both of these cutoffs involve a sum over all particles 
present in a given state. As we shall see, this leads to second-order corrections that are 
spectator-dependent. Spectator dependence may not drastically complicate a perturbative 
renormalization group analysis. 

The boost-invariant transformation that utilizes the cutoff in Eq. (3.23) avoids the 
problems found for T in second-order. Using this transformation, and starting with a 
Hamiltonian that contains the interaction in Eq. (6.3), the correction in figure 3b becomes, 



^ 4 = T 



d?s dx 



167T 3 ;r(l — x) Lj/(1 — y) x(l — x) 



-l 



eikl 



Mifr^- A ?M A ?-j(fe))- 



(6.11) 



x(l — x) J \x(l — x) 



This is identical to Eq. (6.5) except for the cutoffs. There is an important assumption that 
is made in writing Eq. (6.11). I have assumed that there are no spectators, so that the 
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interactions displayed in figure 3 occur without any disconnected lines present. I return to 
this issue later, and discuss spectator effects. At this point we can expand the denominator 
exactly as was done above, but now the cutoffs guarantee that 



> 



x(l-x) y(l-y) 
Following the same steps that led to Eq. (6.7) we find 
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5V4 = — 



4 



d 2 s- L dx 1 



1 v^ pU- 



x r 



l 6yr 3 s ±2^Ly(l-y) s ±2j 



9 



32^ 



ln(2) + 



(1 — x)J \x( 
1 - (l/2) 2n 



< A ° " -TT-^rW-TT-^ " A ') 
V x 1 — x / \x(l — x) / 



n=l 



2n 



y(i-y)A? 



(6.13) 



I have not displayed the step function cutoffs on the incoming energy. To complete the 
transformation we need to rescale according to Eq. (3.10), and as we found above in Eq. 
(6.8), the only effect this has on Su^ is to change Ai in Eq. (6.13) into Ao 

This result is similar to the result for T^, and to see that problems do not arise when 
we add these irrelevant operators and apply the transformation a second time we need 
to follow the steps leading to Eq. (6.10). Starting with a Hamiltonian that contains the 
interaction given in Eq. (6.9), and applying an invariant- mass transformation, we derive 



8v4 = — — 



h 2 



(fis^dx 



4 J 167r 3 a;(l-x) A^ ^Ly(l - y) s ± 



^ 00 
Al^J 



x(l — x) r 



1 2 



9[Al 



3 ln(2) 
+ 



2 I 2 

6 



x(l — x) J \x(l — x) 



A?) 



32tt 2 I 8 4 y(l-y)Af 



1 - (1/2) 



2n-2 



n=2 



2n-2 Vy{l-y)A{ 



(6.14) 



Here again the transformation must be completed by rescaling the momenta, which again 
replaces Ai with Ao- Not only is every operator in the sum bounded by cutoffs on external 
momenta, every coefficient in the expansion is small when h is small. Since these irrelevant 
operators occur at 0(g 2 ) : the corrections in Eq. (6.14) are 0{g 4: ). This should give us 
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hope that irrelevant operators do indeed lead to small corrections in an analysis that starts 
by simply discarding such operators at each stage, at least while the Hamiltonian remains 
near the Gaussian fixed point. 

Unfortunately, it is not sufficient to show that the coefficients of the irrelevant operators 
are small, and that each subsequent correction produced by each irrelevant operator is 
small. One should also worry about the convergence of the corrections produced by the 
entire sum of irrelevant operators. For example, the reader should be concerned with the 
convergence of the sum in Eq. (6.13). Completing the transformation by rescaling i- 1 , and 
completing the sum leads to 



Su 4 = -#J/n(2) + l - ln( 4y{1 - y)A °- r±2 2 ) Wa 2 , - f \ . (6.15) 
32tt2| {) 2 V 4y(1 _ y)A 2_ 4r ±2;j V° y{l-y)) 

I have restored the cutoff on the incoming energy, and one can clearly see a potential 
problem coming from the fact that the logarithm diverges when the incoming energy 
reaches its maximum value. This divergence occurs because the energy denominator in 
Eq. (6.11) vanishes on a surface when the incoming energy equals k\. However, this 
divergence is simply being buried in the irrelevant operators. How can this make sense? 

There are two issues one must address. First, one must determine whether this di- 
vergence actually shows up in observables if the Hamiltonian containing the interaction 
in Eq. (6.15) is solved. While I do not go through the details, if one studies two-particle 
scattering in perturbation theory with such a Hamiltonian, the logarithmic divergence in 
Eq. (6.15) enters first-order perturbation theory with a strength of 0(g 2 ). In second- 
order perturbation theory there is an additional contribution of 0(g 2 ) coming from the 
interaction in Eq. (6.3). This contribution also has a logarithmic divergence that tends 
to cancel the divergence from the interaction in Eq. (6.15). If one drops the irrelevant 
operators, one finds that this latter divergence is not canceled and a large error is made if 
the invariant-mass of the scattering state is near the cutoff. Of course, one expects large 
errors near the cutoff; but well below the cutoff one hopes that errors are small even when 
irrelevant operators are dropped, and that the results are systematically improved as the 
leading irrelevant operators are retained. At least to leading orders in the coupling it is 
easy to verify that this happens. 

The second issue is directly relevant to the renormalization group. One must determine 
whether the divergence in Eq. (6.15) causes the irrelevant operators to have a significant 
effect on the next Hamiltonian produced by a transformation. To determine this, include 
the entire logarithm in a vertex and study the second-order correction in the first Hamil- 
tonian diagram of figure 3b, using this logarithmic vertex in combination with the original 
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vertex from Eq. (6.3). Keeping only the marginal part of Sv4, and not worrying about 
combinatoric factors, we obtain 




/ 



x(l — x) I 



Zn(2) + £ 



1 ~ (1/2) 
2n 




167r 3 x(l — x) 
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x(l - x)Aq 



(6.16) 



x(l — x) 



I have restored the sum to facilitate the discussion of convergence. The integrals are easily 
completed, leading to 



Clearly this sum converges. It is also clear that it converges rather slowly, with an 
error that falls off as the inverse power of the number of terms included. If one wants 
to include corrections of O(g^), such as the third-order corrections to the transformation 
studied in the next Section, without making a large {e.g., 25%) error in the coefficient 
of these corrections, it is also necessary to include some irrelevant operators. However, if 
g is small, it is apparent that these irrelevant operators produce small corrections to the 
second-order analysis, despite the fact that the logarithm arising in Eq. (6.13) diverges. 
This is possible because the divergence is integrable. 

Next I want to show that the transformation generates a mass term in U2, and that the 
mass counterterm required when we try to keep the physical mass zero is rather unusual. 
In figure 3a I show the second-order correction to the Hamiltonian that affects U2, and in 
Appendix B I analyze this correction for arbitrary 114. I use the interaction in Eq. (6.3). 
To proceed, define the Jacobi variables 




(6.17) 




0.480 + 0.141 + 0.055 + 0.027 + • • • 




k x i) = (xP+, xP ± + q X ) , k X 2 ) = (yP+, yP 1 + r x ) , 

(k+,k L 3 ) = (zP + ,zP L + s L ). 



(6.18) 



Using these variables the second-order correction to u<i is 
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3! J 



j2„_Lj„, j2„±j„ r „_L 2 „_L 2 „_L 2 



r s- 
— + + - 

x y z J 



d q dx d r dy d s dz 
16tv 3 x 167r 3 y lQn^z 
(16tt 3 )5(1 -x-y- z)5 2 (q L + r 1 + s 1 ) (6.19) 

I 2 I 2 I 2 I 2 I 2 I 2 

»(a»-S__£ !_WS_ + £_ + ! a?). 

V x y z / V x y z J 

This expression does not depend on P, so it leads only to a mass shift and does not change 
the marginal operators in u 2 or produce any irrelevant operators in u 2 . This conclusion 
changes when spectators are included. 

It should be obvious that the integral does not vanish, and therefore u 2 develops a 
negative mass-squared term if we start with a Hamiltonian that has no mass term. 

It is important to note that the mass shift is not infinite, despite the fact that particles 
with small transverse momenta are removed. If transverse infrared divergences appear, 
transverse nonlocalities may follow. Consider a Euclidean field theory and place a cutoff on 
the four- momentum squared, q 2 . When states with a given large range of q 2 are removed, 
it is not required that every component of the momentum be small. If q$ is large, for 
example, states with small values of q% are removed, and one might naively worry that this 
produces long-range forces in the x 3 direction. Of course this does not happen, because 
of rotational invariance. In light-front field theory the invariant-mass cutoff violates strict 
rotational invariance; however, it retains some features of this symmetry and apparently 
allows one to remove states with small transverse momenta without producing long-range 
transverse interactions, at least in low orders of perturbation theory. The most interesting 
examples of this principle are discussed later in this Section. 

As discussed in Section II, we want to choose the Hamiltonian so that gives 
reasonable results when it is diagonalized. Since there is no inverse transformation, this 
process is typically trial and error in a nonperturbative analysis. One first tries a particular 
H^°, and constructs H^.. If does not produce reasonable observables, must be 
altered. While there are an infinite number of operators that can be adjusted in P^°? only 
the relevant and marginal operators are expected to produce effects that survive after many 
applications of T; therefore one hopes that it is necessary to control only a finite number of 
readily identified terms in to produce desired results in . One also hopes that the 
irrelevant operators in H A ° are less important than the relevant and marginal operators, 
when is diagonalized, but there is no guarantee of this and it is not essential to 

the renormalization group analysis. If this happens, one is led to consider the relationship 
between the coefficients of the relevant and marginal operators in to the corresponding 

coefficients in H^.. There is no guarantee that this relationship is simple. It could even 
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be chaotic, in which case one may want to find a new problem. However, when one 
evolves Hamiltonians near a Gaussian fixed point, the relationship between coefficients 
at the beginning and end of a long trajectory should not be highly nonlinear, as each 
transformation introduces only small nonlinearities. Wilson has given a general discussion 
of this issue [13], and the main new ingredient in the light-front renormalization group is 
the appearance of arbitrary functions of longitudinal momenta. 

The mass terms in ui provide a simple example. I refer to any terms in 112 that do not 
depend on transverse momenta, regardless of their dependence on longitudinal momenta, 
as mass terms. For the purposes of this discussion I simply regard the mass term that 
appears with no dependence on longitudinal momentum in 112 in as the physical mass. 

In reality one must still solve the Schrodinger equation with H^°. to relate the mass in 

to an experimentally observed mass, and it is even possible that there is no experimentally 

observed mass that directly corresponds to the mass in H^°. 

While the mass terms in H^° N are affected by all of the terms in H^°, if one uses a 
second-order approximation of T to construct a Hamiltonian trajectory, any mass term 
in appears directly in H^.. Thus, in a second-order analysis, it is trivial to control 
the mass term in H^.. If we want to produce a Hamiltonian in which the physical 
mass is zero for example, we add a mass 'counterterm' to and adjust it to cancel the 
entire trajectory of mass shifts that begins with the shift in Eq. (6.19). This simplicity is 
lost in third- and higher-order analyses, as we will see in Section VII. On the other hand, 
we will find that the coupling constant coherence conditions fix the mass counterterm to 
each order in perturbation theory to a value that yields a massless physical particle in that 
order of perturbation theory [39]. 

Let me return to the issue of spectator-dependence. Figure 4 is identical to figure 3a, 
but I have added a double line to represent all spectators. The calculation proceeds almost 
exactly as before, except we must include the effects of the spectator momentum. I use the 
same variables given in Eq. (6.18), and I assume that the incoming particle momentum 
and the spectator momentum are respectively, 



(wP+.^ + t 1 ), ((1 - w)P + , (1 - w)P- L - t x ) . (6.20) 

We also need to know the invariant mass-squared of the spectator state, which I assume 
to be M 2 . Even if the particles are massless, we cannot assume that the invariant mass of 
the spectators is zero, of course. The spectator energy cancels in all energy denominators 
when no interactions occur between spectators. Such interactions produce disconnected 
marginal and irrelevant operators that I do not discuss, because they are not important 
in low orders of a perturbative analysis. The second-order correction to ui from the 
Hamiltonian diagrams in figure 4 is 
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^ g 2 fd 2 q- L dx d 2 r ± dy d 2 s ± dz "t ±z q ±2 r ±2 s^i- 1 

3! J 167T 3 X 167T 3 y 167T 3 Z I. w x y z J 

(16tt 3 )5(«; - x - y - ^(t 1 - q x - r 1 - s x ) 

t ±2 + M2 q± 2 ^2 g± 2 (6 . 21) 
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One can compare this result with Eq. (6.19) to see the effect spectators have on 5v2- 
In Eq. (6.19) 8v2 does not depend on the longitudinal momentum of the incoming boson, 
or the transverse momentum of the incoming boson. In Eq. (6.21) 5v2 depends on the 
longitudinal momentum fraction of the incoming boson, w; and it depends on the relative 
transverse momentum of this boson with respect to the rest of the system. 

When faced with any correction to the Hamiltonian, we are supposed to expand the 
correction in terms of relevant, marginal and irrelevant variables. In this case, this means 
we are supposed to expand Sv2 in powers of the transverse momentum t- 1 . In addition, 
we are supposed to expand in powers of M 2 . In the fixed point Hamiltonian the particles 
are massless, and M 2 is a function of the transverse momenta of the spectator particles 
that goes to zero when these momenta go to zero. Remember that masses are treated as 
perturbations. The mass terms in Sv2 are found by setting all transverse momenta to zero, 
leading to the relevant part of Sv2, 
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This expression is identical to the similar limit for Eq. (6.19), except for one very important 
difference. The longitudinal momentum fractions in the integrand add to w instead of 1. 

We need to understand how the mass terms that arise in U2 after repeated application 
of T depend on longitudinal momenta, and here this problem is equivalent to understanding 
how the mass depends on w. This dependence is easily worked out by changing variables, 
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Using these variables, and using the fact that Ai = Aq/2, we obtain 
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The remaining integral is a finite number that has no dependence on any momenta, so 
we have exactly determined the dependence of the new mass term on the longitudinal 
momentum fraction. The new mass term is proportional to the longitudinal momentum 
fraction, unlike the physical mass which is independent of this fraction! 

The next step in the analysis is to complete the calculation of marginal and irrelevant 
operators that arise in Eq. (6.21). The mass term is the only relevant operator occurring 
in Eq. (6.21), and the only marginal operator is the piece of Sv2 that is quadratic in the 
external transverse momenta, including M 2 . Using the new variables in Eq. (6.23), we 
find 
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The factor of t 1 - /(wAq) in the energy denominator and the factor of t _L /( v /«;Ao) in the 
momentum conserving delta function each produce a quadratic term in 5v2 proportional to 

2 

t , with no dependence on longitudinal momentum. This produces a term in the energy 
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proportional to t 1 / (wP + ). The terms proportional to t 1 - + M 2 in the cutoffs produce a 
term in 5v2 that is proportional to w^ 1 - + M 2 )/(l — w), which then produces a term in 
the energy proportional to (t^ +M 2 )/[(1 — w)P + ]. The factor M 2 is itself a sum of terms 
that are each of the form fx, with x being a longitudinal momentum fraction and 
being a relative transverse momentum, if there are no physical masses in the theory. These 
corrections are similar to the fixed point u\, but with an important difference; they do not 
depend on the total transverse momentum. 

Assuming U2(q) = q -1- in the Hamiltonian in Eq. (3.7) and using Jacobi variables 
(q i,q^) = {xiP + ,XiP- L + r i), we can write the terms involving 112 using projection 
operators and obtain 
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The corrections to U2 that are of Oft 1 - , M 2 ) from Eq. (6.25) alter the coefficient of each 

factor r^j/(xjP + ) in Eq. (6.26). These corrections differ from standard wave function 
renormalization, because the coefficient of each factor P- l2 /P + in Eq. (6.26) remains 
1, and wave function renormalization would alter this coefficient also. These corrections 
are marginal operators, and one can include them in a fixed point Hamiltonian if U2 is 
allowed to be spectator-dependent, so that it can depend not only on the momentum of 
a single particle, but also on the total momentum of the state. Eq. (6.24) proves that U2 
must be spectator-dependent if an invariant-mass transformation is used. In the original 
discussion of Hamiltonians in Section III, it was simply assumed that U2 depends only on 
one momentum, and now we are finding an example in which the transformation forces 
us to expand the original definition of the space of Hamiltonians. This point is clarified 
further below. 

In order to control these spectator-dependent corrections to U2 in the final Hamiltonian, 
H^., we must allow such terms to appear in H^°. Since these counterterms are part of 
5H[, they do not modify the second-order behavior of the transformation, but they do enter 
at third order. If the scalar particles appear as asymptotic particles (i.e., are not confined), 
we need to precisely cancel these corrections to U2 to obtain the appropriate dispersion 
relation for the physical scalar particles, as we are forced to introduce mass counterterms 
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to cancel the corrections found in Eq. (6.24). Just as it is trivial to cancel any mass that 
arises in a second-order analysis, it is trivial to cancel these marginal corrections to it 2 . 
Again, these corrections may have nontrivial effects in a third-order analysis. 

The analysis of 114 has not been completed, and we must evaluate the remaining Hamil- 
tonian diagrams in figure 3b. It is convenient to use the variables 



(p+ p\) = ( Xi P + , Xi P L + r\) , (k+, k\) = ( Vi P+, y.P 1 - + s\) . (6.27) 

I am only interested in the correction to the marginal part of 114, so these diagrams are 
evaluated with the external transverse momenta and the total transverse momentum set to 
zero. It is only necessary to explicitly evaluate the second and sixth Hamiltonian diagrams 
in figure 3b, as others are simply related to these two. 

The second Hamiltonian diagram in figure 3b, combining both second-order terms in 
Eq. (4.14) involving incoming and outgoing energy, leads to a correction of the marginal 
part of 5v4, 
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To simplify this expression and remove dependence on external momenta from the inte- 
grand, we can change variables, 



Vi = (x2 ~ X4)Zi , s ± i = ^x 2 - x 4 A q- L i . 
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This leads to 
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This correction is of exactly the same form as the correction to the marginal part of 
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U4 from the first Hamiltonian diagram in figure 3b. The correction from the third diagram 
in figure 3b is identical to the correction in Eq. (6.30), the only difference being that the 
third diagram survives when X4> X2- The fourth and fifth diagrams contribute the same 
amount to the marginal part of 8V4 as the second and third diagrams. The total correction 
to the 2-particle — > 2-particle part of U4 from the first five diagrams in figure 3b leads to 
the second-order transformation, 



3/n(2) 2 

9^9--^9 2 . (6.31) 

It is a straightforward exercise to compute changes in the irrelevant parts of U4 by 
allowing the external transverse momenta to be nonzero and expanding in powers of these 
momenta. It is also straightforward to determine the effects of spectators on the change 
of 114. I return to the discussion of irrelevant parts of U4 below. Note that spectators 
apparently have no effect on the marginal part of U4 for the simple vertex in Eq. (6.3). 
The transverse momenta of the spectators are set to zero when computing the marginal 
part of 114, so the cutoffs in Eq. (6.30) are not affected. The only effect is in the momentum 
conserving delta functions, and using variables similar to those in Eq. (6.23) one explicitly 
recovers Eq. (6.30) unchanged. 

We will see later that the marginal part of 114 is inevitably spectator-dependent. Higher- 
order corrections produce spectator-dependence, and Lorentz covariance and cluster de- 
composition require spectator-dependence. We must adjust the Hamiltonian so that all 
observers in frames related to one another by rotations obtain covariant results, and so 
that systems of particles that are not causally connected do not affect one another. 

The last six Hamiltonian diagrams in figure 3b lead to a correction of the marginal 
1-particle — > 3-particle and 3-particle — > 1-particle parts of U4 identical to Eq. (6.31). 

The first diagram in figure 3c does not occur, because the intermediate state always 
has less energy than the external states, which is not allowed in a second-order correction. 
The second and third Hamiltonian diagrams are allowed and lead to the irrelevant operator 
uq. There are additional contributions to uq that I do not show. The cutoffs require 
that the external states have an energy below the cutoff while the intermediate state has 
an energy above the cutoff. Moreover, momentum conservation completely determines 
the momentum of the single internal particle line. Using the coordinates (pf,P i) = 
(xjP + , X{P L + r i) and ignoring spectator effects, one part of the contribution to 5uq from 
the second diagram in figure 3c is 
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(6.32) 



Here I have included the cutoffs associated with the external lines, because they are im- 
portant for further analysis of this term. 

At this point we should try to proceed by expanding this operator in terms of increas- 
ingly irrelevant operators. All terms in uq are irrelevant if inverse powers of transverse 
momenta do not arise, but this correction seems problematic. We should obtain the lead- 
ing correction by letting the external transverse momenta go to zero. However in this limit 
the energy denominator vanishes and the cutoffs go to zero, forcing us to analyze oo x 0. 
I assume that all external longitudinal momenta remain finite as the external transverse 
momenta approach zero. The leading correction to uq in this limit becomes 

^^(Ag-^^-A?), (6.33) 

where — > 0, and y = k + /P + . Here k- 1 is the transverse momentum carried by the 
internal boson line, and y is its longitudinal momentum fraction. 

To analyze this distribution, let us integrate it with a smooth function of y, f{y). This 
leads to the integral 
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Of course I have used Ai = Ao/2 again, and in the last line I have finally completed the 
limit in which all external transverse momenta are taken to zero. Thus we see that the 
distribution that corresponds to the leading correction to uq is a delta function in the 
longitudinal momentum transfer, and as expected the coefficient of the leading correction 
to uq is inversely proportional to Ajj, as an irrelevant operator should be. No long range 
transverse interactions are produced by the elimination of high energy states, and inverse 
powers of transverse momenta do not appear when one expands the second-order transfor- 
mation in terms of relevant, marginal and irrelevant operators. The appearance of delta 
functions of longitudinal momentum transfer may have interesting consequences, but they 
are not important in low orders of a perturbative analysis. 

In figure 3d two of the disconnected Hamiltonian diagrams that affect us are displayed. 
While such disconnected diagrams cancel when there are no cutoffs on the intermediate 
energies, this cancellation does not occur when the cutoffs are in place. It is possible for the 
incoming and outgoing states to have energies just below the cutoff, while the intermediate 
states have energies just above the cutoff. I do not go through a detailed analysis, because 
the correction is again irrelevant and local in the transverse direction, in this case being 
proportional to 1/Ag. 

This completes the initial analysis of the second-order transformation for a massless 
Hamiltonian with the simple interaction given in Eq. (6.3). There are two remaining topics 
in this Section. First, I want to study repeated second-order transformations. Second, I 
want to discuss Lorentz covariance and cluster decomposition in second-order perturbation 
theory, and introduce coupling coherence [39,94-100]. 

Let us assume that the second-order behavior of the invariant-mass transformation 
is a reasonable approximation of the complete transformation. To compute an example 
trajectory of Hamiltonians, H^°, we can choose to be, 
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Note that in the critical Gaussian fixed point Hamiltonian, £, and 5 are all zero. This 
Hamiltonian contains a complete set of relevant and marginal operators required for an 
approximate second-order analysis. In this example there is no sector-dependence in the 
marginal part of 114, so we are able to write the marginal interaction without using the 
awkward projection operators that are required to display u 2 - When the second-order 
invariant-mass transformation is applied to this Hamiltonian, and irrelevant operators are 
dropped at every stage, the resultant Hamiltonian contains no new relevant or marginal 
operators. In this case the only effect on the relevant and marginal operators is to change 
the constants £0? ^0? an d go- A complete second-order analysis would require us to consider 
more general interactions than shown in Eq. (6.3), and to retain irrelevant operators. 
Irrelevant operators typically produce new relevant and marginal operators. Dropping 
irrelevant operators, we can write the relevant and marginal operators in as in Eq. 
(6.35), using the constants £ n , fj, n , and g n . 

The approximate second-order transformation can be summarized by the equations, 



g n +i = g n -c g g n , 



(6.36) 



£n+l =£,n + C£ g. 



(6.37) 
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(6.38) 
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Note that if g n is small, the initial assumption that the second-order corrections are small 
in comparison to the linear corrections is consistent. I have shown that the constants c g 
and are positive, but have not analyzed c^. It is fairly easy to see from Eq. (6.25) that 
is also positive. If g n is sufficiently small, these equations should be a reasonable crude 
approximation to the entire transformation. 

Wilson has discussed how one must solve such equations so that errors are controlled 
[13], and I have already repeated some of his discussion in Section II. Eqs. (6.36)-(6.38) are 
much simpler than the general case discussed by Wilson, and they are readily solved for 
large J\f. For a massless scalar theory, [Mq and £o can be adjusted so that the physical particle 
has the dispersion relation of a massless particle. If one insists on actually computing the 
value of /Uq required to obtain a specific value of /rj^-, one finds that fiQ must be controlled to 
an accuracy of about 1 part in . In practice, one is never interested in ,Uq directly, and 
the rest of the calculation should be adjusted so that no equations depend on the precise 
value of Hq. It is easiest to fix and solve Eq. (6.38) towards decreasing n. £yv should 
be close to 0, and £o is adjusted to achieve this result. Through second-order in 5H 2 , the 
equation for g n decouples from the equations for £ n and fx n ; and the first step in solving the 
complete set is to solve Eq. (6.36). As successive transformations are applied, Eq. (6.36) 
shows that g decreases. We must start with a sufficiently small value of g for perturbation 
theory to be reasonable. An accurate solution of Eq. (6.36) can easily be constructed by 
iteration. A reasonable approximation over any finite segment of the trajectory is 



9m 9m 9m ,„ oa n 

9n = --, 7 \ = 7 x = 7 x , (6.39) 

l + c g (n-m)g m 1 + J ^ln(^)g m l + I fcZn(^ ra 

where in the last step I have used the result for c g in Eq. (6.31). The error in this 
approximation grows as \n — m\ becomes large. It is interesting to note that the factor of 
ln(2) in c g , which came from the choice Ai = Ao/2, drops out of this final result. The result 
is well-known. In the renormalized Hamiltonian, which as discussed above is obtained by 
allowing M — > oo, g = 0. This analysis is not complete, of course, because I have only 
considered the case where go is small; however, I am only interested in showing that a 
perturbative analysis may be possible. 

This brings me to the final subject for this Section. Up to this point, there has 
been no discussion of how Lorentz covariance and cluster decomposition, both of which 
are violated by an invariant-mass cutoff, are restored in physical predictions. Perhaps 
the most important observation is that if one succeeds in restoring these properties in 
predictions with one cutoff, they are restored for all cutoffs, because the renormalization 
group is designed to preserve the matrix elements of observables as cutoffs are changed. 
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I give three simple examples in second-order perturbation theory. The first example 
is the dispersion relation for a single boson in the presence of spectators. I compute the 
second-order contribution to the Green's function with the external propagators removed; 
which is also the second-order shift in the invariant-mass-squared of the state when e (see 
Appendix A) is chosen to be the state's on-shell free energy. If this invariant-mass-squared 
shift does not transform as a scalar under Lorentz transformations, and/or it depends 
on the spectators, we must add counterterms to restore Lorentz covariance and cluster 
decomposition. The second example is the boson-boson scattering amplitude; for which 
the second-order correction comes from diagrams identical to the first five in figure 3b, 
with spectators added. This amplitude should be manifestly covariant, depending only on 
the invariant-mass of the two bosons that scatter from one another, with no dependence on 
spectators. If these conditions are not satisfied, we must again add counterterms. Finally, 
I list the correction to the one-boson to three-boson Green's function, given by the fifth, 
sixth, and seventh diagrams in figure 3b. 

The main differences between second-order perturbation theory diagrams and second- 
order Hamiltonian diagrams are the ranges of integration and the energy denominators. In 
Hamiltonian diagrams the range of integration is bounded by an upper and lower cutoff, 
and the energy denominator always involves the on-mass-shell energies of incoming, out- 
going, or intermediate states. In time-ordered perturbation theory only the upper cutoff 
appears, and the energy in the denominator is arbitrary. We have already seen that a 
'mass' term appears in 112 with the wrong dispersion relation when there are spectators, 
Eq. (6.24); and we will find a mass with the wrong dispersion relation when we study the 
invariant-mass of a boson in time-ordered perturbation theory. From the point of view of 
perturbation theory, we need to add a counterterm that exactly cancels this mass term or 
we do not obtain an invariant-mass shift that transforms like a scalar. The counterterm 
must completely cancel the mass shift, and I show that one obtains the same counterterm 
from the renormalization group Eqs. (6.36) and (6.38) using the condition that [i n is a 
function of g n with no further dependence on n. This is a coupling constant coherence 
condition. 

The second-order contribution to the boson invariant-mass corresponding to figure 4 
is readily determined using the diagrammatic rules in Appendix A. I assume that the 
total longitudinal momentum is P + and that the total transverse momentum is P , with 
e = P~ = ^p+- for massless bosons. I assume that there are spectators, with all relevant 
momenta given in Eq. (6.20). The invariant-mass of a state with a boson in the presence 
of a massless spectator consists of a kinetic energy term and a mass term for the boson 
itself. At this point I am only interested in the violation of Lorentz covariance coming 
from the mass term that appears in the shift, so I set the relative transverse momentum 

= and assume that the invariant-mass of the spectators is zero. In other words, I 
set all relative transverse momenta equal to zero. In this case, using the Jacobi variables 
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defined in Eq. (6.18), the boson mass shift is 
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Lorentz covariance and cluster decomposition require that this shift be a constant, inde- 
pendent of all longitudinal momenta. 

Following the same steps leading to Eq. (6.24), we find that this mass shift can be 
written as 
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This can be related to the constant c M that appears in the renormalization group Eq. 
(6.38), and one finds that the mass shift is 
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where we can use Eq. (6.24), to show that 
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(6.43) 



To obtain this result for one must remember that 8114 contains an extra factor of four 
not found in Eq. (6.24) that results from the rescaling part of the transformation. To go 
from Eq. (6.41) to Eq. (6.42), note that the integral in Eq. (6.41) can be written as an 
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infinite sum of integrals with upper and lower cutoffs, with each successive cutoff being 
1/4 the previous cutoff. In each of these integrals one can rescale momenta in exactly the 
manner used to obtain Eq. (6.24), leading to a sum 1 + 1/4+1/16 + • • • = 4/3. 

The factor of w that appears in Eq. (6.42) shows that the mass shift is neither covariant 
nor spectator-independent, and a mass counterterm must be added to exactly cancel this 
shift. 

Return to Eqs. (6.36) and (6.38) and ask whether it is possible for ji\ to be a pertur- 
bative function of g n with no further dependence on n. In general we can choose any initial 

value for hq and solve Eq. (6.38) to find jji n ; but if we assume that //^ = (ag n + f3g^-\ )Aq, 

and substitute this into Eq. (6.38), using Eq. (6.36), we find a = 0, j3 = c^/3. To this 
order we have a unique result, ^ = (c^/3) g\ Aq. Higher order terms are not determined, 
because Eqs. (6.36) and (6.38) are altered at 0(g 3 ). For this choice of //^, is uniquely 
determined, 



= C Jf A% + O(gl) . (6.44) 

The factor of w found in Eq. (6.42) is included in the definition of the term in which \i n 
appears, as seen in Eq. (6.35); therefore, this value of /j/q precisely cancels the entire mass 
shift found in second-order perturbation theory and acts to restore covariance and cluster 
decomposition. In other words, with no direct reference to these properties, one can use 
the renormalization group and the coupling constant coherence conditions to remove the 
violations caused by the invariant-mass cutoff. 

Let me next consider the boson-boson scattering amplitude corresponding to the first 
diagram in figure 3b, with spectators added. I assume that the total longitudinal and 
transverse momenta of the two bosons that scatter are given by the first momentum in 
Eq. (6.20), with the spectator momentum being the second momentum. For simplicity I 
assume that the invariant-mass of the spectators is zero. For the bosons that scatter I use 
the Jacobi variables, 



(k+, k\) = (xwP+, xiwP 1 + t 1 ) + s x ), 

(k+, k x 2 ) = ((1 - x)wP+, (1 - x)(wP ± + t 1 ) - S" 1 ) . 



Letting e = (P + M )/P + , the scattering amplitude is 
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The real part of this amplitude is 
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Remember that M 2 is the invariant-mass of the entire state, including the spectators. 
The invariant-mass of the two-boson subsystem is 
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The boson-boson scattering amplitude is neither covariant nor spectator-independent. For 
massless bosons with 



(Pt, p\) = (ywP + , y(wP ± + t x ) + r x ), 

(p+ p x 2 ) = ((1 - (1 - y)(«;P ± + t x ) - r x ) , 

it is also easily seen that 
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The on-shell scattering amplitude blows up in perturbation theory as the relative transverse 
momentum between the scattering bosons goes to zero, but the important observation is 
that counterterms must be added to restore covariance and cluster decomposition. To 
discuss the counterterms, I set t 1 - = for simplicity. To determine what counterterms are 
required, expand the amplitude in Eq. (6.47), obtaining 
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The first term is covariant and spectator-independent, and after the subtraction of a 
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constant associated with coupling renormalization it yields the correct result. All remaining 
terms must be canceled by counterterms. The counterterms are part of U4, and it is only 
the second term in the series that affects the marginal part of u±. In this case we find 
that covariance and cluster decomposition requires us to depart from the simple marginal 
interaction in Eq. (6.3), and use 

Not only are functions of longitudinal momentum allowed in the marginal part of 114, 
they are required by covariance and cluster decomposition when we use an invariant-mass 
cutoff. Calculations of the remaining diagrams in figure 3b reveal additional logarithmic 
corrections to the vertex, but nothing qualitatively new is revealed. Letting p^ = zwP + 
and p\ = (1 — z)wP + , the real part of the complete scattering amplitude can be written 
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The additional counterterms required to restore covariance and cluster decomposition to 
the entire amplitude are easily determined. 

We can think of the corrections to 114 as new marginal and irrelevant variables, and 
as such we expect to find renormalization group equations that show how their strengths 
change with the cutoff. On the other hand, from the point of view of perturbation theory 
the strengths of these counterterms are determined by g; and they change with the cutoff 
only because g changes with the cutoff. 

It is straightforward to compute the one-boson to three-boson Green's functions cor- 
responding to the diagrams in figure 3b, and the real part is, 
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Here I have used = wP + , p\ = wxP + , p% = wyP + , and p\ = wzP + . Again, we find 
that covariance and cluster decomposition are violated, and that counterterms must be 
added to 114 to subtract these violations. 

Next I want to use the coupling constant coherence conditions to compute the com- 
plete set of 0(g 2 ) counterterms. Let us first consider the generic problem of determining 
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the strength of an irrelevant variable from its renormalization group equation, using the 
condition that it can depend on the cutoff only through its perturbative dependence on 
g n . The generic equation for an irrelevant variable can be written, 



/ ]_\n w 

Wn+1= \4J w n~ c w g n + 0(g n ) , (6.55) 

where n w is an integer determined by the transverse dimension of the operator. Note 
that the 0(g%) corrections include corrections to the Hamiltonian coming from diagrams 
that are third-order in the original interaction in Eq. (6.3), as well as corrections coming 
from diagrams that include one original vertex and one counterterm vertex. The fact 
that tadpoles are eliminated when zero-modes are dropped is used here, because an 0(g 2 ) 
counterterm produces O(g^) corrections when zero-modes are removed. The assumption 
that all counterterms are at least 0(g 2 ) must be justified a posteriori. It is easy to see that 
Eq (6.55) implies that an expansion of w n in powers of g n should start at second order. 
Assuming that w n = ug 2 + 0(g^) : and using Eq. (6.36), we find, 



u = — jTp • ( 6 - 56 ) 

Thus, we find that the second-order transformation fixes the strength of all irrelevant 
operators when we insist that they run only with the coupling. To illustrate the conse- 
quences of Eq. (6.56), consider 6114 given by Eq. (6.15). Eq. (6.15) yields the values 
of c w for an infinite number of irrelevant operators. Using Eq. (6.56) we find that the 
Hamiltonian must contain the set of irrelevant operators, 
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(6.57) 



x(l — x) - 



x(l — x) 



4). 



All momenta are defined in Eq. (6.4). There is another counterterm coming from the first 
Hamiltonian diagram in figure 3b, with the momenta in Eq. (6.57) being replaced by the 
momenta of the outgoing particles; as well as additional counterterms from the remaining 
Hamiltonian diagrams in figure 3b. However, one can easily determine the integral form 
of the counterterms in each case. Notice that the integrand in Eq. (6.57) is similar to the 
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integrand in Eq. (6.11). There is a subtraction of the latter integrand with the external 
transverse momenta set to zero to remove the marginal piece, and the step functions are 
altered so that it is intermediate energies above the upper cutoff that are included, rather 
than energies between the upper and lower cutoffs. In addition to having a simple universal 
form, the integral representation of the irrelevant 'counterterms' will prove convenient in 
Section VII, when I group a second-order contribution that includes this counterterm and 
the original interaction in Eq. (6.3) with a third-order contribution from the original 
interaction alone. In diagrammatic terms, I group the one-loop correction that contains 
a vertex counterterm with a two-loop contribution to the running Hamiltonian; and the 
integral representation of the counterterm allows this regrouping to be performed directly 
in the two-loop integrand. 

The coupling constant coherence conditions placed on the irrelevant operators require 
that to 0(g 2 ) they must be invariant under the action of the full transformation. In the 
first part of a transformation one lowers the cutoff, and generates new 0(g 2 ) irrelevant op- 
erators. When these new irrelevant operators are added to the old irrelevant operators, the 
resultant operators must be identical in form to the old irrelevant operators, with the only 
change being the replacement of Ao with Ai. After the scaling part of the transformation 
is completed, the complete set of irrelevant operators returns to exactly its original form; 
in this case, the logarithm in Eq. (6.57) is exactly reproduced. Higher order corrections 
to the irrelevant operators insure that the g\ coefficient runs correctly, so that it is indeed 
the correct running coupling that appears, and they generate new irrelevant operators of 
0(gf l ) and higher. 

As a final example of an irrelevant counterterm, consider the correction to uq resulting 
from Eq. (6.32). All operators in uq are irrelevant, and if uq runs only because the coupling 
in Eq. (6.3) runs, the same reasoning used above implies that uq must contain, 



Comparison of this term with the correction to uq resulting from the second Hamiltonian 
diagram in figure 3c, Eq. (6.32), shows that the second-order tree level counterterms are 
easily determined. 



x 2 "r" x 3 x 6 




(6.58) 
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Return to the calculation of the boson-boson scattering amplitude, and add the above 
irrelevant counterterms to 144. The real part of the scattering amplitude becomes, 
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(6.59) 



It should be obvious from this expression that covariance and cluster decomposition are 
still violated, because of the longitudinal momentum fractions appearing in the logarithms; 
however, these properties can now be restored by marginal operators alone. In other 
words, coupling coherence leads to irrelevant operators that restore covariance and cluster 
decomposition to the irrelevant part of the scattering amplitude. Adding the appropriate 
irrelevant operator contributions to the one-boson to three-boson Green's function yields, 
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(6.60) 



Violations of Lorentz covariance and cluster decomposition in the irrelevant part of this 
Green's function are also removed. 

Let us now consider the generic problem of determining the strength of a marginal 
variable from its renormalization group equation, using the condition that it can depend on 
the cutoff only through its dependence on g n . To simplify the presentation let me simply 
state that we must simultaneously consider dependent marginal operators of O(g^), to 
which I collectively refer as h n , and dependent marginal operators of 0(gf l ), to which I 
collectively refer as j n . The reason that both are required will become apparent below. We 
have already seen that no new marginal operators are produced in the 0(g 2 ) Hamiltonian 
diagrams, so the generic equations for these marginal variables can be written, 



h n +i = K- c h g n h n - d h gl + 0(g%) , 



(6.61) 



jn+i = jn - Cjg n h n - djgl + 0{g^) 



(6.62) 



In principle, there should be terms of 0(g n w n ) in each of these equations; which result 
from second-order (i.e., one-loop) corrections that contain an irrelevant vertex in addition 
to the original vertex in Eq. (6.3). Since the irrelevant counterterms have already been 
uniquely determined to O(g^), I group these one-loop corrections with corrections that 
are third-order (i.e., two-loop) in the original vertex for simplicity. Thus the final terms 
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in these equations, d^gi^ and djg%, result from a sum of two- loop diagrams and irrelevant 
counterterm insertions in one-loop diagrams. It is straightforward, but tedious, to sepa- 
rately display these terms in Eqs. (6.61) and (6.62). It is far more tedious to separately 
compute these one-loop and two-loop diagrams in closed form. Fortunately, there is no 
need to do so. 

By assumption, 

K = agl + (3gl , (6.63) 



jn = <jg 3 n , (6.64) 

with higher order terms being suppressed. It is straightforward to confirm that there are 
two types of solution to Eq. (6.61). If dh = and a ^ 0, we must have 

c h = 2c g ; (6.65) 
where c g is shown in Eq. (6.36). If dh ^ 0, then 

= d h 

2Cg ~C h ^ ' 

(3 is not determined by Eq. (6.61), but it is clear that one must perform a third-order renor- 
malization group calculation and determine dh to fix a, even though this is the coefficient 
of the O(g^) piece of the marginal operator. 

7 is not fixed by Eq. (6.62), but if a marginal operator arises in the third-order 
behavior of the transformation, Eq. (6.62) shows that the strength of this operator is 
0(g3)i as opposed to O(g^), only if 

acj + dj = 0. (6.67) 

This condition must be met if an operator appears in the third-order analysis, but is not 
0{g 1 n ). However, since 7 is not fixed, there is no guarantee that 7 7^ 0; so the appearance of 
the operator in the third-order behavior of the transformation does not insure its appear- 
ance in the Hamiltonian. The important point to observe here is that Eq. (6.67) requires 
a 7^ if dj 0. This means that if we find any new marginal operators in the third-order 
analysis, at least one new marginal operator is O(g^). 
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This generic analysis shows that there are two types of marginal operator that can 
arise with strength of O(g^). One type is explicitly produced by the third-order behavior 
of the transformation, after the irrelevant counterterms have been properly included to 
this order, and this type cannot be missed. The second type of marginal operator does 
not appear in the subtracted third-order analysis; but it must exactly reproduce itself in a 
second-order analysis when combined with the interaction in Eq. (6.3), with the strength 
fixed by the condition = 2c g . I proceed no further with this analysis in this Section. In 
the next Section I show that the logarithmic functions of longitudinal momentum fractions 
required to restore Lorentz covariance and cluster decomposition are indeed solutions to 
Eqs. (6.61) and (6.62). 

VII. Third- and Higher-Order Behavior Near Critical Gaussian Fixed Points 

Third-order behavior of a transformation near the Gaussian fixed point can be com- 
puted using the final terms in Eq. (4.14), and the additional terms required to study 
higher-order behavior are readily computed, with rapidly increasing algebraic complexity. 

There are two issues I want to address in this Section. First, I want to complete the 
calculation of corrections to the marginal part of 114 initiated in the last Section; and show 
that the coupling constant coherence conditions lead to the marginal counterterms required 
to restore Lorentz covariance and cluster decomposition to the four-point functions in 
second-order perturbation theory. Second, I want to discuss the errors one encounters when 
various approximate analyses are performed. The latter discussion is not intended to be 
rigorous or complete. The renormalization group offers the possibility of overcoming some 
serious logical flaws in 'old-fashioned perturbation theory' [112,105,113], a fact that does 
not yet receive sufficient attention in text books. The perturbative renormalization group 
may allow one to effect renormalization so that the perturbative expansions encountered 
at every stage of a calculation involve only small coupling constants. The caveat is that the 
running coupling (s) must remain small over the entire range of scales directly encountered 
in the perturbative portion of the calculation. In 'old-fashioned' perturbation theory the 
bare parameters always diverge, and one must simply follow renormalization recipes [112] 
without being distracted by intermediate expansions in powers of divergent constants. We 
will find that the canonical running coupling constant in the scalar theory remains small if 
it is small in H^°; however, there are logarithms of longitudinal momentum fractions that 
appear in the Hamiltonian, leading to effective couplings that diverge at small longitudinal 
momentum transfer. 

There are severe limitations to how much can be learned from a perturbative renor- 
malization group study of scalar field theory, because it is not asymptotically free and 
there is no possibility of the perturbative analysis being complete [12]. It is not possible to 
construct a trajectory of renormalized Hamiltonians using a perturbative renormalization 
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group if the theory is not asymptotically free. For the examples in this Section, I simply as- 
sume that the coupling constant in H^°, go, is small; in which case Eq. (6.36) reveals that 
it should remain small in every Hamiltonian on a renormalization group trajectory. This 
analysis is naive; however, my interest is the development of the perturbative formalism. 

Let me turn now to the third order corrections to the part of that governs one-boson 
to three-boson transitions. This part of is chosen rather than the two-boson to two- 
boson transition because fewer Hamiltonian diagrams must be computed. The requisite 
Hamiltonian diagrams are shown in figure 5, where I have suppressed the arrows that indi- 
cate the energies appearing in the energy denominators. Each diagram in figure 5 actually 
corresponds to four Hamiltonian diagrams, which differ in their energy denominators and 
in the distribution of cutoffs, as seen in the four third-order terms in Eq. (4.14). As dis- 
cussed in the last Section, I want to group all 0(5H 2 ) terms that are 0(g 3 ) because they 
contain one 0{g) marginal vertex and one 0(g 2 ) irrelevant vertex, with the 0(5H 3 ) terms 
that are 0(g 3 ) because they contain three 0{g) marginal vertices. This comment is most 
easily understood after studying the results below. The bracketed expressions in figure 5 
indicate diagrams in which the momenta of the outgoing particles are permuted, with the 
permutation being listed in the brackets. 

I only want to compute the third-order corrections to the marginal part of U4, so 
I set all incoming and outgoing transverse momenta to zero. Since I am studying the 
critical theory, this means that the initial and final energies are zero, which considerably 
simplifies most of the energy denominators. I do not display spectators because they have 
no direct effect on 114. One can always choose variables in which Sv^ is independent of the 
longitudinal momenta of the spectators, even though cluster decomposition is restored by 
marginal counterterms that explicitly depend on these momenta. 

Let me begin with the Hamiltonian diagrams in figure 5a. There are four diagrams, 
and each makes an identical contribution to the marginal part of 8v±. Using Eq. (4.14) 
and the interaction in Eqs. (6.2)-(6.3), one finds, 



where dk is defined in Eq. (3.4), k~ = k- 1 /k + , and p + is the longitudinal momentum 
entering the loops. This last factor results from the presence of an internal line that is 
not part of a loop. Only one of the four terms in Eq. (4.14) survives here, because only 
the intermediate state that includes the boson loops can be a high energy state, since the 
initial and final states must both be low energy states. The first factor of —1/2 is seen in 
Eq. (4.14); while 3! is a symmetry factor. I have introduced, 




(7.1) 



(k 1 + fc 2 +^3) ®high(ki +^2 ^3 ) ' 
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OM 9h (K-) =<>($:- K-) 0(K- - , (7.2) 
which projects onto states of high energy; and later I also need, 

Q l0W {K-) = e(^-K-), (7.3) 

which projects onto states of low energy. P + is the total longitudinal momentum. I have 
not displayed the low energy cutoffs associated with external lines. 

The integrals in Eq. (7.1) are most easily evaluated by introducing Jacobi coordinates; 
for example, p + = wP + , k^ = xyp + , k^i = y/wAo(xs + r), k£ = (1 — x)yp + , k _L 2 = 
a/m;Ao((1 — x)s — r), k% = (1 — y)p + , = — y/wAos. In these coordinates, 



l l 

7 3 [ d 2 r d 2 s 



Sv ' M = - g ulT6^^J dx J dy 



where now, 



x(l - x)y(l - y) 
o o 

( 1-2 4 ^ ^ ~ 2 O ( r2 4 ^ 1 
\yx(l-x) y{l-y)) high \yx(l - x) y(l-y)/' 

(7.4) 



M7) = f7) J (7-')' (7 - 5) 



and, 



MtH("-t)- (7 ' 6) 

I have introduced rj = A^/Aq because this ratio appears repeatedly. The appropriate 
definition of Qhigh an d ©Zow is always apparent from context. All integrals are evaluated 
using Jacobi coordinates, so the definitions in Eqs. (7.5) and (7.6) are always used to 
explicitly evaluate a correction. 

The remaining integrals are readily completed, and one finds that each Hamiltonian 
diagram in figure 5a contributes, 
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5v 4 m = 



24 (16tt 2 ) 2 V?? 




(7.7) 



Of course, the entire correction 5v 4 is extremely complicated, but the marginal part of 5v4 
is simple. The rescaling does not change this marginal term, so Eq. (7.7) immediately 
yields the correction to the marginal part of u\ from each diagram in figure 5a. 

Next consider the diagrams in figure 5b. The first diagram in figure 5b leads to, 



The first two of the four third-order terms in Eq. (4.14) combine to give the first term in 
the integrand. The third term in the integrand is zero because it is not possible for the 
energy of the first intermediate state in figure 5b to have higher energy than the second 
intermediate state. I do not display terms that vanish in this manner below. 

While it is possible to evaluate this integral, it is convenient to group it with the one- 
loop diagrams that contain all irrelevant counterterms that result from sub-diagrams (e.g., 
nested loops) in figure 5b. In figure 6a I show the first diagram in figure 5b added to a 
new one-loop diagram in which there is a new four-boson vertex. The new vertex is a 
sum of irrelevant operators that must be added to the Hamiltonian to satisfy the coupling 
constant coherence conditions, and figure 6b shows the original diagram that led to their 
addition. The derivation of irrelevant operators is discussed in the last Section. Using an 
integral representation for the appropriate irrelevant counterterms, the one-loop correction 
to the Hamiltonian shown in figure 6a is, 




(7.8) 
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^(16ir 3 )d 3 (p 2 -k 1 -k 2 -h) 

(k\ I ^2 ' ^3 ) (^1 ^2 ^4 ) ® super (^1 "I - ^2 ^3 ) 

- (167V 3 )5 2 (k ± 2 + k ± 3 )8{p+ -k+- k+ - k+) + fc^")- 1 e^perC^r + fc 2 ) j • 

(7.9) 

The first two terms in the integrand have identical energy denominators to the second-order 
Hamiltonian diagram in figure 6b that leads to these counterterms, while the third term 
cancels the marginal part of the counterterm and insures that it is composed of irrelevant 
operators only. I have introduced a new cutoff function, 



Q super {K-) = e{K- --^j . (7.10) 

This cutoff is not really associated with an intermediate state, because it is part of the 
integral representation of an irrelevant operator that is added to the Hamiltonian, not 
part of an operator that is induced by a transformation. However, one can think of an 
intermediate state in which the energy lies above all cutoffs; and say that the counterterm 
directly implements a complete set of irrelevant interactions that would have been provided 
by the exchange of particles whose energy is above the cutoffs placed on the Hamiltonian. 

I want to emphasize that the four-boson vertex nested in the one-loop diagram in figure 
6a, while complicated, is uniquely determined by coupling coherence. If this vertex is not 
in the Hamiltonian, there are variables that explicitly run with the cutoff other than the 
coupling in Eq. (6.3). 

To proceed we need to define Jacobi coordinates that satisfy the delta function con- 
straints, and then combine the integrals remaining in Eq. (7.8) with those in Eq. (7.9). 
The calculation from this point is complicated only because of tedious calculus, which I 
do not detail. Perhaps the most difficult part of the calculation is keeping track of the 
various step function cutoffs, and finding how to regroup terms at appropriate stages of 
the calculation so that nested integrals lead to simple analytical results. These problems 
make the calculation more difficult than a simple two-loop Feynman diagram calculation, 
but this should be no surprise. 

If we let P2 — w(p~l — p±), so that p% = (1 — w)(pl — p\ ), the complete result from 
the diagrams in figure 6a is, 
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This result needs to be added to those from the remaining diagrams in figure 5b, with the 
one-loop corrections that are analogous to Eq. (7.9) being added. The second diagram in 
figure 5b leads to a correction that is identical in form to Eq. (7.11), but with w and 1 — w 
interchanged. Letting p 2 = xp\ , = W\-, an d v\ = Z V\ \ the complete sum of terms in 
figure 5b, with the appropriate one-loop counterterm diagrams added, is, 
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(7.12) 

The first diagram in figure 5c is redisplayed in figure 7a, with the one-loop correction 
that contains the appropriate irrelevant counterterm. Figure 7b displays the sub-diagram 
from which the irrelevant counterterm results, and one sees that in this case the countert- 
erm is part of uq. The sum of both diagrams is, 
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+ 2^1" + ^ + [(K + k 2 + KY 1 + (K - KV 1 

®high(k\ + k 2 + ^3 ) © ' super {k 1 + k 2 + ^4 ) ^ • 

(7.13) 

The first two terms in the integrand result from the third-order corrections to the Hamil- 
tonian in Eq. (4.14), while the last two terms result from the second-order corrections in 
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Eq. (4.14) with one of the interactions being an irrelevant operator. The steps required to 
evaluate this integral are identical to those required above, and the result is, 



SV4M = ~ 2(16tt 2 ) 2 K^) /n(1 " X) ' (7 ' 14) 

where = xp\ . 

The complete set of diagrams in figure 5c, with the accompanying one-loop corrections, 
yield, 
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where I have again chosen = xpf, p% = yp\ , and p\ = zp\ . 

The first diagram in figure 5d is redisplayed in figure 8a, along with two one-loop 
corrections with which it must be grouped. In this case there are two sub-diagrams that 
lead to irrelevant counterterms, as shown in figures 8b and 8c. After using coupling 
coherence to uniquely determine the irrelevant vertices in the one-loop diagrams, one finds 
that the diagrams in figure 8a yield, 
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The complete set of diagrams in figure 5d, with the accompanying one-loop corrections, 
yield 
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The final set of two- loop contributions, along with the accompanying 0(g 3 ) one-loop 
contributions, to the marginal part of Sv^ are shown in figure 9a. These diagrams yield 
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At this point we have all O(g^) contributions to the renormalization of the marginal 
one-boson to three-boson part of U4, except for one- loop contributions that involve one of 
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the 0(g 2 ) marginal operators that we must find using coupling coherence. For the benefit 
of later calculations it is convenient to list the complete result using coordinates in which 
= X{P + . The complete set of two-loop contributions to the marginal one-boson to 
three-boson part of Su^, along with all O(g^) one-loop contributions that result from the 
0(g 2 ) irrelevant operators determined by coupling coherence, yields: 




This result includes four separate functions of longitudinal momenta that must be 
considered in the renormalization group analysis, each symmetric under the interchange 
of outgoing momenta. Any one or all of these may appear in 114 at 0(g 2 ). At least one of 
them must appear at this order if the Hamiltonian satisfies the coupling constant coherence 
conditions, as was demonstrated at the end of Section VI, so there are several possibilities 
that should be studied. 



A similar calculation that is slightly more tedious leads to the two-loop contributions 
to the two-boson to two-boson marginal part of 6114. The complete set of two- loop contri- 
butions to this part of 8U4, along with all 0(g 3 ) one-loop contributions that result from 
the 0(g 2 ) irrelevant operators determined by coupling coherence, yields: 
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Here I use the momenta in figure 3b, again choosing = X{P + . There are four functions 
of longitudinal momenta that appear in this correction to the marginal operator, all of 
which must be included in the renormalization group analysis. 

A complete renormalization group analysis of the marginal operator requires us to at 
least introduce all of the functions of longitudinal momentum fractions appearing above in 
the marginal part of which I indicate as g(x\, x 2 , X3, X4). We should distinguish be- 
tween the two-boson to two-boson and the one-boson to three-boson parts of this operator. 
The three-boson to one-boson interaction is not independent because the Hamiltonian is 
Hermitian. Thus, we are led to consider, 
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(7.21) 
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I have simply assumed that the constant, g, appearing in each of these terms is the same. 
This issue does not need to be resolved until one computes the O(g^) corrections to the 
renormalization group equation for the running coupling, so I do not pursue it further. I 
have also simply assumed that some of the operators appearing in g^ 2 ^ 2 ^ and g*- 1 ^ 3 ** are 
0(g 2 ), with their couplings being foW; while other operators are assumed to be at least 
0(g 3 ), with their couplings being j'W. These assumptions are justified a posteriori, and I 
do not try to prove that this solution is unique. I believe that the solution is unique, but 
a proof exceeds my patience. 

To proceed further we must generalize Eqs. (6.61) and (6.62) to allow for additional 
operators. The complete equations are, 



h 



(k) 

n+l 



1=1 



(7.23) 



J 



;(*) _„■(*) 



A k) - E c t l) 9n $ } - df gl + 0{gi) 



n+l 



1=1 



Eqs. (7.19) and (7.20) indicate, 
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The calculation of the various coefficients c^'^ and c^ k ' 1 ^ requires us to evaluate a set of 

second-order Hamiltonian diagrams in which one vertex comes from the constant g, and 

(1 2) 

the second vertex comes from the appropriate marginal operator. Thus, to compute c- , 
for example, we need to include the vertex corresponding to the operator whose strength 
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is and determine the part of the resultant one-loop integral that leads to a function 
of longitudinal momenta that exactly matches the function in the operator with strength 
These calculations are straightforward, and the results are, 
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All others are zero. 

When one solves Eqs. (7.23) and (7.24) using these coefficients, it is found that each 



These results correspond exactly with the counterterms required to restore Lorentz co- 
variance and cluster decomposition to the boson-boson scattering amplitude and to the 
one-boson to three-boson Green's function. Thus, the coupling constant coherence con- 
ditions lead to a solution of the complete renormalization group equations that yields 
covariant results with cluster decomposition despite the fact that the cutoff violates these 
conditions. 

An exact renormalization group analysis apparently leads to cutoff-independent results 
with a cutoff Hamiltonian, and coupling constant coherence apparently allows one to find 
Hamiltonians that produce Lorentz covariant results. This is useless, however, unless one 
can make approximations with bounded errors. Moreover, if one simply wants to use 
renormalization group improved perturbation theory in powers of the canonical coupling, 
it should be clear from the above calculations that the light-front renormalization group 
is not the simplest tool at one's disposal. 

It is not obvious how one should estimate 'errors'. One might consider fixing the initial 
Hamiltonian, and measuring the errors by computing differences in the results produced by 
the final Hamiltonian on a trajectory in comparison to those from the initial Hamiltonian. 
This is not a good measure, because in general we do not know the appropriate initial 
Hamiltonian nor can we compute with it, and one of the goals of the renormalization 
group is to formulate physical problems so that it is never necessary to explicitly deal 
with the initial Hamiltonian. In order to produce a meaningful discussion of errors, we 
should consider how a 'typical' renormalization group calculation proceeds, and study how 
observables change as the approximation is systematically improved. This discussion was 
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initiated at the end of Section II, and some basic issues were discussed in Section VI; 
however, the results computed in this Section have a dramatic effect on the analysis of 
errors in the light-front renormalization group. 

In Wilson's perturbative Euclidean renormalization group calculations, [12] it is nec- 
essary to fix all relevant and marginal couplings at the lowest cutoff, and any irrelevant 
couplings at an upper cutoff. The upper cutoff is not chosen to be infinite, but should be 
sufficiently large that all final results are insensitive to the boundary values chosen for the 
irrelevant variables. It should not be so large that intolerable errors accumulate over the 
trajectory. One must approximate the marginal couplings over the entire trajectory, and 
use an iterative algorithm to compute the trajectory. The output is the irrelevant couplings 
at the lower cutoff, because all relevant and marginal couplings at this cutoff must be in- 
put. This is an extremely powerful procedure when there are a finite number of relevant 
and marginal couplings, because it allows one to obtain an accurate approximation to the 
endpoint of a complete renormalized trajectory by inputting a finite number of boundary 
values [13]. However, we have seen that the light-front renormalization group requires an 
infinite number of relevant and marginal operators. This means that a light-front renor- 
malization group calculation requires a new type of approximation not considered in the 
Euclidean renormalization group. One must approximate the boundary conditions placed 
on the relevant and marginal variables at the lowest cutoff. 

There are three obvious approximations one should consider. First, one can approxi- 
mate the transformation itself by dropping terms at a given order in 5H . Second, one can 
approximate the trajectory by discarding specific operators; e.g., all or some of the irrel- 
evant operators. Third, one can employ coupling coherence to compute the Hamiltonian 
to a given order in the running canonical coupling and drop higher orders. One can also 
employ combinations of these approximations. 

I assume that Feynman perturbation theory is accurate, and estimate the errors from 
each approximation by identifying the order in Feynman perturbation theory at which an 
error first arises and then discussing the magnitude of this error. Such an analysis is of 
limited use if Feynman perturbation theory is not valid for the computation of low energy 
observables, which is exactly the case for which the light-front renormalization group is 
being developed; but this should give a rough guide to the problems we should study. Let 
us begin by approximating the transformation by dropping terms starting at some given 
order in 5H. The simplest comparison we can make that is of any interest is between the 
0(5H) analysis and the 0(5H 2 ) analysis. The linear analysis is trivial to complete. We 
set the irrelevant variables to arbitrary strengths, and find that they go to zero at the lower 
cutoff as the upper cutoff is taken to infinity. This leads to errors in Feynman perturbation 
theory that are 0(g 2 ) 1 at least in the irrelevant variables. This means for example, that 
the real part of the boson-boson scattering amplitude, computed in Eq. (6.53), contains 
errors of at least 0{g 2 M 2 /A 2 V ). As long as one is studying scattering for states whose 
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invariant mass is much less than the cutoff, these errors are small. This assumes that one 
inputs the correct marginal and relevant variables at the lower cutoff, Ay\y\ For example, 
if one does not input the marginal operators computed above to 0(g 2 ) using coupling 
coherence, there are logarithmic errors shown in Eq. (6.59) that are arbitrarily large. As 
w — > 0, these errors diverge like g 2 log(w); and as the longitudinal momenta of the outgoing 
bosons approaches those of the incoming bosons, there are comparable logarithmic errors. 
For small g these errors are suppressed relative to the leading term by one power of g, but 
the relative error can be arbitrarily large for scattering states of any invariant mass. 

The lesson here is that arbitrarily large errors arise when functions of longitudinal 
momenta diverge. Longitudinal divergences forced us to discard several candidate trans- 
formations, and we are seeing that the invariant-mass transformations do not completely 
control the spectrum of the longitudinal operators. This problem is so severe that its so- 
lution may require a completely different renormalization group than has been developed 
in this paper. 

If we improve the analysis by keeping 0(SH 2 ) corrections to the trajectory, we auto- 
matically generate the correct irrelevant operators to 0(g 2 ) 1 even if we do not include the 
correct 0(g 2 ) relevant and marginal counterterms at the lower cutoff. This means that 
we obtain Eq. (6.59) instead of Eq. (6.53) for the real part of the boson-boson scattering 
amplitude, for example. If the correct marginal operators are included in the Hamiltonian 
at the lower cutoff, we obtain Feynman results to 0(g 2 ). Therefore, if we approximate the 
transformation at 0(5H 2 ), and in addition we make a perturbative expansion in terms 
of the canonical coupling to 0{g 2 ), and we impose the correct boundary conditions on 
the marginal and relevant operators at the lower cutoff, which means fixing functions of 
longitudinal momenta, we obtain the Feynman results to 0(g 2 ). Of course, we do not need 
to make the additional perturbative expansion in powers of g, but we must make some 
additional approximations, because the second-order transformation generates an infinite 
number of vertices and each contains an entire function of momenta. 

When the transformation is approximated to 0(5H 2 ) : U2 does not affect any of the 
other functions in the Hamiltonian, as discussed in Section VI. This means that the relevant 
operators produced by the transformation have no effect on the other operators and can be 
studied separately. There is a large host of additional approximations one might consider 
that are not perturbative in the canonical coupling and may be of interest. All of them 
produce 'errors' at O(g^); so if Feynman perturbation theory in powers of the canonical 
coupling is accurate, most of the additional approximations one can make offer little or 
no improvement to the perturbative approximation discussed in the last paragraph. One 
interesting approximation is to keep only the marginal part of «4, and complete an analysis 
that is a generalization of the analysis leading up to Eq. (6.36) for the running canonical 
coupling. Since the boundary condition on this marginal operator includes functions of 
longitudinal momenta, the analysis in Section VI must be generalized to include such 
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functions. This analysis reduces to the study of coupled one-dimensional integral equations, 
because the dependence of the marginal operator on transverse momentum is fixed, making 
all transverse integrals the same. Perhaps the most important question one can address 
with such an investigation is whether functions of longitudinal momenta with nonintegrable 
singularities are generated, since the boundary functions include logarithmic divergences. 
I believe that it is relatively easy to show that in the 0(5H 2 ) analysis no such singularities 
arise. One can solve the coupled equations for the marginal part of 114 'backwards', towards 
larger cutoffs, using the functions in Eqs. (7.21) and (7.22) as a boundary condition. In 
each step one simply convolutes the functions produced from the previous step, which 
means that after any finite number of steps one effectively considers a multidimensional 
convolution of logarithms, and these are always finite. 

Any estimation of errors requires one to place bounds on the longitudinal functions. I 
have ignored this issue, except where the integrals over longitudinal momentum produce 
divergences. I have shown that this does not happen to lowest orders in the running canon- 
ical coupling when one uses the invariant-mass renormalization group; however, I have not 
shown that this does not happen in higher orders. Even if a perturbative expansion in 
powers of the canonical coupling is finite to every order, it is possible for some of the 
running variables in the invariant-mass renormalization group to diverge; as long as all 
such divergences cancel against one another when one re-expresses any physical result as 
a power series in terms of the canonical coupling. If one performs only an expansion in 
powers of the canonical coupling, and makes no further approximations, it is even possible 
to use the transformations that run cutoffs on the transverse and longitudinal momenta 
of individual particles. One should be able to use coupling coherence in this case to again 
obtain the Feynman results to any order one desires. The problem is that these results 
depend on huge cancellations that must be precisely maintained. I have not found any rea- 
son to believe that nonintegrable singularities arise in the invariant-mass analysis; however, 
this is far from satisfactory. 

Suppose one uses coupling coherence and computes the Hamiltonian exactly to a given 
order in g, and then computes low energy results using this Hamiltonian. How large are 
the errors in perturbation theory? Obviously the results are exactly those of Feynman 
perturbation theory up to the order to which the Hamiltonian is computed, but beyond 
this order one encounters errors. For example, the canonical <^> 4 Hamiltonian is accurate to 
0(g); but if we compute the boson-boson scattering amplitude to 0(g 2 ) we obtain the result 
in Eq. (6.53), with the same errors discussed above. The errors in perturbation theory 
are of 0(g 2 ), but even when g is small the errors can be arbitrarily large. This same type 
of error arises no matter how many orders in g are included in the calculation of H. This 
problem may best be understood by thinking of the coupling as running with longitudinal 
momentum. The second order corrections to the Hamiltonian in Eqs. (7.21) and (7.22) 
show that the coupling decreases as the longitudinal momentum transferred through the 
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vertex decreases. In fact, according to the perturbative analysis, the coupling actually 
becomes negative at sufficiently small longitudinal momentum transfer. The perturbative 
analysis breaks down when such logarithmic divergences arise, and one must find a method 
of re-summing these corrections. In this case, we are finding that even if the coupling is 
small for all longitudinal momentum transfers, we cannot expand the coupling for one 
momentum transfer in powers of the coupling at some drastically different momentum 
transfer. This is the sort of problem the Euclidean renormalization group manages to 
avoid, but the invariant-mass renormalization group does not treat all components of the 
momentum on an equal footing, and this is a price that must be paid. 

VIII. Conclusion 

After defining several renormalization groups that might be of interest in the study of 
light-front field theory, I have found the Gaussian fixed points and completed the linear 
analysis about the Gaussian fixed point of interest for relativistic field theory. The subse- 
quent second-order analyses of these transformations show that a perturbative expansion 
of a transformation about the Gaussian fixed point, in which some of the irrelevant opera- 
tors are discarded, can only converge for transformations that remove states of higher 'free' 
energy than all states retained, where the free energy is determined by the Gaussian fixed 
point. This constraint naturally leads to transformations that employ invariant-mass cut- 
offs. The linear analysis of invariant-mass transformations reveals that there are an infinite 
number of relevant and marginal variables, because functions of longitudinal momenta do 
not affect the invariant-mass scaling dimension of operators. 

While the linear analysis of a light-front transformation about the Gaussian fixed 
point is readily completed for arbitrary Hamiltonians, the second-order analysis has a 
complicated dependence on the interactions. One can easily write a general expression for 
the transformed Hamiltonian using Eq. (4.14). The relevant and marginal operators in 
U2 do not affect other operators in the second-order light-front analysis when zero-modes 
are dropped, and this simplifies the analysis considerably. The only marginal interaction 
that directly enters the second-order analysis is the marginal part of U4. If all irrelevant 
operators are dropped, only the relevant and marginal parts of U2 and U4 survive. Since the 
marginal part of w.4 contains an arbitrary function of longitudinal momenta, the second- 
order correction to the marginal part of W4 involves a new function of longitudinal momenta 
and a complete analysis requires one to compute trajectories of such functions, in general. 
While an analysis that allows arbitrary functions is not extremely complicated, the simplest 
example is the trajectory generated when the initial value of the marginal part of w.4 is a 
constant, as suggested by canonical field theory. In this case the second-order correction 
leads to a new constant and not to a function of longitudinal momentum. This case was 
studied in detail, and it was shown that the canonical coupling decreases as the cutoff is 
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lowered. This analysis is readily improved by retaining leading irrelevant operators, and 
allowing functions to appear in the marginal operator, but this was not done in this paper. 

The invariant-mass cutoff violates explicit Lorentz covariance and cluster decomposi- 
tion, so the Hamiltonians one must investigate do not display these properties manifestly. 
If one uses the canonical </> 4 Hamiltonian with an invariant-mass cutoff to compute Green's 
functions, it is not surprising to find that covariance and cluster decomposition are vio- 
lated. If one arbitrarily chooses a more complicated renormalized Hamiltonian, there is no 
reason to expect that these properties are restored. One must select the correct functions 
of longitudinal momentum in the marginal and relevant operators at the lowest cutoff, 
and use the perturbative renormalization group to estimate the irrelevant operators at 
the lowest cutoff, to restore Lorentz covariance and cluster decomposition to observables. 
While it is possible to adjust these functions until correct results are obtained in pertur- 
bation theory, I showed that one can also use coupling constant coherence [39] to restore 
covariance and cluster decomposition. By insisting that the canonical coupling is the only 
variable that explicitly runs with the cutoff, and that all other variables depend on the 
cutoff through perturbative dependence on the canonical coupling, one uniquely fixes the 
relevant and irrelevant operators at the lowest cutoff to 0(g 2 ). The marginal variables are 
not uniquely fixed until one completes a third-order calculation, and the correct results are 
not obtained unless 0(5H 3 ) terms are kept in the transformation. When these third-order 
terms are retained, it is possible to use coupling coherence to fix the functions appearing in 
the marginal part of 114 to 0(g 2 ), and these functions are exactly those required to restore 
Lorentz covariance and cluster decomposition to observables computed in second-order 
perturbation theory. 

The appearance of entire functions in the relevant and marginal operators severely 
complicates the development of a perturbative light-front renormalization group. While 
it is encouraging to find that coupling coherence apparently fixes these functions without 
direct reference to Lorentz covariance and cluster decomposition in the output observables, 
the light-front renormalization group does not offer a convenient method for performing 
perturbative calculations. The Euclidean renormalization group is far more convenient 
for perturbative calculations, because there are a small number of relevant and marginal 
operators [12]. In gauge theories Euclidean cutoffs typically violate gauge invariance, and 
thereby covariance (with the lattice providing an example of how one can take advantage 
of irrelevant operators to maintain gauge-invariance [114]); but one should be able to 
use coupling coherence to restore these properties in a Euclidean renormalization group 
analysis, as was done for the light-front renormalization group in this paper. Thus, the 
appearance of functions in the relevant and marginal operators leads to computational 
complication of a perturbative light-front analysis in comparison to a Euclidean analysis 
at best. 

Many of the basic problems one encounters in trying to apply the light-front renormal- 
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ization group to QCD have been discussed in this paper; however, several very important 
problems have been completely avoided. The most immediate problems result from di- 
vergences in QCD that are not regulated by the invariant-mass cutoff. In second-order 
perturbation theory, one finds that ui for both quarks and gluons diverges if one uses 
only an invariant-mass cutoff. There are logarithmic divergences that remain in one-loop 
longitudinal momentum integrals after the invariant-mass cutoff is imposed, and these re- 
quire an additional regulator. These infrared longitudinal divergences do not appear in 
the one-loop corrections to the three-particle and four-particle interactions in QCD if one 
uses the canonical Hamiltonian, but they are avoided only because there are precise can- 
cellations maintained by the canonical Hamiltonian [116,59]. These cancellations require 
sectors of Fock space with different parton content to contribute with relative strengths 
given by perturbation theory, and there is no reason to believe that such cancellations are 
maintained nonperturbatively. 

In addition to new infrared divergences, a perturbative analysis using coupling co- 
herence in QCD may reveal that the canonical coupling increases in strength for small 
longitudinal momentum transfer, rather than decreasing as found in the scalar theory. 
This would mean that one cannot simply re-sum the logarithmic corrections to the canon- 
ical coupling and produce a reasonable coupling that runs with longitudinal momentum 
transfer, because the coupling one obtains in this fashion diverges for a finite longitudinal 
momentum transfer. 

At this point it is not even obvious that these problems are unwelcome. When dealing 
with QCD one must balance seemingly contradictory goals. Because of asymptotic free- 
dom a perturbative renormalization group analysis may enable one to lower the cutoff on 
invariant masses to a few GeV. At the same time, the theory should confine quarks and 
gluons, and confinement certainly cannot result from any interaction that can be treated 
perturbatively at all energies. Ideally only a few operators are required to accurately ap- 
proximate confinement effects. If a large number of operators are required it is probably 
quite difficult to find them in an approximate analysis. Moreover, hopefully these interac- 
tions do not change particle number. If operators that change particle number diverge in 
strength, it is difficult to approximate the states. A renormalization group analysis that 
produces an accurate effective QCD Hamiltonian with a cutoff of a few GeV must either 
be able to treat the confining interactions nonperturbatively when eliminating high energy 
states, or these interactions must be perturbative in strength for high energy states and 
diverge in strength only for low energy states. 

Another problem that has been avoided in this paper is symmetry breaking. Unless 
the cutoff violates the symmetry of interest, one never finds symmetry breaking terms in 
a perturbative analysis. They must simply be added to the Hamiltonian by hand, and 
their strength must be tuned to reproduce an observable that is computed using the final 
Hamiltonian, or one must use coupling coherence. It is also possible for terms that do not 
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violate any symmetry to arise in this fashion. Any effective interaction associated with 
a vacuum condensate, or more generally any effective interaction that does not depend 
on the canonical coupling analytically, must be introduced directly in the Hamiltonian if 
one performs a light-front calculation with the zero modes removed, because the vacuum 
is forced to be trivial in this case. Once again, there are two attitudes one can take 
concerning vacuum-induced interactions. One can emphasize the fact that the light-front 
offers little or no advantage to anyone interested in solving the QCD vacuum problem; or 
one can emphasize the fact that light-front field theory may offer anyone who is primarily 
interested in building hadrons a way to reformulate the vacuum problem [117]. 

The best way to check whether a cutoff QCD Hamiltonian resulting from a perturbative 
light-front renormalization group analysis is even crudely accurate is to diagonalize it and 
determine whether it produces even a crude description of low-lying hadrons. One can 
use a trial wave function analysis once the Hamiltonian is selected, which is a powerful 
nonperturbative tool that is not usually available in a field theoretic analysis. It seems 
highly unlikely that a perturbative QCD Hamiltonian will produce reasonable results, and 
one knows that it will not produce the mass splittings associated with chiral symmetry 
breaking if zero modes are discarded. It is necessary to introduce new operators to produce 
symmetry breaking in the spectrum, and it is almost certainly necessary to introduce 
additional vacuum-induced interactions to produce a reasonable description of the mass 
splittings associated with confinement. 

Many of the problems one expects to encounter in such an endeavor are familiar from 
the histories of the constituent quark model [118] and lattice QCD [21]. The constituent 
quark model indicates that it should be possible to obtain accurate models of all low-lying 
hadrons with cutoff Hamiltonians, and a few degrees of freedom. A constituent picture 
may arise naturally from light-front QCD if the parton number-conserving interactions in 
the cutoff Hamiltonians produce mass gaps between sectors with different parton content 
that are reasonably large in comparison to the cutoff. The easiest Hamiltonians to use in 
producing simple constituent hadrons may employ fairly low invariant-mass cutoffs, making 
them analogous to strong-coupling lattice actions [114]. To obtain reasonable results with 
'strong-coupling' Hamiltonians, it may be necessary to simply tune the strength of various 
operators by hand, using phenomenology as a guide when selecting candidate operators. 
It is essential to show that if one makes such uncontrolled approximations it is possible 
to reproduce low energy hadronic observables. On the other hand, if one wants to work 
with large cutoffs for which the canonical coupling remains small, it may be possible to 
use somewhat simpler Hamiltonians; but then the Fock space required to diagonalize the 
Hamiltonian is large and there is no reason to expect that a few constituents will yield 
even crude approximations. Hopefully an intermediate ground exists in which low-lying 
hadrons are adequately approximated as few parton states, while the Hamiltonian is not 
forced to be absurdly complicated by the low cutoff. 
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The light-front renormalization group may eventually lead to the solution of some of 
the most important and interesting problems encountered in the study of low-energy QCD. 
At this point it merely offers a new perspective on these problems. This perspective differs 
radically from those offered by Euclidean field theory, and it is not surprising that difficult 
new challenges appear. While the primary accomplishment of this article is to illustrate 
how coupling coherence may allow one to obtain unique answers in perturbation theory 
without maintaining manifest covariance and gauge invariance, I hope that this article 
also elucidates some of the challenges we must meet to use light-front field theory to solve 
QCD. 
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Appendix A: Canonical Light-Front Scalar Field Theory 



Canonical light-front scalar field theory is discussed by Chang, Root and Yan [48,49], 
who derive many of the results below. This Appendix is not intended as an introduction to 
light-front field theory, but merely collects some useful formulae and establishes notation. 
In this paper I choose the light-front time variable to be 

x + = x° + x 3 , (A.l) 

and the light-front longitudinal space variable to be 

x- = x°-x s . (A.2) 

With these choices the scalar product is 

a ■ b = a% = \a + b~ + \a~b + - a 1 ■ , (A3) 
the derivative operators are 

and the four-dimensional volume element is 

d A x = ^dx + dx~d 2 x ± . (A5) 
The Lagrangian is independent of the choice of variables, 

C = \d^d^ - ^ V - j ( A6) 

for example. The commutation relation for the boson field is 



[<P(x + , x-, x ± ),d + <P(x + , IT, r)] = it\x - y). (A7) 

In order to derive the Hamiltonian and other Poincare generators, one typically begins 
with the energy-momentum tensor [48,49]. While it is certainly possible to derive a formal, 
ill-defined expression for the complete tensor, I merely list the canonical Hamiltonian, 
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dx d 2 x~ L 



(A.8) 



Eq. (A.8) provides a formal definition of the Hamiltonian. The boson field can be 
expanded in terms of a free particle basis, 



a(q)e- iq - x + a\q)e iq - x 



with 



If we use 
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[a( q ),a}(q')] = l^\ + ^(q-q'). 



(A9) 



(A10) 



4>(x) 



d A k 

(2^ 



5{k 2 - ii 2 ) (j){k) e 
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(All) 



we find that, 



(f>{q) = a{q) [q + > 0] , <f>{q) = -a\-q) [q + < 0] 



(A12) 



In the remainder of this Appendix I simply assume that all longitudinal momenta are 
greater than zero. 

The free part of the Hamiltonian in Eq. (A.8) is, 



H = I dx-d 2 x L ^ <f)(x) [ -d^ + n z ) <j)(x) 



The interaction term is more complicated and I do not expand it in terms of creation and 
annihilation operators. Fock space eigenstates of the free Hamiltonian are, 



191,92,...) = a\qi)a\q 2 ) ■ ■ ■ |0) 



(A14) 



with the normalization being, 



{k\q) = 167T 3 q + S(k + - q+^ik 1 - g x ) 



(A15) 



Completeness implies that 
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, rt w rtl fdq + d 2 q ± l , /l 1 [ dq+d 2 q L f dk + d 2 k L . IW 
so we can write the free Hamiltonian as 



#0 = 
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+ 2! 

+ • • 



dq + d q 
167r 3 g+ 
dq + d 2 q- L 



q 12 + [i 2 



q~ 



\q){q\ 



dk+d 2 k~ 



167r 3 g+ J 16n s k + 



q ±2 + [i 2 



+ 



k ±2 + fi 2 
k+ 



\q,k)(q,k\ 



(A.17) 



In order to provide further orientation, let me consider the problem of computing the 
connected Green's functions for this theory. This problem in equal-time field theory is 
discussed in many textbooks, and has been reviewed for light-front field theory [83-85]. To 
find the Green's functions of a theory, consider the overlap of a state |i(0)) at light-front 
time x + = with a second state \f(r)) at time x + = r. We split the Hamiltonian into a 
free part Hq and an interaction V, and find that 



i 



(f(r)m) = (l^)5\P f -P i )G(f^r) = (f\e-^l 2 \ 
■ J ^e-^/ 2 (16 7 r 3 )5 3 (P / - Pi) G(/,i;e). 



— i 



This definition differs slightly from that given in some other places. It is then straightfor- 
ward to demonstrate that 



(K^)6HP, - Pi) G(/,i;e) = (/I e _ g + j0+ 1 



e- H + i0 + e - H + i0+ e - H + i0+ 
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+ V V + • • -i ). 

e - H + i0+ e - H + i0 + e - H + i0 + 

(A19) 

Operator products are evaluated by inserting a complete set of eigenstates of Ho between 
interactions, using Eq. (A. 16), and using 
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H Q \qi,q2,...) = E?, 7 



91,92,-} , 




(A20) 



to replace operators occurring in the denominators with c-numbers. 

Divergences arise from high energy states that are created and annihilated by adja- 
cent Vs, for example. All divergences in perturbation theory come from intermediate 
states (internal lines) that have large free energy. The free energy is given by Eq. (A. 20), 
so divergences occur in diagrams containing internal lines that carry large perpendicu- 
lar momentum ('ultraviolet' divergences) and/or small longitudinal momentum ('infrared' 
divergences) . 

Given a light-front Hamiltonian, H = Hq + V, one can determine the rules for con- 
structing time-ordered perturbation theory diagrams. These diagrams are similar to the 
Hamiltonian diagrams used in the text, but there are important differences related to the 
energy denominators. The diagrammatic rules allow us to evaluate all terms that occur in 
the expansion for the Green's functions given in eq. (A. 19). For ^> 4 canonical field theory 
the diagrammatic rules for time-ordered connected Green's functions are: 

(1) Draw all allowed time-ordered diagrams with the quantum numbers of the specified 
initial and final states on the external legs. Assign a separate momentum k^ to each 
internal and external line, setting k~ = {k\ + [i 2 )/k + for each line. The momenta 
are directed in the direction time flows. 

(2) For each intermediate state there is a factor (e — ^ k~ + z0+) \ where the sum is 
over all particles in the intermediate state. 

(3) Integrate J d 1 fc 6 ' 7r 3 fc + ± 0(k + ) for each internal momentum. 

(4) For each vertex associate a factor of A &{Kf n — K^ ut ) S 2 (K^ — K^ ut ), where Kf n is 
the sum of momenta entering a vertex, etc. 

(5) Multiply the contribution of each time-ordered diagram by a symmetry factor 1/S, 
where S is the order of the permutation group of the internal lines and vertices 
leaving the diagram unchanged with the external lines fixed. 

To obtain the Green's function defined in Eq. (A. 19), propagators for the incoming and 
outgoing states must be added, and one must divide by an overall factor of (167r 3 )5 3 (Py — 



Pi)- 
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Appendix B: Example Calculation of u 2 in the Effective Hamiltonian 



In this Appendix the second-order change in the effective Hamiltonian shown in figure 
3a is computed. Let the first term in the Hamiltonian be 

H = j dqx dq 2 (16tt 3 )5 3 (qi - q 2 ) u 2 {-q u q 2 ) a\ qi )a{q 2 ) + ■■■. (B.l) 
Then the matrix element of this operator between single particle states is 

(p'\H\p) = (0| a(p')Ha\p) |0) = (16n 3 )5 3 (p - p) u 2 (-p,p) . (B.2) 

Thus, we easily determine u 2 from the matrix element. It is easy to compute matrix 
elements between other states. Eq. (4.14) gives us the matrix elements of the effective 
Hamiltonian generated when the cutoff is lowered, and we are interested in a second-order 
term generated by the interactions in Eq. (3.6), 

v= \ I d ^ d & d & d ^ (16tt 3 )5 3 (qi + q 2 + q 3 - q 4 ) 

o J (B.3) 

u 4 (-qi, -q 2 , 94) a\q 1 )a ] (q 2 )a\q 3 )a(q4) , 

=^ J dqi dq 2 dq 3 dq 4 (16n 3 )5 3 (qi - q 2 - q 3 - g 4 ) ^ 

u 4 {-qi, q 2 , <?3, 94) a) (qi)a(q 2 )a(q 3 )a(qA) • 

Since we will find that the incoming and outgoing momenta must be the same, implying 
that e a and e& in Eq. (4.14) are the same in this case, I combine the two second-order 
terms and obtain, 



(16n 3 )5 3 (p' -p) 6v 2 (-p,p) = 

I (7 /I„t f d ~ kld ~ k2 dh e(fa, fa, k3 ) ^ k \ h) { t k2 t l - \p) ; m) 

O. J p K-y 

where p~ = p 1 ~ 2 /p + , etc. Q(k\, k 2 , k 3 ) is the appropriate cutoff, which is determined by 
the transformation employed. One can readily verify that 
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(p'\v*\fa, k 2 , fa) = (16ir 3 )5 3 (p' -fa- fa- fa) u±(-p' ', fa, fa, fa) , 



(5.6) 



(ki, fa, fa\v\p) = (167r 3 )5 3 (p -ki-k 2 - fa) u 4 (-fa, -k 2 , 



h,p) 



(5.7) 



Substituting these results in Eq. (B.5) leads to the final result, 



Sv 2 (-p,p) = — 




dfadfadfa Q(fa, k 2 , fa) (16tt 3 )5 3 (p - fa - k 2 

U4(~P, fa, fa, fa) u 4 (-fa, -fa, ■ 
j) k-^ k^ 



-fa) 
fa,p) 



(5.8) 



To obtain Su 2 from Sv 2 we must complete the set of re-scalings appropriate to the trans- 
formation. If we want to study spectator dependence, we simply need to compute matrix 
elements between states containing additional particles, which directly affect the cutoff 
function © for some transformations. 

Appendix C: Physical Masses in the Light-Front Renormalization Group 

What happens if the physical particles are massive? In this case one expects to find 
at least one mass parameter in the relevant mass operator that is an independent function 
of the cutoff. I assume that this is the mass term that produces the correct relativistic 
dispersion relation for free particles (i.e., the part of u 2 that does not depend on either 
transverse or longitudinal momenta), and I call this the 'physical mass', even though it is a 
running mass in the Hamiltonian that should not be confused with the mass of a physical 
particle. The complete 'mass' operator contains an infinite number of relevant operators 
(i.e., functions of longitudinal momentum fractions that produce a complicated dispersion 
relation) even in the critical theory, and one should expect new functions to appear in 
the massive theory. In the massive theory, one can again use the conditions Wilson and 
I have proposed, selecting a single coupling and a single mass to explicitly run with the 
cutoff, while all other operators depend on the cutoff only through their dependence on this 
coupling and mass. In this Appendix I briefly study a few of the consequences of adding 
a physical mass. This study is both incomplete and preliminary. I explicitly consider 
only the second-order behavior of the transformation, and I focus on the portion of the 
Hamiltonian trajectory near the critical Gaussian fixed point where the physical mass is 
exponentially small in comparison to the cutoff. 
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The physical mass is proportional to some Ajv", which sets the scale at which the 
physical mass can no longer be regarded as small and must be treated nonperturbatively. 
Since = 2~ Ar A Q , as we construct a renormalized trajectory we expect the physical 
mass in H^°, mo, to be exponentially small in comparison to Ao- In other words, on 
a renormalized trajectory we should find an infinite number of Hamiltonians near the 
critical fixed point in an asymptotically free theory. In a scalar theory, if we let Ao — > oo 
while keeping the initial coupling small, we expect to find the trajectory approach the 
critical fixed point and then asymptotically approach the line between the critical and 
trivial Gaussian fixed points, with the strength of the interaction going to zero. For an 
infinitely long trajectory there are interactions only over an infinitesimal initial portion of 
the trajectory, after which the trajectory misses the critical fixed point by an infinitesimal 
amount and then follow infinitesimally close to the line of Gaussian Hamiltonians joining 
the critical and trivial Gaussian fixed points. 

After adding a physical mass to the Hamiltonian, the analysis that led to Eq. (6.11) 
yields 



5v4 ■ 



2 f d 2 S ± dx ™-L 2 ^ m 2 ,±2^_ 2 



r 



16tt 3 x(1 - x) L y(l - y) x(l - x) J 



a 2 , , ^i 2 



V x(l — x) J \x(l — x) J 

Note that I have placed the physical mass in the denominator. If we are considering 
perturbations about the critical Gaussian fixed point, we should include the mass as an 
interaction; but we can find the results of such a treatment here by expanding in powers 
of the physical mass. When the physical mass is an independent parameter in the Hamil- 
tonian, and it is small, we expect to find a well-defined perturbative expansion in powers 
of the running coupling and in powers of the running mass. Here I seek a perturbative 
renormalization group in which such expansions exist. 

In Eq. (C.l) I have used an invariant-mass cutoff that contains no mass, and the 
same problem that plagued T 1 - is encountered; the denominator can change sign inside the 
integral. There is nothing to prevent m 2 /y/(l — y) — m 2 /x/(l — x) from becoming large 
and positive for small y. If we expand the denominator in powers of r 1 - , and complete the 
integrals, the coefficients in the expansion become arbitrarily large, producing divergences 
in subsequent transformations. 

The only simple solution that I have found is to abandon the first invariant-mass 
transformation and use the transformation that employs the cutoffs in Eq. (3.24). Some 
of the peculiarities of this transformation have already been discussed. Here I want to focus 
on the changes in the analysis of the massless theory wrought by the physical mass terms 
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in energy denominators and cutoffs. Perhaps the most important observation to make at 
this point is that the addition of a small mass may lead to small changes in the analysis, 
because all integrals are finite before the mass is added. To see that this is apparently the 
case, return to the correction to in Eq. (C.l), and use the new invariant mass cutoffs 
to get 



5v 4 = 144^ 



d 2 s^dx 



I 2 o 



s- 1 +?rr-i- 1 



167r 3 x(l — x) L y(l — y) x(l — x) 
V u x(l-x)J \x(l-x) V V 1 y(l-y) ) 



(C.2) 



Here I display the cutoff associated with the incoming state, which prevents the energy 
denominator from going through zero. Let me show the analytic result and use it to 
analyze 6114 when there is a physical mass. Defining 



E 



+ m z 



2/(1 — 2/) 



(C.3) 



we obtain, 



Sv 4 = — 



64tt 2 



0(A?-£W In 



1 + 



4m 2 /hi 



1 - 4m 2 /A 2 



1-^1- Am 2 /h 2 Q 1 + ^1- Am 2 /h\ 



+ 



sJl-4m 2 /Eln 



yfl - 4m 2 /A 2 - y/l - 4m 2 /E ^1 - 4m 2 / h\ + y/l - 4m 2 /E 
yjl - 4m 2 1 h 2 + y/l - 4m 2 1 E yjl - 4m 2 /h\ - ^/\-4m 2 /E, 



(C.4) 



Without the cutoff on the external energy, we would encounter negative arguments in the 
logarithm. With the cutoff the argument can still go to zero; however, I have discussed this 
same basic issue following Eq. (6.15), where I argued that it is not necessarily a serious 
problem for a perturbative analysis. 

The marginal part of this correction is obtained by taking the limit r 1 - — > with y and 
m fixed, which leads to 
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5V4 




"1 + ,/l -Am 2 /hi l-,/l-4m 2 /Af 



.l-^l-4m 2 /Ag 1+ -y/l -4m 2 /A 2 . 



+ v 7 ! - 4y(l - ?/) /n 



y / l-4m 2 /A 2 - VI - 42/(1 - 2/) 



(C.5) 



y / l-4m 2 /A 2 + v/1 - 42/(1 - 2/) 



y / l-4m 2 /A 2 + v /l- 42/(1-2/) 



y / l-4m 2 /A 2 - v/1 -42/(1 -2/). 



Unlike the 0((/ 2 ) correction encountered in the theory with no physical mass, this marginal 
operator has a complicated dependence on y. One can continue to compute transforma- 
tions, using this entire marginal operator, or one can make an expansion of this operator 
in powers of m 2 . Such an expansion should converge, as long as m 2 is exponentially small 
in comparison to A 2 ,. After setting Ai = Aq/2, the expansion yields 



The first term in the expansion is the result obtained for the massless theory, while 
the second term must be considered further. There are cutoffs on the remaining states 
that prevent m 2 /y/(l — y) from becoming larger than Ajj/4, so the first correction to Zn(4) 
is at most 3/4. Moreover, the cutoff actually involves r 1 + m 2 , and not just m 2 . Thus, 
no individual term in the expansion of Sv4 given in Eq. (C.5) diverges within the limits 
imposed by the cutoffs, although the entire sum can diverge. I do not repeat the discussion 
following Eq. (6.15). 

As long as the mass remains small, we can make a perturbative expansion of every term 
in powers of the mass divided by the cutoff. In this case, the mass is treated as if it were 
a transverse momentum. This must be done with some care. As we have already seen in 
the critical theory, the step function cutoffs sometimes lead to singular distributions when 
limits involving their arguments are taken. Just as taking transverse momenta to zero 
in the critical theory led to delta functions involving longitudinal momentum fractions, 
taking the mass to zero may produce such singular distributions, and it is essential that 
this limit be taken exactly. Thus, as long as the mass is small, we can continue to expand 
the Hamiltonian about the critical fixed point provided we carefully evaluate distributions. 
We expand every operator in powers of transverse momenta, because the linear analysis 
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reveals that such powers lead to increasingly irrelevant eigenoperators of the linearized 
transformation. While powers of the physical mass do not lead to increasingly irrelevant 
operators, they do lead to increasingly large powers of ?ti 2 /Aq. The subsequent rescaling 
replaces each power with 4 m 2 /Ag; however, the expansion should remain reasonable until 
we run the effective cutoff down to the point where m 2 /A 2 is not small. 

To this point I have ignored spectators in the discussion of 8V4 . Spectators have a 
drastic, potentially disastrous, effect on the analysis. When we studied the marginal part 
of Su^ in the massless theory, spectators had no effect, because their transverse momenta 
were set to zero to find the marginal operator and they dropped out of the cutoffs as a 
result. However, if the spectators are massive, they do not drop out of the cutoffs when 
their transverse momenta are zero. Their presence effectively lowers the cutoff. If there 
are n spectators, the cutoffs are effectively shifted, 

1=1 1=1 

where X{ are the longitudinal momentum fractions of the spectators. If the total longitu- 
dinal momentum fraction of the spectators is 1 — w, the smallest shift that this sum can 
introduce is nm 2 / (1 — w). For sufficiently large n, this always becomes larger than A 2 ,, and 
we find that the invariant-mass cutoff is also a cutoff on particle number. Moreover, in the 
highest sectors of Fock space that survive this cutoff, the effective cutoff is arbitrarily small 
and we must expect that any change in this cutoff produces nonperturbative effects in the 
exact results. I see no way around this conclusion; however, the error made by treating 
the highest sectors of Fock space perturbatively as the cutoff is changed may still be small 
if one is interested only in states of much lower energy than the states removed, which is 
always the case in a renormalization group analysis. 

When the effective cutoff becomes very small, it is because the spectator state is a 
high energy state. Thus, the nonperturbative problem occurs when we need to accurately 
approximate the part of the Hamiltonian that governs states of high energy. In this case, 
the intermediate state has high energy because there are a large number of massive par- 
ticles. If our ultimate interest is to study states with such high energy, we do not expect 
to be able to lower the cutoff to this scale without encountering a nonperturbative prob- 
lem. However, if our ultimate interest is to study low energy states, we need only concern 
ourselves with the error made in the matrix elements that ultimately govern the mixing 
of these many-body high energy intermediate states with the few-body low energy states 
of interest. It may be possible to accurately approximate these matrix elements using 
perturbation theory even when large errors are made in the complete wave function for 
the many-body state. Since it is far easier to study such spectator effects quantitatively 
when computing Su2, I turn to this. 
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In order to compute Su2 for the massive theory, we can follow the same steps that led 
to Eq. (6.25) and simply add a mass to obtain 



■2 r , 72„±', / ,./,/->,. U a,,i a^X'aJ r t -L 2 q ±' 2 r ±' 2 s ±' 2 2m 2 T- ! 



c g l2 f d ? dx' dr r 1 - dy' d 2 s ± dz' 
8v 2 = T7«)A • 



3! u 7 16ttV 16ttV 16ttV LwA 2 x' y' z> wA 2 \ 
(16tt 3 )5(1 -x'-y'- z')8\^L- - q-L' - - s^) 



± 2 , 1*2 o_2 „_L /2 „_L /2 „_L' 2 



. +M 2 3m 2 q 



(1 - w)Aq wA^ x' y' z' 

-i) 



±'2 r± /2 s± /2 t ±2 + M2 3m 

"^ + ^ + ^ + 7^ — + 

T 



y' z' (1 - w)Kl wAl 



+ m 2 t ±2 + M 2 



wA 2 (l-w)Al 

(C.8) 

Here the invariant-mass-squared of the spectators, M 2 , is dependent on the individual 
longitudinal momentum fractions of the spectators, and it no longer vanishes when all 
spectator transverse momenta vanish. I have displayed the cutoff on the external state 
to make it clear that an expansion in powers of m 2 may be reasonable even though one 
ht fear that M 2 , 1/(1 — w), or 1/w could become large. When 5v2 is expanded in 
powers of m 2 , an infinite number of functions of longitudinal momentum fractions appear. 
Since this happens even when all external transverse momenta are taken to zero, this 
means that an infinite number of relevant operators appear. Relevant operators usually 
must be precisely controlled because they grow at an exponential rate. The exception is 
when their initial strength is exponentially small, and remains exponentially small over 
all but the final part of a trajectory of Hamiltonians. In this small error in the 

coefficient of a relevant operator remains small over the entire trajectory. This error 
becomes exponentially large only when the mass becomes exponentially large, and I am 
interested in approximating the portion of the renormalization group trajectory over which 
the mass is small. There is a separate important question of how precisely one must 
approximate this trajectory if one wants accurate predictions for low energy observables 
to result from an exact diagonalization of the final Hamiltonian on the trajectory. I do not 
address this issue here. 

If we expand every term in the transformed Hamiltonian in powers of m 2 , we find that 
the second-order transformation is almost identical to the massless case. We can compute a 
trajectory of massive Hamiltonians, using the second-order transformation and expanding 
every term in powers of m 2 / A 2 ,. We can choose to be 
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d 2 p ± dp+ 

16tt 3 P+ 



7 lD7r J Xl 



P 



_L2 



+ (1+6)- 



+ mr 



+ Mo 



+ 

p_L2 



167r 3 xi 7 167r 3 a;2 



+ (i+e )^^ + (i+eo) 1 2 



(lGTT 3 )^^! + Y L 2 )5{1 - Xi - x 2 ) 

9i»92)(?i,92| + 



+ m. 



<7o 
6 
<7o 
4 



+ 



go 
6 



c% (167t 3 )5 3 (qi + q 2 + qs - q±) (J (qi)a< (q 2 )a< (q s )a(q4) 

dqi dq 2 dq s dq 4 (167r 3 )5 3 (<?i + q 2 - <?3 - Qa) a\qi)a\q 2 )a(qs)a(q4) 
dqi dq 2 dqs dq 4 (167r 3 )5 3 (<?i - q 2 - q s - g 4 ) a\qi)a(q 2 )a(q 3 )a(q4) . 



(C.9) 

This Hamiltonian contains a complete set of relevant and marginal operators through 
0(m 2 1 'Aq). When the second-order invariant-mass transformation is applied to the Hamil- 
tonian, and all irrelevant operators are dropped, the resultant Hamiltonian contains no 
new relevant or marginal operators to these orders, and the only effect is to change the 
constants mo, £o> A*o 5 an d go- Thus we can write the relevant and marginal operators in 
as in Eq. (6.35), using the constants m n , £ n , /z n , and <7 n . 

The complete second-order equations for the evolution of £ n , /u n , and g n are identical 
to Eqs. (6. 36)- (6. 38). Moreover the equation for the evolution of m to this order is trivial, 



m 2 n+1 =4m 2 n . (CIO) 

It is clear that the presence of a physical mass severely complicates the higher-order anal- 
yses, but the qualitative features of the analysis for the massless theory should survive 
when an additional expansion in powers of m 2 / Aq is made. 
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Figure Captions 

1. Wilson's triangle of renormalization. The units are chosen so that Ayv = 1. 

2. Examples of Hamiltonian diagrams, (a) in which a single energy denominator ap- 
pears, and (b) in which two energy denominators appear. The arrows indicate the 
energy differences found in these denominators. 

3. Second-order corrections to: (a) U2, (b) W4, (c) uq, and (d) u$. 

4. Second-order correction to U2 with spectators. 

5. Third-order corrections to the one-boson to three-boson part of 114. 

6. (a) A two-loop correction to u\ paired with the appropriate counterterm insertion in 
a one-loop correction to U4. (b) The source of the counterterm. 

7. (a) A two-loop correction to 114 paired with the appropriate counterterm insertion in 
a one-loop correction to 114. (b) The source of the counterterm. 

8. (a) A two-loop correction to paired with two appropriate counterterm insertions 
in one- loop corrections to 114. The sources of these counterterms are shown in (b) 
and (c). 

9. (a) A two-loop correction to paired with two appropriate counterterm insertions 
in one- loop corrections to u\. The sources of these counterterms are shown in (b) 
and (c). 

10. One-loop corrections to the two-boson to two-boson part of with one marginal 
counterterm vertex. 

11. One-loop corrections to the one-boson to three-boson part of w.4 with one marginal 
counterterm vertex. 
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